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Preface 


This  tutorial  is  aimed  at  the  nonlinear  novice  trying  to  get  a 
handle  on  pho tore tractive  processes.  Often,  the  engineering  or  physics 
published  article  or  text  leaps  from  one  point  to  the  next  with  only  a 
cleverly  placed  "it  can  be  easily  shown  that..."  -type  expression  as  the 
connector.  The  author  can  hardly  be  blamed  since  these  intermediate 
steps  may  seem  redundant  after  being  through  them  a  good  number  of 
times.  He  is  also,  undoubtedly,  working  under  a  publisher's  directive 
for  economy  of  space. 

The  typical  student,  on  the  other  hand,  may  not  view  this  missing 
material  as  nonchalantly.  Nevertheless,  he  often  accepts  the  stated 
points  on  faith  because  the  usual  pressure  of  a  time  constraint  prevents 
nim  from  going  through  the  tedium  of  a  derivation.  Valuable  insights, 
implicit  assumptions  and  boundaries  of  applicability  remain  hidden  along 
with  the  true  understanding  of  the  concept  and  a  "turn  the  crank" 
approach  is  accepted. 

The  purpose  of  this  tutorial  is  to  gather  pertinent  information  from 
articles  and  texts  into  one  place  and  fill  in  some  of  the  blanks  the 
original  authors  left  behind.  It  is  specifically  designed  to  lay  the 
foundation  for  another  Master's  thesis  which  would  pick  up  where  this 
tutorial  leaves  off  and  apply  these  theories  to  observed  phenomena. 

Of  course,  being  a  student,  myself,  I  could  n'ot  have  put  this  paper 
together  alone.  First,  I'd  like  to  sincerely  thank  my  advisor,  Dr.  Ted 
Luke,  for  lending  me  his  expertise  in  the  area  and  for  his  patience  with 
me  as  I  visited  his  office  what  seemed  to  be  several  times  daily.  I'd 
also  like  to  thank  Dr.  Won  Roh  for  valuable  conversations  on  these 
topics.  Finally,  my  unending  gratitude  goes  to  my  wife,  Melony,  who  was 

ii 


an  "AFIT  widow"  for  most  of  this  project  then  became  the  best  typist  I 
could  have  hoped  to  find. 
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Abstract 

This  thesis  presents  a  detailed  interpretation  of  published  theory 
of  nonlinear  optical  phase  conjugation  by  degenerate  four  wave  mixing 
and  by  the  photorefractive  effect.  Photorefractive  phase  conjugation  is 
shown  to  be  a  low  incident  intensity  alternative  to  degenerate  four  wave 
mixing  (which  requires  high  intensity  incident  light  to  achieve  the 
nonlinear  polarization  of  a  medium). 

The  derivations  of  two  models  for  the  photorefractive  effect, 
Feinberg's  Hopping  Model  and  Kukhtarev's  Solid  State  Model,  are 
presented.  The  significance  of  Kukhtarev's  model  is  his  development  of 
criteria  for  which  spatially  sinusoidal  incident  light  (interference 
fringes)  produces  a  spatially  sinusoidal  electro-static  field  in  the 
crystal  -  an  assumption  casually  made  in  much  photorefractive  work.  Both 
models  culminate  in  expressions  for  the  magnitude  of  the  space-charge 
field  and  the  spatial  phase  shift  between  the  field  and  the  incident 
interference  pattern  for  small  modulation  conditions. 

Feinberg's  model  is  extended  to  develop  expressions  for  two-beam 
coupling,  diffraction  efficiency  and  phase  conjugate  reflectivity.  An 
alternative  development  for  these  expressions  assumes  the  electro-static 
field  exists  in  the  crystal  and  modulates  the  crystal's  refractive  index 
by  the  linear  electro-optic  (Pockel's)  effect. 

The  functional  dependencies  of  Feinberg's  expressions  agree  with 
(published)  empirical  data.  However,  specific  crystallographic 
orientations  exist  for  which  the  absolute  magnitude  of  Feinberg's 
expressions  or  even  the  relative  magnitudes  between  different 
orientations  are  inaccurate.  Accurate  relative  magnitude  expressions  are 


made  by  assuming  an  electro-static  field  exists  in  the  crystal  and  using 
an  index-ellipsoid  approach  to  determine  refractive  indices. 
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NONLINEAR  OPTICAL  PRINCIPLES  AND  THE  PHOTOREFRACTIVE  EFFECT  APPLIED  TO 


OPTICAL  PHASE  CONJUGATION 

I.  Introduction 

Optical  phase  conjugation  occurs  when  more  than  one  coherent  light 
beam  is  incident  on  a  nonlinear  material  oriented  properly.  A  new  light 
beam  is  produced  which  is  the  phase  conjugate  of  one  of  the  incident 
beams,  i.e.,  it  propagates  back  in  the  direction  of  the  incident  light, 
and  if  the  incident  light  is  diverging,  the  phase  conjugate  will 
converge  or  vice  versa.  The  phase  conjugate,  then,  appears  to  be  the 
time  reversal  of  the  incident  beam.  Because  phase  conjugation  is  often 
compared  to  reflection  from  a  mirror,  the  nonlinear  material  producing 
this  effect  is  referred  to  as  a  "phase  conjugate  mirror"  or  PCM 
(Figure  1)  (1:157). 


Incident  wavefront 


Reflected  wavefront 


Figure  1.  Illustration  of  phase  conjugate  reflection  as  compared  to 
reflection  from  an  ordinary  mirror  (1:157) 

The  practical  usefulness  of  PGMs  is  evident  in  these  few  examples. 
First,  since  a  PCM  conjugates  the  phase  of  incident  light,  aberrations 
given  to  a  wavefront  as  it  passes  through  a  medi’jm  will  be  removed  if 
the  wave  back-propagates  through  that  same  medium.  Because  of  this, 
glint  from  the  illumination  of  a  target  may  be  amplified,  reflected  from 


Ordinary  mirror  Phase  conjugate  mirror 
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a  PCM  and  passed  back  through  the  amplifier  and  aberrating  medium  to  be 
focused  (with  much  higher  energy)  to  a  diffration-limited  spot  on  the 
target.  Because  the  direction  of  a  phase  conjugate  wave  is  the  reverse 
of  the  incident  wave,  this  is  an  automatic  point  and  track  system 
(Figure  2)  (1:175). 


Figure  2.  PCM  application  of  diffraction-limited  focusing  and  automatic 
point-and- track  (1:175) 


Another  example  is  the  use  of  a  PCM  in  a  laser  cavity.  Here,  the 
phase  conjugate's  "time  reversed"  reflection  removes  the  ordinary  mirror 
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stability  criteria  of  curvature  and  spacing,  and  allows  only  one  more 
highly  monochromatic  mode.  (Figure  3)  (1:177-179). 


Figure  3.  PCM  application  as  one  of  the  mirrors  in  a  laser  cavity 
(1:178) 

The  first  PCMs  were  produced  by  degenerate  four  wave  mixing  (all 
waves  are  of  the  same  frequency),  a  third  order  nonlinear  process.  The 
geometry  for  this  effect  consists  of  two  counterpropagating,  highly 
intense  pump  beams  and  a  probe  beam  which  are  incident  on  and  cause  the 
polarization  of  a  nonlinear  material.  There  are  several  orders  of 
polarization  produced  but  the  third  order  polarization  wave  is  of 
interest  because  it  stimulates  an  electro-magnetic  wave  which  is  the 
conjugate  of  the  probe.  Both  the  conjugate  and  probe  beams  then  gain 
intensity  by  extracting  energy  from  the  pump  beams. 

The  electric  field  which  drives  the  third  order  polarization  in  this 
configuration  is  due  to  the  superposition  of  the  incident  light  waves. 
Since  the  third  order  optical  susceptibility  which  "links"  this  E-field 
to  the  medium  is  very  small  in  magnitude,  very  high  intensity  incident 
light  is  required  to  produce  the  phase  conjugate  effect. 

A  technique  for  producing  phase  conjugation  with  lower  intensity 
incident  light  is  through  the  use  of  photorefractive  crystals.  The 
interference  of  the  (coherent)  probe  and  one  of  the  pump  beams  "writes" 
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a  refractive  index  grating  in  the  crystal.  Then,  the  counterpropagating 
pump  beam  is  diffracted  off  this  grating  and  produces  the  conjugate  of 
the  probe. 

Hie  index  grating  is  written  through  the  photorefractive  effect.  A 
photorefractive  crystal  contains  various  sites  in  its  lattice  which 
accept  and  hold  electrons  with  very  low  thermal  probability  of  freedom. 
Under  illumination,  however,  these  charges  are  freed  and  the  movement  of 
either  electrons  or  holes  (depending  on  the  crystal)  occurs  until  the 
charge  is  retrapped  in  an  area  of  the  crystal  which  is  not  illuminated. 
In  the  case  of  interference  fringes,  then,  interference  maxima 
correspond  to  the  location  the  charges  left  and  minima  to  the  locations 
where  they  are  retrapped. 

The  separation  of  opposite  charges,  either  electrons  or  holes 
(whichever  is  the  carrier)  and  the  ions  left  behind,  causes  an 
electro-static  field  to  form.  This  field  modulates  the  index  of 
refraction  of  the  crystal  through  the  linear  electro-optic  (Pockel's) 
effect.  Since  the  interference  fringes  are  sinusoidal  in  intensity,  the 
charge  distribution,  E-field  and,  thus,  the  index  grating  are  all 
sinusoidal,  as  well. 

For  this  type  of  effect,  the  E-field  driving  the  nonlinearity  of  the 
crystal  is  the  created  space-charge  field.  The  magnitude  of  this  field 
is  proportional  to  the  modulation  index  of  the  interfering  pump  and 
probe  beams.  Since  the  modulation  index  is  dependent  on  the  relative 
intensities  of  these  beams  rather  than  the  incident  intensity  of  either 
of  them,  phase  conjugation  can  be  produced  from  lower  intensity  incident 
light. 

In  order  to  theoretically  describe  light  incident  on  the  index 
grating  that  is  formed,  a  model  describing  the  grating,  itself,  is  first 
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necessary.  Two  models  of  the  E-field  forming  the  grating  are  Feinberg's 
Hopping  Model  and  Kukhtarev's  Solid  State  Model. 

Feinberg's  model  is  based  on  the  idea  that  a  charge  may  occupy  any 
of  a  large  number  of  sites  in  the  photorefractive  crystal  and,  under 
illumination,  "hops"  to  any  of  its  nearest  neighbor  sites.  In  this 
model,  then,  the  probability  of  charge  movement  at  any  site  is 
proportional  to  the  probability  that  the  site  and  its  nearest  neighbors 
are  occupied  and  to  the  illumination  at  the  site  and  at  its  nearest 
neighbors.  The  solution  to  the  model  for  assumed  spatially  sinusiodal 
illumination  yields  spatial  charge  distribution  (and,  therefore, 
potential  distribution  and  electro-static  field  distribution). 

Kuhktarev's  model  is  based  on  more  standard  solid  state  arguments 
such  as  charge  generation  rates  and  transport  lengths.  His  solution  (as 
well  as  Feinberg's)  for  the  space-charge  field,  given  sinusoidal 
illumination  with  small  modulation  index,  is  a  sinusoidal  field  of  the 
same  spatial  frequency  as  the  illumination  but  which  may  be  out  of  phase 
with  it.  This  phase  shift  will  be  shown  to  be  an  important  factor  later. 

It  has  been  observed  experimentally  that  when  two  beams  write  an 
index  grating  in  a  photorefractive  crystal  and  propagate  through  the 
grating,  one  beam  may  experience  a  gain  in  intensity  while  the  other 
experiences  a  corresponding  loss.  Two  parallel  theories  describing  this 
"two-beam  coupling"  are  presented.  One  is  developed  from  Feinberg's 
Hopping  Model.  The  other  assumes  a  sinusoidal  index  grating  exists  and 
can  be  described  by  the  electro-static  field  through  an  "index 
ellipsoid"-type  argument.  Both  conclude  the  key  factor  in  two-beam 
coupling  is  the  phase  shift  between  the  electro-static  field  and  the 
incident  interference  fringes.  Without  such  a  phase  shift,  no  energy 
transfer  could  occur. 
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Two-beam  coupling  can  be  considered  a  figure  of  merit  for  the 
photorefractive  effect.  Other  such  figures  of  merit  are  diffraction 
efficiency  and  phase  conjugate  reflectivity.  Each  of  these  figures  of 
merit  is  important  because  they  allow  quantification  of  those  aspects  of 
the  photorefractive  effect  which  play  a  significant  role  in  optical 
phase  conjugation. 

Diffraction  efficiency  is  the  amount  of  power  in  a  wave  diffracted 
at  the  Bragg  angle  off  a  grating  compared  to  the  amount  of  power  in  the 
wave  incident  on  the  grating.  Phase  conjugate  reflectivity,  in  the 
photorefractive  phase  conjugation  geometry  described  earlier,  is  the 
power  in  the  conjugate  wave  diffracted  off  the  grating  compared  to  the 
power  in  the  incident  probe  wave  which  had  a  part  in  writing  the 
grating.  Again,  parallel  theories  from  the  extension  of  Feinberg's  model 
and  the  sinusoidal  grating  assunption  are  presented.  The  mutual 
conclusion  here  is  the  magnitude  of  the  (space-charge  field  and, 
therefore,  the)  index  grating  is  proportional  to  the  diffraction 
efficiency  and  also  to  the  phase  conjugate  reflectivity. 


II.  Second  and  Third  Order  Nonlinear  Optical  Principles  -  Background  to 

Optical  Phase  Conjugation 

Two  types  of  optical  phase  conjugation  are  presented  in  this  paper  - 
degenerate  four  wave  mixing  and  photorefractive.  The  first,  degenerate 
four  wave  mixing,  occurs  because  the  superposition  of  light  waves  in  a 
nonlinear  medium  causes  the  nonlinear  polarization  of  that  medium.  This 
polarization  gives  rise  to  an  electro-magnetic  wave  as  described  by 
Maxwell's  equations.  In  the  proper  configuration,  this  E-M  wave  can  be 
the  phase  conjugate  of  one  of  the  input  waves. 

This  chapter  presents  a  general  development  of  the  nonlinear 
processes  mentioned  above  with  the  following  chapter  describing  the 
special  configuration  for  phase  conjugation.  Although  four  wave  mixing 
is  a  third  order  nonlinear  process,  insight  to  nonlinear  processes,  in 
general,  can  be  gained  through  the  more  simple  case  of  the  second  order. 


Second  Order  Nonlinear  Principles 

Second  Order  Nonlinear  Polarization 

When  electric  fields  superimpose  in  a  medium,  the  linear  and 

nonlinear  polarization  of  that  medium  may  be  written  in  the  form: 

P .  =  e  rX  •  y  •  •  E  •  +2i-z,  d  •  • ,  E.  E,  +  ... 
i  xiJ  J  J  k  ljk  j  k 

where 


=  i^1  component  of  polarization  (i  =  x,  y  and  z) 

E^  =  1  n  component  of  the  instantaneous  field 
X—  =  Linear  susceptibility 
d^k  =  Second  order  nonlinear  susceptibility 
and  the  susceptibilities  are  invariant  under  shuffling  of  subscripts, 

e-8‘  Xij  =  Xji  and  dijk  =  dikj’  etc*  (2:504)- 

The  linear  polarization  from  equation  (1)  is: 
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(2) 


For  example,  the  component  of  linear  polarization  in  the  x  direction 
would  be: 

P  ^  =  eA  (  y  E+v  E+y  E)  (3) 

x  0  Axx  x  Axy  y  Axz  z 

The  second  order  polarization  is: 

h™- '  “jS  dijk  Ej  h 

for  degenerate  cases  or,  in  general,  for  cases  with  E-fields  of  more 
than  one  frequency  present,  it  is: 

PiNL<t>  ’  2£j  Wm  dijk  (4> 

For  example,  the  x-component  would  take  the  form: 

P>)  =  2  {dxxx  [Exa(t)  +  Exu2(t)  +  ...  +  E^Ct)]2 

+  2  dxxy  [Exa(t)  +  Exuj2( t)  +  ...  +  Ex^(t)] 

[Ey^(t)  +  Ey^U)  +  ...  +  E^t)] 

+  2  dxxz  [Ex^(t)  +  Ex^(t)  +  ...  +  E^Ct)] 

[EzL)1(t)  +  Ezuj2(t)  +  ...  +  EzlJ1(t)] 

+  d  [E  ^(t)  +  E  ^(t)  +...  +  EvlJ1(t)]2 

+  2  dxyz  [Ey^(t)  +  Ey^Ct)  +  ...  +  EylJ1(t)] 

[Ezu)1(t)  +  EzLj2(t)  +  ...  +  Ezuri(  t)  ] 

+  dxzz  [Ezwl(t)  +  Ez(j2(t)  +  ...  +  E^U)]2}  (5) 

Expansion  of  the  E.^(t)  E,  ^(t)  terms  gives: 

J  * 

EVd)  Ek‘jn(t)  =  (1/2)  [A.^  expdu.t'  +  mV5*  expC-i^t)] 
(1/2)  [A^^  exp(i^t)  +  Ak('lJT1')"  exp(-i^t)] 

=  (l/4){[Aj(wl)  A^1^  exp[i(Ul  +  u^)tl 
+  Aj<'Ul''  4k('(JT1'>'‘  exp[i(Wl  -  u^)t] 

+  A.(^>VAk(jTl)  exp[-i(Wl  -  ^t] 

+  a/^"  A^^”  exp[-i((j1  +  ) t ] }  (6) 

Thus,  this  combination  produces  frequency  te^s  in  the  nonlinear 
polarization  which  are  sum  and  difference  frequencies  of  the  driving 
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E-fields.  The  frequency  terms  in  the  previous  example  would  be: 

P  ^(t)  =  P  ^(t)  +  P  ^(t)  +  P  ul  +  w2(t)  +  P  ul  “  u2  (t)  +  ... 

X  X  X  X  X 

+  p^1  +  ^(t)  +  Px  wn  (t)  +  Px2<j2(t)  +  ... 

+  Pxw2  +  un(t)  +  pjw2  "  unl(t)  +  ...  +  P  ^"(t)  (7) 

In  general,  the  sum  frequency  terms  are  in  the  form: 
p_ul  +  ^2(t)  =  exp[i((l)1  +  U)2)t]  IjZk  dijk  A^1^  A^2'* 

+  exp  [-i(Ul  +  u)2)t]  Ejj;k  d^k  a/u1^‘  Ar^2^  (g) 

where  oj.  and  are  different  frequencies.  The  form  of  the  difference 
terns  is: 

p_|(jl  -(j2|(t)  =  -  W2)t]  Zjlk  A^1^  Ak^2^“ 

+  exp  [-iCujj  -  u,2)t]  Zjik  Ak^2'>  (9) 

Again,  cjj  and  u ^  are  different.  The  general  expression  for  the  second 
harmonic  case  *  u2  =  w)  is: 

P>(t)  =  (1/2)  [exp(i^t)  ^  dijR  A.^  Ak^} 

+  exp(-i2jt)  Zj  L^.  dijk  A^^'  A^-*"]  (10) 

The  factor  of  2  difference  between  expressions  (8)  and  (10)  follows 
directly  from  the  expansion  of  equation  (5)  (2:506,507). 


-9- 


jk  where  xx  =  1,  yy  =  2,  zz  =  3,  yz  =  zy  =  4,  xz  =  zx  =  5  and  xy  =  yx  = 

6  (2:507).  Certain  elements  of  the  second  order  nonlinear  optical 
susceptibility  tensor  are  always  zero  for  specific  crystal  point  groups. 
(Refer  to  Reference  2  (pages  508-511)  for  further  details.) 

Summary 

An  electro-magnetic  wave  may  be  defined  by  a  traveling  electric 
field  oscillating  at  some  frequency,  When  a  number  of  such  E-fields 
superimpose  in  some  medium,  they  may  be  mixed  via  the  nonlinear 
susceptibility  of  the  medium  to  produce  a  polarization  of  the  medium. 

For  second  order  polarization,  no  matter  how  many  fields  are  present, 
simultaneously,  the  frequency  components  of  the  polarization  can  be 
defined  by  taking  the  fields  two  at  a  time.  Such  a  component  is  then 
either  the  sum  or  difference  of  the  frequencies  of  the  two  fields 
considered.  The  amplitude  associated  with  each  frequency  component  is 
given  by  the  preceding  equations  and  is  dependent  on  the  second  order 
nonlinear  optical  susceptibilities,  d^^* 

Nonlinear  Polarization  in  a  Medium  Gives  Rise  to  an  Electro-Magnetic 

Wave 

The  nonlinear  polarization  produced  by  the  interaction  of  light 
waves  through  the  nonlinear  susceptibility  of  a  medium  will  be  shown 
through  Maxwell's  equations,  to  be  the  stimulus  which  generates  a  new 
light  wave. 

The  Wave  Equation 

The  first  of  Maxwell's  equations  to  be  used  is 

V  x  H  =  J  +  a  D/at  (12) 

where 


H  *  Magnetic  field  vector 


J  =  Electric  current  density 


o  =  Conductivity 

E  =  Electric  field  vector 

D  =  Electric  displacement  vector 


cn  E  +  P 


=  Permittivity  of  free  space 
P  =  Electric  polarization 

=  e0  *L  E  +  PNL 
=  Linear  susceptibility 

=  Nonlinear  polarization 


then: 


where 


V  x  H  =  aE  +  e  3E/3t  +  3P^/3t 


e  =  Dielectric  tensor 
=  E0  (1  +  *L} 

The  second  Maxwell’s  equation  used  is: 


V  x  E  =  -3/3t  (UqH) 


where 


Uq  =  Permeability  of  free  space 
If  the  curl  of  both  sides  of  equation  (14)  is  taken: 

V  x  V  x  E  =  -3/3t  (mqV  x  H) 
using  the  vector  identity: 

VxVxE=V-V-E-V2E 

with  the  assumption  that  the  E-fields  are  perpendicular  to  their 
propagation  (k)  vectors  (3:8)  so: 

V  •  V  =  o 

and  equation  (14),  equation  (13)  becomes: 
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(15) 


v2e  =  Mo  o3Eyat  +  Mq  c  aVat2  +  UQ  a2PNL/at2 

th 

If  the  1  component  (i  =  x,  y,  or  z)  of  E  at  frequency,  cj3,  is 
written  as: 

E/^r.t)  =  (1/2)  [A>*3(r)  exp  [i^t  -  kyr)] 

+  Aioj3"(r)  exp[-i(U)3t  -  kyr)]]  (16) 

then: 

32E.iw3(r,t)/ax2  =  (1/2)  {[a2Aiw3(r)/sx2]  expfiC^t  -  k^r)] 

+  2  [aAilj3(r)/ax]  (-ik3x)  exp[i(co3t  -  k3*r)] 

+  Aiu3(r)  (~k3x2)  exp[iu3t  -  kyr)] 

+  [a2Aiu3"(r)/ax2]  exp[-i(W3t  -  k3*r)] 

+  2  [aAiw3“(r)/ax]  (ik3x)  exp[-i(cj3t  -  k3»r)] 

+  Aia)3'(r)  (~k3x2)  exp[-i(W3t  -  k-^r)]}  (17) 

2  2  2  2 
with  similar  expressions  for  a  /ay  and  a  /az  . 

When  the  complex  envelope  varies  little  over  one  optical  period,  a 

"Slowly  Varying  Envelope  Approximation"  (SVEA)  can  be  made: 

a2Aitj3(r)/ax2  <<a  Aiu,3(r)/ax  (18) 

2  2  2  2 

likewise  for  a  /ay  and  a  /az  .  This  approximation  is  necessary  to  make 
the  math  manageable  as  the  development  continues. 

The  Laplacian  of  E/^  (r,t)  now  becomes: 

V2Eiw3(r,t)  =  { -i[k3x  3Aiuj3(r)/ax  +  k3y  aAiu)3(r)/ay 
+  k3z  aAiw3(r)/az] 

-  (1/2 ) |  k3 [  2  Aiu,3(r)}  exp[i(W3t  -  k3-r)J 
+  {i[k3x  aAiu3"(r)/ax  +  k3y  aAiw3“(r)/ay 
+  k3z  3Aiu3',(r)/az] 

-  ( 1/2 > | k3 1  2  A/^'Xr)}  exp[-i(oj3t  -  kyr)]  (19) 
The  first  and  second  time  derivatives  of  E.w  (r,t)  are  in  the  forms: 

3Ei(j3(r,t)/at  =  iu3  EiU)3(r,t) 

and: 
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(20) 


32Eiuj3(r,t)/3t2  =  ^32  Eiu3(r,t) 

The  wave  equation  now  becomes: 

{ -i[ k3x  9Aiu3(r)/ax  +  k3y  3A./3(r)/3y  +  k3z  3Aiw3(r)/3z] 

-  (1/2)  |  k3|  2  AiLj3(r)}  exp[i(u3t  -  kyr)] 

+  { i[k3x  aAiU)3“(r)/3x  +  k3y  aA/3"(r)/ay  +  k3z  3A/3''(r)/3z] 

-  (1/2)  |  k3 1  2  A-w3''(r))  exp  [-i(u3t  -  kyr)] 

=  (i  mq  o  d3/2)  (Aitj3(r)  exp[i(u3t  -  kyr)] 

-  AiU)3'(r)  exp[-i(cj3t  -  k3*r)]} 

-  (u0  e  oj32/2)  {Aiw3(r)  exp[i(u3t  -  kyr)] 

+  A>’3"(r)  exp[-i(oj3t  -  kyr)]} 

+  uo  32p  itj3(r,t)/at2  (21) 

Difference  Frequency  Components 
Consider  the  component  of  nonlinear  polarization  which  is  the 
difference  of  two  interacting  light  frequencies.  With  the  "conservation 
of  energy"  condition,  u3  =  tpie  nor|linear  polarization  term 

takes  the  form: 

=  ul  w2(r>t) 

=  expii[(u»1  -  c^)t  -  (kj  -  k2)-r]}ljCk  dijk  Ajujl(r)Ak‘Jj2*(r) 

+  exp{-i[(u1  -  ut^t  -  (kj  -  k2)*r])ljlk  d^k  A^w1' (r)  Aku2(r) 

sV^3  =  ul  -  u2(r>t)/3t2  .  .(^  .  ug)2?."3  =  ul  "  ^(r.t)  (22) 

Equating  terms  synchronous  with  exp(iw3t)  from  either  side  of  the 
wave  equation  (21): 

{-i[ k3x 9  \iw3(r)/3x  +  k3y  ^^(rVdy  +  k3z  ^iu3(r)/3z] 

-  (1/2)  |k3  |2  Aiu3(r)}  exp[-i(k3*r)] 

=  (iu  o°  ‘V2)  exp[-i(k3*r)] 

-  (i^e  <^2/2)  Aiu3(r)  exp[-i(ky  r)] 

-  u^2  d^k  A^ul(r)Aku2“(r)  expL-id^  -  k^)  -r]  (23) 

2  2 

Since  |k3  |  =  ^  e,  then: 
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k3x  ^^(rVax  +  k3y  3Aitj3(r)/3y  +  k3z  dA^3(r)/dz 
=  (-p0o  u/2)  A.w3(r) 

-  i  w32  u0IjSk  dijk  Ajljl(r)  Aku2"(r)  expOKt^  -  -  k3)-r]  (24) 

Similarly,  equating  terms  synchronous  with  exp(-i  3t)  gives: 
k3x  3Aiu)3'>(r)/3x  +  k3y  3Aitj3“(r)/3y  +  k3z  3Ai°j3"(r)/3z 
=  (“Pq0  w3^3^  A^u3  (r) 

+  i  u0  w32  ZjZk  dijk  A^'Xr)  Aku2(r)  expfiO^  -  ^  -  k3)-r]  (25) 

Sum  Frequency  Components 

Consider  the  nonlinear  polarization  component  which  is  the  sum  of 
two  interacting  frequencies  where  these  frequencies  are  different,  cj3  ¥ 
u-y-  The  conservation  of  energy  condition  is  =  +  w2  an3  tkie 

polarization  is  in  the  form: 

Piu3m<~'1  +w2(r,t) 

=  exp;  i[ (u^  +  w2)t  -  (kj  +  k^-r]}!^  d^  AjUl(r)Aku2(r) 

*-  exp  (-i[(w1  +  w2)t  -  (kj  +  l^)*  1.]} Z k  d^k  AjWl"(r)Akw2"(r) 

3 2?.^ 3  =  ul  +  uj2(r,t)/3t2  =  -(uj  +w2)2  P^3  =ujl  +w2(r,t)  (26) 

Equating  synchronous  terms  from  either  side  of  equation  (21)  yields,  for 
exp(ic3t): 

k3x  ^A|U)3(r)/3x  +  k3y  dA^w3(r)/3y  +  k3z  3APj2(r)/3z 

=  (Uq  0u) j/2)  AiUi3(r) 

-  i  w32  Po^-j^k  di j k  AjLal(r)Akw2(r)  exp[-i(k1  +  -  k3)-r]  (27) 

and  for  exp(-iw3t): 

k3x  3AiU3' (r)/3x  ‘  k3y  3AiLi3'“(r)/3y  +  k3z  3Ailj3"(r)/3z 

=  (-MqO  u>3/2)  Ai‘jl3"(r) 

+  i  p0  w32  jZk  dijk  Ajljl“(r)Akw2"(r)  exp[i(kj  +  -  k3)-r]  (28) 

Second  Harmonic  Case 

Finally,  consider  the  component  of  nonlinear  polarization  which  is 
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the  sum  of  two  interacting  frequencies  and  these  frequencies  are  the 
same.  The  conservation  of  energy  condition  is  =  2oj  and  the  form  of 
the  nonlinear  polarization  is: 

PiUi3  ~  2u’(r,t)  =  (1/2)  |  [exp[  i(2ut  -  2k-r)]ljZk  d.^k  A^( r )  Aku(r) 

+  exp[-i(2ut  -  2k-r)]ljlk  d^  AjU)"(r)AkL)''(r)) 

32p^3  =  2cj(r}t)/at2  =  _^2  p.<j3  *  2u(rjt)  (29) 

Again,  equating  synchronous  terms  yields,  for  exp(i2(jjt): 
k3x  3Ai2u)(r)/3x  +  k3y  3Ai2u(r)/ay  +  k^z  aAp^rVgz 

=  -MQOoAp^Cr) 

-  i  ufl^  d^k  AjH(r)AkH(r)  exp[-i(2k  -  k^).  r]  (30) 
and  for  exp(-i  2uJt): 

k3x  3Ai2j*(r)/Sx  +  k3  3A>"(r)/3y  +  k32  3A>*(r)/az 

=  ~MqO  ojAp^^'(r) 

+  i  2u2  u0Zj\  dijk  Ai^r)  Ak£J>,('r')  exPfi(2k  "  ^3)  *rI  (31.) 

Surmary 

The  nonlinear  polarization  wave  produced  by  the  mixing  of  electric 
fields  in  a  nonlinear  medium  gives  rise  to  another  E-field  (and  thus  a 
new  light  wave).  The  amplitude  of  this  field  can  be  described  by 
solution  to  Maxwell's  wave  equation  and,  under  conservation  of  energy, 
the  frequency  of  the  response  is  the  same  as  the  frequency  of  the 
stimulus.  Therefore,  the  frequency  of  this  new  light  wave  is  the  same  as 
the  frequency  of  the  polarization. 

Phase  Matching's  Role  in  the  Selection  of  the  Wave  of  Highest  Energy 

Transfer  (A  Nondepleted  Input  ExampleT 

The  amplitude  of  the  light  wave  produced  via  the  nonlinear 
polarization  of  a  medium  was  derived  in  the  previous  section.  The 
propagation  vector  of  this  wave  will  now  be  derived. 

Consider  the  sum-frequency  conservation  of  energy  condition, 
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r 

* 


► 


¥ 


(^3  =  ^1  +  u>2 .  The  amplitude  relationships  are  given  by  equations  (27) 
and  (28). 


First  assumption:  the  medium  is  assumed  to  be  transparent  at  u3, 

then  o=0  and  equation  (27)  reduces  to: 

k3x  3AiW3(r)/3x  +  k3y  3Aito3(r)/3y  +  k^z  8Ai‘j3(r)/3z 

=  -i  U32M0XjIk  AjUl(r)Aktj2(r)  exp  (-iAk*r)  (32) 

where  Ak  =  +  1^  -  k3  and  with  similar  expressions  for  equation  (28). 

Second  assumption:  The  depletion  of  the  mixing  beams  (those  at  oj^ 

and  103)  is  assumed  negligible: 

3Aul(r)/3x  =  3A,Lj2(r)/3x  =  0  (33) 

J  K 

with  similar  expressions  for  x,  y  and  the  conjugates.  Then: 

A  .ul(r)  =  A  ul 
J  J 

A^r)  -  Aj"1" 

Ak“2(r)  -  Ak“2 

Ak“2*(r)  -  Ak"2"  (34) 

Equation  (59)  now  becomes: 

k3x  3AiW3<r)/3x  +  k3y  3AiUj3(r)/3y  +  k3z  3AiU)3(r)/3z 

=  -iw32  M0XjXk  dijk  A^1  Aku2  exp(-iAk-r)  (35) 

with  a  similar  expression  for  the  conjugate. 

If  the  solution  to  equation  (35)  is  assumed  to  be  in  the  form: 

A^u3(r)  =  K  exp(-iAk*r)  +  C  (36) 


then: 

3A^u3(r)/3x  =  -iK  Ak^  exp(-iAk*r)  (37) 

with  similar  expressions  for  y  and  z. 

Equation  (35)  becomes: 

-i  K  (Ak*k3)  exp(-iAk*r)  =  -i  w32  UQ^j^k  d^jk  AjUl  Aku2  exp(-iAk* r)(38) 
so  the  constant,  K,  is  in  the  form: 
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(39) 


2  _  ,  .  oil  .  u>2 

CjO  q  p  Z »  d .  .  *  A  •  Ai 

K  =  3___0_j_k__ijk _ j _ k _ 


If  there  is  no  interaction  length  within  the  medium,  no  nonlinear 
polarization  is  produced  to  drive  a  new  E-M  field.  Therefore,  the 
amplitude  of  the  new  wave  is  zero.  Thus,  the  boundary  condition: 


A^to.O.O)  =  0  =  K  +  C 


C  =  -K 


The  solution  to  equation  (35)  is: 


A/tr)  =  K[exp(-i  Ak*r)  -  1] 
By  similar  development,  the  conjugate  solution  is: 

A^u3-'(r)  _  [eXp(i  Ak*r)  -  1] 


where 


2  _  ,  .  wl*  .  lj2* 

u0  uAZ  .1,  d.  A-  A, 

K"  _  3  0  j  k  ljk  j  k 


From  Poynting’s  theorem  (2:13),  the  i1"  mode  power/unit  area  at  u^, 
at  some  interaction  length,  k^L  (in  the  direction  of  the  wave,  k^),  is 
given  by: 

P^3/ area  =  (l/2)(e/uQ)1/2  A.^k^L,  k3yL,  k^L) 

\“3*<k3xL’  k3yL-  k3zL> 

.  n  exp(-iAk*k:!  L/2)  -  exp(iAk*ko  L/2) 

=  (l/2)(e/pn)1/2  K[ - - - - . . . ] 

exp(-i  Ak*k3  L/2) 

j  exp(i  Ak*k^  L/2)  -  exp(-i  Ak*k^  L/2) 

K  [ . . . . . . ] 

exp(-i  Ak*k3  L/2) 

-  (L  /2)  (c/u0)  u3  pQ  (j  R  di-k  Aj  Ak  ) 

(j  k  dijk  AjWl*  tsin2<Ak*k3  L/2)/(Ak*k3  L/2)2]  (44) 

The  power  in  the  produced  wave  (at  u3)  is  maximized  when  the 
argument  of  the  (sin  x)/x  term  in  equation  (44)  approaches  zero.  This 
implies  Ak* k3  =  0  and,  therefore,  that: 
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Ak  =  0 


(45) 


k 


w  y 


t 


This  is  the  "conservation  of  momentum"  Condition- 

Similar  arguments  can  be  made  for  the  frequency-difference  and 
second  harmonic  cases  with  similar  results:  Ak  =  0  for  significant  power 
transfer. 

Summary 

The  light  wave  produced  from  the  second  order  polarization  is  of  a 
frequency  which  results  from  the  mixing  (either  sum  or  difference)  of 
two  input  optical  frequencies.  There  may  be  an  infinite  number  of 
input  frequencies  available  (either  by  the  overlapping  of  a  large  number 
of  beams  or  the  generation  of  additional  frequencies  via  the  mixing 
process)  but  a  specific  frequency,  is  selected  by  meeting  the 
conservation  of  energy  requirement  (m^  =  )  and  the 

conservation  of  momentum  requirement  (Ak  =  0)  so  the  only  significant 
energy  transfer  occurs  to  the  newly  created  beam  at  and  k^. 

Third  Order  Nonlinear  Principles 

The  concepts  of  nonlinear  polarization,  Maxwell's  equations  and 
phase  matching  in  the  second  order  will  now  be  extended  to  the  third 
order.  Phase  conjugation  by  degenerate  four  wave  mixing  is  then  a 
special  case  of  these  third  order  concepts. 

Third  Order  Nonlinear  Polarization 

Again,  when  E-fields  superimpose  in  a  nonlinear  medium,  the 
polarization  of  that  medium  may  be  written  as  in  equation  (1): 

Pi  *  *  •  5  Ej  +  2  ]  l  'iW  Ej 

The  third  order  nonlinear  polarization  for  the  degenerate  case  is: 

piNL  ■  4  '  “  *ijki  Ei  %  Ei  <46> 

or  for  cases  when  more  than  one  frequency  E-fields  are  present: 
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(47) 


p  NL  =  4  I  I  I  Z  Z  I 
ri  j  k  1  m  p  q 


Y  £  ^  F  p 

xijkl  j  \  H 


(jq 


Expansion  of  the  E^0111  E^P  Ep0^  term  gives: 

E^Ct)  EkuP(t)  E^Ct)  =  (l/8){AjtJn  AkuP  A^  exp[i(um  +  up  +  uq)t] 


+  A.  A,  WP  A1aK’"  exp[i(u  +  u  -  u  )t] 
j  k  1  FLVm  p  q 

+  A A,UP  A,0^  exp[i(u  -  u  +  u  )t] 
j  k  1  KL  v  m  p  q 

+  A  um  .up* 
l  k 


+  A. 
J 


unr-' 


A 

Ak 


uq* 

T 

exp[i(um  - 

UP  ” 

Vt] 

^  expt-i(wm  + 

UP  " 

A  uq 
A1 

exp[-i(co  +  co 
rL  v  m  p 

-  “q)t] 

uq'* 

1 

exp[-i(um  ■ 

-  U) 

p 

*  “q)tl 

&* 

3 

< 

exp[-i(um 

+  GO 

P 

+  Uq)t 

(48) 

Thus,  the  frequency  components  of  the  third  order  nonlinear  polarization 
are: 

P  NL/t\  =  n  I  p  (JH  +  up  +  w<^(t)  +  P  !um  +  UP  “  uql(t)  (49) 
i  v  '  m  p  q  i  i 

The  expansion  of  the  first  component  yields: 

P.um  +  up  +  wq(t)  =  K  1  f  Z.  x.  ...  [A.^A^PA,^  exp[i(u  +  u  +  u  )t] 
l  N  j  !•.  1  *ijkl  L  j  k  1  m  p  q 

+  A.^'A^P  A,wq  exp[-i(u  +  u  +  u>  )t] 
j  k  1  KL  m  p  q 

(50) 

where 


K  =  3 


if 


U  ,  Cl  _  _ 

nr  p  q 


and  u  are  different 


-  3/2  if  two  of  u  ,  u  ,  or  u  are  same,  one  different 
m’  p’  q  ’ 

=  1/2  if  u  ,  u  ,  and  u  are  all  the  same 
nr  p’  q 

the  values  of  K  result  from  an  expansion  of  equation  (47)  similar  to 
that  for  the  sum-frequency  and  second  harmonic  cases  in  the  previous 

chapter.  The  amplitudes  are: 

p  (urn  +  up  +  uq)  .  2|/  H  Z  A  urn  A  up  A  uq 

i  j  k  1  xijkl  j  k  1 

p  (urn  +  up  +  uq)*  _  2k  ^  ^  unT'>'  A  up"  A  uq*"  (51) 

i  j  k  1  xijkl  j  k  1  v  J 

The  other  components  yield: 
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p  um  +  up  -  uq  /  v 
ri  '  ' 

■  *  ]i  i  xijki  [*r  vp  v"*  + -p  -  vt] 

+  AJ“"r"  Ak“P"  Aj“q  expt-iLm  +  Up  -  Uq)t] 

(52) 

where  the  amplitudes  are: 

p  [ com  +  up  —  uq|  — .  ok  ^  ^  ^  v  A  a  a 

ri  j  k  1  xijkl  j  k  1 

p  jum  +  up  -  uq[*  =  2R  I  Z  Z  A  u™*  A  up*  A  uq  (53) 

1  J  k  1  *ijkl  j  k  1  v 

Just  as  with  the  second  order  nonlinear  polarization,  any  number  of 
electric  fields  may  superimpose  in  a  medium  and  their  frequencies  may  be 
mixed  via  the  nonlinear  susceptibility  of  the  medium.  For  third  order 
nonlinear  polarization,  however,  the  fields  are  considered  three  at  a 
time  to  produce  a  polarization  frequency  component  which  is  some 
combination  of  sums  and  differences  of  those  three  frequencies.  The 
amplitude  of  such  a  component  is  given  by  the  preceding  equations  and 
dependent  on  the  third  nonlinear  susceptibility,  X-jj^l’ 


The  Wave  Equation 

Just  as  in  the  second  order  case,  the  nonlinear  polarization  will  be 
a  stimulus  for  which  the  response  is  an  electromagnetic  wave.  The  wave 
equation  developed  from  Maxwell's  equations,  equation  (21),  is  in  the 
following  form: 

1  -i[k3x  9Aiw'3(r)/3x  +  k3y  dA^(r)/dy  +  k^z  aA/^rVaz] 

-  (1/2)|  k3|  2  Aiu2(r)}exp  [i(u3t  -  ky  r)] 

+  { i[k3x  3Aiuj3*(r)/ax  +  k3y  3Aiaj3"’(r)/ay  +  k3z  a'Aiu3*(r)/az] 

-  (1/2) |  k3|  2  A.w3'*(r)}  exp[-i(u3t  -  kyr)] 

=  (i  u0  ou32/2)  (Aiuj3(r)  exp  [i(u3t  -  kyr)]  -  A^'Xr)  exp[-i(u3t  -  k^r)]} 
-  (uQ  eu32/2)  (Aiti)3(r)  exp[i(u3t  -ky  r)  +  AiU)3'>(r)  exp[-i(u3t  -  ky  r) ]} 

+  u0  a2  PiW3(r,t)/9t2 

The  third  order  nonlinear  polarization  with  conservation  of  energy 
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k 


condition,  ^3  =  +  ^2  ~  u4’  ^s: 

pi  (T’C)  3  j  k  ixijki  'Aj  (r)  Ak  (r)  Ai  <r) 

exp{i[(wj  +  ^2  -  -  (k^  +  -  k^)*r]} 

+  ?k.  (r)  Ak  (r)  Ax  (r) 

exp  { -i[(w^  +  l>2  ~  “(kj  +  k2  -  k^)'  r]l } 


then: 


*(r,  t)/3  t 

_  /  ,  x2  n  u3  =  col  +  o>2  -  u>4/  .  N 

=  -(0^  +  (J2  -  u4)  Pi  (r, t) 


Equating  terms  synchronous  with  expCiw^t)  from  either  side  of  the 
wave  equation  gives: 

{ _i[ k3x  3Aiu3(r)/3x  +  k3y  3Aiu3(r)/3y  +  3Aiw3(r)/3z] 

-  (1/2)| k^j  2  Aiu3(r)}  exp(-ik3-r) 

=  (i  UqOw3/2  -  Uq  e  m32/2)  Aiu3(r)  exp(-ik3’r) 

"  3  U0U3  j  k  1  *ijkl  Aj  (r)  Ak  (r)  A1  (r) 

exp[ -i(kj  +  k2  -  k4)-r]  (56) 

2  2 

Since  |k3|  =  u3  Uq  e ,  then: 

k3x  ^/^(rVSx  +  k3y  3Aiw3y(r)/3y  +  k3z 
=  (-Mq  ou3/2)  Aiu3(r) 

■7  2  Z  Z  Z  .  o>l,  x  x  .  m4‘V  x 

-  1  3  u0  u3  j  k  1  xijkl  Aj  (r)  Ak  (r)  A1  (r) 

exp[-i(k^  +  k2  -  k^  -  k3)*r]  (57) 

Similarly  equating  terms  synchronous  with  exp(-iu>3t)  gives: 

k3x  9AiU3"(r)/9x  +  k3y  3Aiu3'‘(r)/3y  +  k3z  3Aiw3"(r)/3z 
-  (-  UQ  om3/2)  A.u3*(r) 

+  1  3  Uq  u3  jkl  Xijkl  AJ  (r)Ak  (r)Aj  (r) 

exp[i(kJ  +  kj  -  k4  -  k3)'r]  (58) 

Similar  developments  can  be  made  for  other  sum  and  difference 
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frequency  combinations. 


Sirmary 

Just  as  in  the  second  order,  the  third  order  nonlinear  polarization 
wave  produced  by  the  mixing  of  electric  fields  in  a  nonlinear  medium 
gives  rise  to  a  new  E-field  (and,  thus,  a  new  light  wave).  The 
amplitude  of  this  new  field  can  be  determined  using  the  wave  equation 
and,  under  conservation  of  energy  conditions,  the  frequency  of  the  light 
wave  is  the  same  as  that  of  the  polarization  wave. 

Phase  Matching 

Again,  just  as  in  the  second  order  development,  a  condition  for  high 
energy  transfer  from  the  input  waves  to  the  produced  wave  will  be 
developed.  The  conservation  of  energy  condition,  u>3  = 
will  again  be  considered.  Equations  (57)  and  (58)  give  the  amplitude 
relationships. 

The  two  assumptions  made  in  the  second  order  development  are  again 
made:  the  medium  is  transparent  at  therefore,  o=0;  and  the  input 
beams  are  not  significantly  depleted  by  the  mixing  process  so: 

9Ajtil'*‘(r)/3x  =  3A^u^(r)/3x  =  3A^tjZ*(r)/3x  =  0  (59) 

with  similar  expressions  for  y,  z  and  the  conjugates.  Then: 

A“‘(r)  -  A"11  A.^-tr)  -A."1" 

A^2(r)  -  A^“2  A^V)  -  A^2* 

A1  ^  =  A1  A1  ^  =  A.l  (60) 

Equation  (57)  becomes 

k3x  3AjU3(r)/9x  +  k3y  3Aiu3(r)/3y  +  k3z  3Ai‘j3(r)/3z 

m  -1  3  U0  u3  j  k  1  Xijkl  Aj  Ak  A1  exp(-iAk-r)  (61) 
where  Ak  =  +  kc,  ~  The  solution  to  equation  (61)  is  assumed 

to  be  in  the  form: 
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then: 


> 


u>3 


(r)  =  K  exp(-i  Ak*r)  +  C 


3A^^(r)/3x  =  -iAkx  K  exp(-iAk'r) 
with  similar  expressions  for  y  and  z.  Equation  (61)  becomes: 

-i  K  (Ak*k3)  exp(-iAk-r)  =  'i  3  ^  ^  j  k  1  xijkl 

AjW^  A^3  A exp(-iAk'r)  (62) 

so  the  constant,  K,  is  in  the  form: 


0  2  I  I  I  ,  ul  ,  u2  ,  u4* 

K  =  _ 

Ak-  k^ 


(63) 


Again,  as  in  the  second  order,  if  there  is  no  interaction  length 
within  the  medium,  no  nonlinear  polarization  occurs  and  no  new  wave  can 
be  produced: 

Aiu3(0,0,0)  =  0  =  K  +  C 
C  =  -K 

The  solution  to  equation  (61)  is: 

A3u3(r)  =  K[exp(-i  Ak'r)  -  l]  (64) 

By  similar  development,  the  conjugate  solution  is: 

A^3  "(r)  =  «"  [exp(iAk*r)  -  1]  (65) 


where 


K  = 


J0  W1 


2  l  I  l 

j  k_l__Xijkl 

Ak*  k0 


A.ul*  A, w 2* 
1  k 


A  u 
A1 


(66) 


Again,  from  Poyn ting's  Theorem  (2:13),  the  i^  mode  power/unit  area 
at  cj3,  at  some  interaction  length,  k^L  (in  the  direction  of  the  wave, 
k^),  is  given  by: 
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P^/area  =  (l/2)(e/u())1//2Aiu3(k3xL,  k3yL,  k^L) 

Ai  ('k3xL’  k3yL’  k3zL') 

=  (L2/2)  (e/u0)1/2  M02  o32 

C1  Z  Z  y  A  wl  A  u2  A  w4*') 

'-j  k  1  xijkl  j  Ak  A1  ' 

(111  ul*  .  m2* 

k  1  xijkl  Aj  Ak  1  ; 

[sin2(Ak-k3  L/2)/  (Ak-k3  L/2)2]  (67) 

Again,  the  power  in  the  response  wave  (at  ^3)  is  maximized  when  the 
argument  of  the  (sin  x)/x  term  approaches  zero.  This  implies  Ak*k3  =  0 
and,  therefore,  the  conservation  of  momentum  condition,  Ak  =  0. 

Similar  developments  may  be  made  for  other  sum  and  difference  input 
frequencies  with  similar  results. 

Summary 


As  in  the  second  order  case,  the  light  wave  produced  via  the  third 
order  nonlinear  polarization  is  of  a  frequency  which  is  the  result  of 
mixing  three  input  frequencies  (in  some  sum  or  difference  combination) 
meeting  the  conservation  of  energy  condition.  The  direction  of 
propagation  of  this  wave,  k3  is  determined  by  meeting  the  conservation 
of  momentum  condition,  Ak  =  0,  so  that  energy  transferred  to  this  newly 
created  beam  is  significant. 
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III.  Optical  Phase  Conjugation  by  Degenerate  Four  Wave  Mixing  (A  Third 

Order  Nonlinear  Process) 


In  the  previous  chapter,  the  third  order  concepts  of  nonlinear 
polarization,  Maxwell's  wave  equation  and  phase  matching  were  developed. 
A  special  case  of  these  third  order  concepts,  optical  phase  conjugation 
by  degenerate  four  wave  mixing,  will  now  be  described. 

Figure  4  shows  the  geometry  to  be  considered.  EL^  and  are 
counterpropagating,  non-depleted  pump  beams.  E^  is  the  probe  beam  and 
is  the  conjugate  of  the  probe.  It  is  shown  that  in  order  to  meet 
phase  match  conditions,  the  conjugate  and  the  probe  must  also 
counterpropagate.  Expressions  for  the  amplitudes  of  the  conjugate  and 
probe,  which  extract  energy  from  the  pumps,  are  also  derived. 


Figure  4.  Geometry  of  phase  conjugation  by  four  wave  mixing  (3:29). 


Third  Order  Polarization 

The  polarization  of  the  nonlinear  medium  due  to  the  superposition  of 

the  four  electric  fields  is  given  by  equation  (1): 

Z  Z  Z  Z  Z  Z 

P.  ■  c  x.  •  E.  +  2  .  ,  d.  E.  E,  +  4  .  ,  ,  X..,,  E.  E,  E,  +  ... 
l  0  j  ij  j  ]k  ijk  j  k  1  ljkl  j  k  1 

The  third  order  nonlinear  polarization,  equation  (46),  is: 

„  z  z  z 

1  j  k  1  ljkl  j  k  1 
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The  total  electric  field  consists  of  the  four  superimposed  fields  so 
that: 


where 


=  4 


ZEE 

j  k  1  xijkl 


E  E  E  E  Z  E 
j  k  1  m  n  q 


(Eul 

J 

(E  ul 

(E^1 

xijkl 


+  E.Ui-  +  E.wp  +  E.1^) 
J  J  J 

+  Eku2  +  Ekwp  +  E^) 

+  E^2  +  ELup  +  EpUC) 
£  urn  g,  on  £  uq 

j  k  1 


(68) 


j ,  k  and  1  =  x ,  y  or  z 
m,  n  and  q  =  1,  2,  p  or  c 

(The  expansion  of  equation  (68)  yields  1728  terms!) 

The  product  of  the  E-fields  in  equation  (68)  can  be  further  expanded 
as  follows: 


um  ^un  E^uq  =  (1/2)  {A  expfi^t  -  k^r)] 


+  A^  "(r)  exp[-i(cjmt  -  k^r)]} 

(1/2)  (AkUin(r)  exp[i(Unt  -  k^.r)] 
ton*/ 


+  Akw“'(r)  exp[-i(unt  -  k^r)]} 

(1/2)  {A1U)q(r)  exp[i(wqt  -  k  *r)] 

+  A-^^'(r)  exp[-i(uqt  -  kq*r)]}  (69) 

(The  expansion  of  equations  (68)  and  (69)  yields  13,824  terms!) 

Thus  the  polarization  contains  many  frequency  components.  However, 
the  conservation  of  energy  condition  of  interest  here  is 
gjc  =  (jj  +  ~  Wp  so  the  frequency  components  of  the  nonlinear 

polarization  at  uc  are  singled  out: 

P ,uc  =  ojI  +  co2  up(  t)  =  3  E  E  E  A  ul  A  d>2  A  «P*(r) 

l  ’  j  k  1  * ljkl  1  j  k  1 

exp[i(u1  +  u2  -Up)t  -  (k^  +  -  kp)*  r] 

+  A.wl*  A*2*  ApU>p(r) 

exp(  - i [  (oo  1  +  u>2  ~  up)t  -  (k^  +  1^  -  kp)-  r]} } 

(70) 
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uc  =  wc  +  <4n  -  iJm 


(r,t)  3  j  k  1  Xijkl 

U.^r)  A^r)  A^Cr) 

k  -r)] 


exp[i(wct 
U XT* 


,  A  UC'V  \  »  UXTPr / _ \  •  *-i4 1 1  /  \ 

AAJ  <r>\  (r)Al  <r) 
exp[-i(uct  -  kc • r ) ] } 

where  m  =  1,  2,  p  or  c,  and  A^*,  ,  A^3  and  A^3  are  not 

functions  of  r  because  of  the  non-depleted  pump  assumption.  The 
amplitudes  of  the  polarization  wave  are  (as  in  equation  (53)): 


P.(uc>(r)  -  6  E  f  f  x....  A  ul  A  “2  A.“P*(r) 
i  v  j  k  1  Aijkl  j  k  1 

z  z  z  z 


,  ,  l  L  L  L  CJC/  \  .  (411/  \  . 

+  6  .  ,  ,  X-  -i  i  A.  (r)  A,  (r)  A-, 

j  k  1  m  ljkl  j  7  k  1 


um* 


(71) 


(r) 


p  (uc)'-/  a  0  1  Z  Z 
i  b  j  k  1  xijkl  J 

L,  nn 

+  6  .  ,  ,  X-  i  A- 
j  k  1  m  ljkl  j 


A,"1*  A,“2*  Aj^tr) 


k 

U)1 « 


(r)  A^Xr)  A1Mm(r)  (72) 


The  Wave  Equation 

The  wave  equation  developed  from  Maxwell's  equations,  equation  (21), 
is  in  the  following  form: 

{-i[kcx  aAiLlX'(r)/3x  +  kcy  3AjL‘JC(r)/ay  +  kcz  aAi^(r)/9z] 

-  (1/2)  |kc|2  ApWC(r) }  exp[i(w,t-  kc*r)] 

+  { i[kcx  3AilJC''(r)/ax  +  kcy  aAiUC“(r)/3y  +  kcz  aA^'Xr/^z] 

-  (1/2)  |kj2  AilJC”(r)}  exp[-i(wct  -  k.‘r)] 

=  (i  ooc/2)  (Aiur(r)  exp[i(wct  -  kc*r)] 

-  A^'Xr)  e  p[i(wct  -  k,*r)]} 

-  (uQ  e  u,2/2)  lAiuc(r)  exp[i(oct  -  kc*r)] 

+  A^'Xr)  exp[-i(uct  -  kc*r)]} 

+  m0  o2PiuC(r,t)/at2  (73) 

From  equations  (60)  and  (71): 

a2Pi(JC(r,t)/at2  =  -wc2Piwc(r,t)  (101) 

The  same  simplifying  assumption  made  in  the  previous  chapter  is  made 
here:  the  medium  is  transparent  at  w  so  that  o  =  0.  Then,  using 
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2  2 

|kj  =  wc  Mq  e ,  terms  on  either  side  of  equation  (73)  synchronous 
with  exp(i  w  t)  can  be  equated: 

-i^kcx  +  kcy  3AiL>C(r)/3y  +  k£z  aAiWC(r)/3z]  exp(-i  kc*r) 


Z.  n  U)C/  .  \ 

-_y0uc  Pi  (r’t} 


kcx  3Aiwc(r)/3x  +  kcy  SA/^r)/^  +  kcz  3AiLJC'(r)/3z 


cy 


■  r  2  J  Z  I  Z  .  Wi  .  wt  A 

=  -l  fa  pA  u  .  ,  ,  v.  ,  A.  A,  A, 
0  c  j  k  1  *nkl  j  k  1 


wi  0)2  up* 


(r) 


exp[-i(k^  +  -  k  -  kc)*r] 


,1111  »  U)C  /  \  «  OJTT  -  \  . 

+  j  k  1  m  Xijkl  Aj  (r)  A1  lr>  A1 


p 


(75) 


Since  Ej(r,t)  and  E2(r,t)  are  counterpropagating,  k^  +  =  0 

(2:555).  (For  an  argument  when  k^  +  /  0,  consult  Reference  6  (pages 

30-33).)  Then,  if  k -  -  k  (the  conjugate  is  counterpropagating  to  the 

c  p 

probe),  the  exponential  in  equation  (75)  may  be  set  to  unity. 

Another  simplification  made  toward  achieving  a  solution  is:  the 
amplitudes  of  the  pumps  are  much  greater  than  that  of  the  probe  or 


conjugate: 

ul  »  0)1* 


Ak  Al  +  V  Ai  »  Akup(r)  A1wp*(r)  +  AwC  (r)  AwC'"(r)  (76) 
Then,  equation  (75)  reduces  to: 

kcx  3AiWC(r)/3x  +  k  3AiwC(r)/3y  +  k,z  3Aiu,c(r)/3z 


•  r  2  E  II  .  0)1  .  0)2  .  up*/  \ 

■l  6  u0uc  j  k  1  Xijkl  Aj  Ak  A1  (r) 


+  Ajuc(r)[Akwl  A^1*  +  Aku2  A^2*] 

(77) 


The  third  simplification  made  is:  since  E^JC(r)  and  E^p(r)  are  the 
only  terms  which  are  a  function  of  position  and  they  are  a 
counterpropagating  pair,  let  both  travel  in  z,  only.  Then,  equation 
(77)  becomes: 

ul  A  u?  d  up* 

Hjkl 
uc 


dAitjC(z)/dz  -i6uc  (u0/c)  ^  j  k  1  ^ijkl  ^AjW  Ak^  A1U)P  ^ 


+  A  ^(z)[Akul  A^1”  +  Ak“z  A^"]}  (78) 


The  final  simplifying  assumption  is:  often,  one  Xjjki  will  be 


much 
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greater  than  any  of  the  others  in  the  tensor.  Assume  the  fields  here 


are  polarized  along  the  appropriate  axes  to  take  advantage  of  that 

(3) 

dominant  susceptibility,  xv  .  Equation  (78)  becomes: 
dAU)C(z)/dz  =  -i  6  (u0/e)1/2  [A^  Aw2  Aup'(z) 

+  AUC(z)  (|Aul|2  +  |AUV)]  (79) 


(3:31) 

By  a  similar  development,  the  equation  for  E^P(z)  is: 
dAup(z)/dz  =  i6  u)(Uq/c  )^2  x^[AW*  Aw2  Auc"(z) 

+  AUP(z)  (|Awl|2  +  |A-2|2)]  (80) 


(3:30) 


Solutions  for  the  Amplitudes  of  the  Probe  and  Conjugate 

Assume  the  forms  of  the  solutions  to  equations  (79)  and  (80)  to  be: 
Awc (z)  =  Auc,(z)  exp[-i  6  U(p0/c)1/2  X(3)  (|Aul|2  +  |Au2|2)z] 
Awp(z)  =  Aw,J,(z)  exp[ i  6  u(uQ/ x^  (|Awl|2  +  |au>2|2)z]  (81) 

then: 

d.AlJC(z)/dz  =  (dAWC' (z)/dz)  exp[-i  (^(Uq/e)^2  x^  (|aW^|2  +  |aw2|2)z] 
-i  6  t j(m0/c)  X  ( I A  i  +  I A  I  )  A  (z) 

d.AU)P(z)/dz  =  (dAup' (z)/dz)  exp[i  6  u(uqA:)^2  x^  (|Aw^|2  +  |AU)2|2)z] 
+  i  6  u)(m0/c)1/2  X(3)  (  I A^1 1  2  +  |Au2|2)  AUP(z)  (82) 

(3:33) 

Inserting  equations  (82)  into  equations  (79)  and  (80)  gives: 

-i  6  o(u0/e )  x(3)  Awl  Au2  Aup*(z) 

=  dAwC'(z)/dz  exp[-i  6  u(uq/c  )1,/2  x^  (lA^1)2  +  |Aw2|2)z] 
dAwC,(z)/dz  =  -i  6  W(uQ/e)1/2  X(3)  Awl  Aw2  Awp,*(z) 

=  i  hAup'*(z)  (83) 

and: 
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i  6  u(n0/e)1/2  x(3)  A^1  A^2  AUc“(z) 

=  dA^P' (z)/dz  exp[i  6  u(Pq/e)1/2  x^3'*  ( I A^1 1  2  +  |aw2|2)z] 
dAuP,(z)/dz  =  i  6  Uu0/c)U2  x(3)  AW1  AW2  A^’dz) 

=  -i  hAUC,"(z)  (84) 

where 

k=  -6  lj(u0/c)1/2  X^3^  AUl  Au2 

Note  in  equation  (81)  that  the  primed  notation  merely  represented  a 
phase  shift.  This  phase  shift  will  be  inconsequential  so  the  primed 
notation  is  dropped  here.  Then,  equations  (83)  and  (84)  become: 

dAuc(z)/dz  =  i  x  A^'iz) 

dAUp(z)/dz  =  -i  h  (z)  (85) 

(2:556,  3:36). 

Yariv  gives  the  solution  to  equations  (85)  (2:556,  3:36)  as: 

Aup(z)  =  Auip(0)  cos[  |  x  |  (z  -  L)]/cos(  [  h  | L) 

AtX(z)  =  i  A^'(O)  h"  sin[|x|  (z  -  L)]/(|*|  cos  Ml)  (86) 

using  as  boundary  conditions  that  E  p(0)  is  finite  and  E  C(L)  =  0.  The 
value  of  the  probe  wave  at  the  output  plane  is  given  by: 

A^l)  =  Aup(0)/cos(  |h|l)  (87) 

whereas  the  value  of  the  conjugate  wave  reflected  at  the  input  plane 
is : 

Auc(0)  =  -i  A^dO)  (h'/M)  tan  (  |h|l)  (88) 

Figure  5  represents  equations  (86),  (87)  and  (88)  graphically. 
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Figure  5.  Intensities  of  the  probe  and  conjugate  waves  in  phase 

conjugation  by  degenerate  four  wave  mixing  (2:559,  3:37). 

Summary 

Optical  phase  conjugation  by  degenerate  four  wave  mixing  was  shown 
to  be  a  special  case  of  the  third  order  nonlinear  processes  discussed  in 
the  previous  chapter.  By  applying  the  concepts  of  nonlinear 
polarization  and  the  wave  equation,  it  was  shown  that  the  probe  and 
conjugate  must  counterpropagate  to  achieve  phase  match  conditions  and 
significant  energy  transfer  from  the  pumps  to  the  probe  and  conjugate 
waves.  Finally,  expressions  for  the  amplitudes  of  the  probe  and 
conjugate  as  functions  of  position  in  the  nonlinear  medium  were  derived. 
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IV.  The  Photorefractive  Effect 


Introduction 


The  previous  chapter  showed  that  optical  phase  conjugation  was 
possible  by  degenerate  four  wave  mixing  in  which  the  conjugate  light 
wave  was  generated  by  the  third  order  nonlinear  polarization  of  the 
medium.  The  amplitude  of  this  polarization  is  directly  proportional  to 
the  third  order  nonlinear  susceptibility  of  the  medium  and  the  amplitude 
of  each  of  the  incident  electric  fields  (equation  (72)). 

Table  1  lists  the  dominant  third  order  susceptibilities  for  some 
exemplary  media.  Because  the  magnitude  of  these  susceptibilities  is  very 
small,  high  intensity  incident  light  is  required  to  drive  the  medium 
into  nonlinearity  and  produce  phase  conjugation.  Therefore,  an 
alternative  technique  for  producing  optical  phase  conjugation  with  lower 
intensity  incident  light  would  be  desireable.  Such  a  technique  is 
through  the  use  of  photorefractive  crystals. 

Table  I 

Examples  of  third  order  nonlinear  susceptibilities,  X  ,  at  a 
wavelength  of  694  nm  (2:554) 


Material 

n 

X., <*>.  -«) 

(10- M  MK.S) 

cs- 

1.612 

441 

CCT„ 

1  454 

6.2 

Fused  quartz 

1.455 

1.5 

YAG 

1.829 

7.41 

Benzene 

1  493 

68  9  . 

LSO  glass 

1.505 

2.26 

ED-4  glass 

1.557 

2.8 

SF-7  glass 

1.631 

9.8 

BK-7  glass 

1.513 

2.26 

LaSF-7  glass 

1.91 

12.4 

The  interference  of  the  (coherent)  pump  and  probe  beams  (a.k.a. 


"write"  beams)  in  a  pho tore tractive  crystal  causes  a  spatially 
sinusoidal  distribution  of  intensity  as  shown  in  Figure  6-a.  A 
photorefractive  crystal  contains  various  sites  in  its  lattice  which 
accept  and  hold  electrons  with  a  low  thermal  probability  of  excitation. 
(The  origin  of  these  sites  is  uncertain.  They  may  be  naturally  occurring 
impurities,  doped  impurities  or  defects  in  the  crystal.)  Under 
illumination,  these  charges  are  freed  and  the  movement  of  either 
electrons  or  holes  (depending  on  the  crystal)  occurs  until  the  charge  is 
retrapped  in  an  area  of  the  crystal  that  is  not  illuminated.  In  Figure 
6-b,  the  location  from  which  charges  left  corresponds  to  the 
interference  maxima  of  Figure  6-a  and  the  location  where  they  are 
retrapped,  to  the  interference  minima. 


a)  Intensity 


b)  Charge  separation 
(holes  are  carriers) 

c)  Electro-static  field 
(positive  referenced 

to  the  left) 

d)  Change  in 
refractive  index 


Figure  6.  Illustration  of  the  photorefractive  effect  (4:422) 


The  separation  of  opposite  charges,  either  electrons  or  holes 
(whichever  is  the  carrier)  and  the  ions  left  behind,  causes  an 
electro-static  field  to  form  (Figure  6-c).  This  field  modulates  the 
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index  of  refraction  of  the  crystal  through  the  linear  electro-optic 
(Pockel's)  effect  (Figure  6-d)  (2:223,226). 

For  this  "photorefractive  effect",  the  electric  field  driving  the 
nonlinearity  of  the  crystal  is  the  created  space  charge  field.  The 
magnitude  of  this  field  is  proportional  to  the  modulation  index  of  the 
interfering  beams.  Since  the  modulation  index  is  dependent  on  relative 
intensities  rather  than  the  incident  intensity  of  either  beam,  phase 
conjugation  can  be  produced  from  lower  intensity  incident  light. 

Phase  conjugation  in  photorefractive  crystals,  then,  differs  from 
that  of  degenerate  four  wave  mixing  in  that  rather  than  being  dependent 
on  the  INTENSITY  of  the  incident  beams,  it  is  dependent  on  the  RE1ATIVE 
intensity  of  the  pump  and  probe  beams. 

Two  models  describing  the  photorefractive  effect  are  Feinberg's 
Hopping  Model  and  Kukhtarev's  Solid  State  Model.  Feinberg's  model  is 
based  on  the  idea  that  a  charge  "hops"  from  site  to  site  in  the  crystal 
and  develops  the  probability  of  this  charge  movement.  Kukhtarev's  model 
is  a  more  standard  solid  state  development  of  charge  transport. 

Both  models'  solutions  for  a  given  sinusoidal  illumination  with 
small  modulation  index  are  a  sinusoidal  space  charge  field  of  the  same 
spatial  frequency  as  the  illumination  but  out  of  phase  with  it.  This 
space  charge  field  is  the  basis  for  the  development  of  "two  beam 
coupling"  and  diffraction  efficiency  expressions  in  later  chapters. 

Feinberg's  Hopping  Model 

The  Model  Itself 

Ihe  "hopping"  model  is  based  on  the  idea  that  a  charge  may  occupy 
any  of  a  large  number  of  sites  in  the  crystal  lattice.  When  the  crystal 
is  not  illuminated,  a  charge  is  bound  to  a  site.  However,  under 
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illumination,  the  charge  is  given  enough  energy  to  free  it  from  any  one 
site  and  it  is  free  to  "hop"  from  site  to  site.  When  the  charge  is 
re-bound  at  another  (dark)  site,  the  separation  of  the  charge  and  the 
ionized  site  left  behind  causes  an  electro-static  field  to  be  set  up 
which  will  modulate  the  index  of  refraction  of  the  crystal. 

Assumptions  which  are  implicit  in  the  model  are: 

>  1.  Typical  photorefractive  material:  charge  is  bound  deep  in  a  site  in 

darkness  but  is  free  to  hop  to  another  site  under  illumination.  (5:1298) 
In  other  words,  the  charge  cannot  be  thermally  excited  enough  to  hop  from 
site  to  site.  The  implication  cf  this  is  that  once  a  grating  is  formed 
in  the  crystal,  if  the  crystal  is  kept  in  darkness  the  grating  will 
remain  indefinitely.  Although  the  grating  does  persist  in  darkness,  it 
does  eventually  decay.  (5:1297) 

2.  The  number  of  sites  that  these  charges  can  occupy  is  much  greater 
than  the  number  of  charges  available  (5:1298).  Thus,  the  model  does  not 
account  for  saturation  effects.  A  charge  is  assumed  to  always  have  a 
place  to  migrate  to. 

3.  The  number  of  charge  carriers  excited  into  the  conduction  band  is 
proportional  to  the  incident  intensity.  Therefore,  the  hopping  rate  is 
also  proportional  to  the  intensity  (5:1298). 

Definitions  necessary  before  setting  up  the  model  are: 

=  probability  that  a  site  at  position  xn  is  filled 
dW^/dt  =  probability  per  second  of  hopping  into  that  site 

Then: 
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oc 


dW^/dt  «  probability  that  the  site  is  filled  and  the  charge  has 

t  h  th 

enough  energy  to  hop  to  the  n  site  (and  the  n  site  is 

tH 

open)  minus  the  probability  that  the  nLr  site  is  filled  and 

th 

the  charge  has  enough  energy  to  leave  (and  the  m  site  is 
open) . 

“  ‘m  lV.  exP[C  P  (tm  '  V/kB  Ti  (1  '  Wn> 

-  In  Wn  exp[C  q  (♦„  -  »m)/kB  T]  (1  -  Km) 

where 

I  ,  \  =  Intensity  at  site  n(m) 
n(m)  } 

C  =  Constant  to  assure  model  obeys  statistical  mechanics 
q  =  Charge  (Note  that  it  is  a  positive  charge  in  this  form.  It  would 
be  -q  for  an  electron.) 


<J>  n(m) 


Potential  at  site  n(m) 


kg  =  Boltzmann's  constant 
T  =  Lattice  temperature 

If  the  number  of  sites  is  much  greater  than  the  number  of  charges, 


then  Wn(  v  <<  1  and  the  rate  equation  becomes: 

dW  /dt  =  -  Z  D  [I  W  exp(C  3  $  )  -  I  W  exp(C  B  *  )] 
n  m  mn  n  n  nm  m  m  r  mn 

where 


(89) 


Dmn  =  Constant  of  proportionality 
B  =  q/(kg  T) 

^nm(mn)  ^n(m)  %(n) 

To  determine  C,  under  uniform  intensity,  the  probability  that  a  site  is 

filled  should  be  the  same  for  all  sites: 

W  exp(8  $  )  =  W  exp(B  $  ) 
n  K  n  m  K  m 


W  /  W  =  exp(B  $  ) 

n  m  v  nm 

(90) 

and  the  probability  of  movement  should  be  zero: 

dW  /dt  =  0 
n 

(91) 
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Then,  from  equations  (89)  and  (91): 


W  exp(C  B  4>  )  =  W  exp(C  B  <J>  ) 

n  K  nm  m  K  mn 

W  /  W  =  exp(2C  B  t  ) 
n  m  r  nm 

Therefore,  from  equation  (90): 

C  =  1/2 

(5:1298) 

The  model  is  expressed  mathematically  as  (5:1298): 

dW  /dt  =  -  ^  D  [W  I  exp( B  <t>  /2)  -  W  I  exp( B  0/2)] 
n  m  mn  L  n  n  F  nm  mm^  rnn 


(92) 


(93) 


Solution  to  the  Model 

More  assumptions  made  to  achieve  a  solution  to  the  model  are: 

4.  Sites  are  equally  rather  than  randomly  spaced.  The  rms  distance 
between  sites  is  1  (5:1299).  This  allows: 

\  \  -  1  <94> 

5.  There  are  many  adjacent  sites  (and  hopping  opportunites)  within  each 

spatial  wavelength  (kl  <<  1).  This  allows: 

Dmn  =  0  except  for  nearest-neighbor  sites 

=  D  for  nearest-neighbor  sites.  (95) 

This  removes  the  summation  from  the  hopping  probability  equation. 
(5:1299) 


The  model  is  now;  in  the  form: 

dW  /dt  =  -D  [W  I  exp(8  <t>  /2)  -  W  I  exp(-B  $  /2 ) ]  (9( 

n  Lnnrnm  mmr  nm 

Consider  this  special  case:  Two  light  waves  interfere  as  in 
Figure  7.  ITie  electric  field  of  each  wave  can  be  written  as  (5:1298): 

Ej  =  e^  exp(ikj’x) 

^2  =  ^2  ^2  exP(i*c2’x)  (9- 

and  the  interference  as  (5:1299): 

rn(m)  =  S  [1  +  m  cos  k'xn(m)]  (9i 
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where 


i 


m  =  modulation  index 
=  2  Ej  E2  e1  e2  /GQ 
=  G/Gq 

k  =  kj  - 

Further  assumptions  toward  a  solution  are: 

6.  Small  modulation  conditions  will  only  be  considered  (5:1299): 

m  <<  1  (99) 


7.  The  form  of  the  solutions  (following  the  form  of  the  interference 
above)  is  (5:1299): 

W  =  UT  +  Re  W  exp(ik*x  ) 
n  On  K  n 

♦  =  0^  +  Re  0  exp(ik*xn)  (100) 

where 

=  Unmodulated  probability  site  n  is  filled 

=  Potential  due  to  external  or  intrinsic  fields  (everything  but 
charge  migration) 

W  =  Charge  distribution  amplitude 
$  *  Potential  distribution  amplitude 
Kukhtarev  will  show  assumption  7  is  good  if  assumption  6  already  exists. 
However,  both  Feinberg's  and  Kukhtarev' s  models  are  often  extended  to 
cases  where  the  modulation  index  approaches  unity  and  the  resulting 
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field  is  still  assimed  to  be  spatially  sinusoidal  following  the 
interference  fringes.  The  field  may  not  be  sinusoidal  under  such 
conditions.  ITiis  problem  is  addressed  to  some  extent  in  the  development 
of  Kukhtarev's  model. 

Equations  (100)  are  not  independent.  From  Poisson's  equation 
(6:164): 

V  •  D  =  pv  q  (101) 

where 

D  =  Electric  displacement  vector 
=  eQ  e-E 

Cq  =  Permittivity  of  free  space 
e  =  Static  dielectric  tensor 
E  =  -  V  <J> 

V  =  ik  (for  plane  waves) 
qo  ,  =  Charge  distribution 
=  p  W  q 

p  =  Average  density  of  sites  (number  per  unit  volume) 

Then,  the  charge  and  potential  distributions  are  related  by: 

-ik-e- ik  $  =  p  q  W/  0 

<t>  =  p  q  W/e  g  q  k^  (102) 

(5:1299) 

Final  simplifying  assumptions  made  toward  achieving  a  solution  are: 
8.  First  order  series  expansion  of  exponential  will,  be  used  (5:1299): 

exp(B  4>  /2)  =  1  +  B  $nm/2  (103) 

run  mil 

This  is  reasonable  since  (3  v  1/40  eV  at  room  temperature,  typical 
internal  fields  are  <  50  kV/cm  and  1  v  1/10^  cm  so  B  <t>__/2  <<  1. 

ion 

The  hopping  probability  equation  (96)  becomes: 


d/dt  [V0  +  Re  W  exp(ik*xn)]  =  -D  {[Wq  +  Re  W  exp(ik*xn)] 

[Gq  +  Re  G  exp(ik*xn)] 

{l  +  (B/2)  [A<J>q  +  Re  9  exp(ik*xn) 

-  Re  0  expfik’Cx^  +  1)]} 

-  {WQ  +  Re  W  exp[ik*(xn  +  1)]} 

{Gq  +  Re  G  exp[ik*(xn  +  1)]} 
jl  -  (B/2)  { A<tg  +  Re  9  exp(ik*xn) 

-  Re  $  exp[ik*(xn  +  1)]))}  (104) 

where 

A*0  =  *0n  ‘  $0m 

Equating  synchronous  terms  [in  exp(ik*xn)]  from  both  sides  (5:1299): 
dW/dt  =  -D  {(W0  Gq  +  W  G/2)(B  9/4)[l  -  exp(ikl)] 

+  (WQ  G/2  +  W  Gq/2)(1  +  BA  9Q/2)  +  (W  G/4)(B  9/4) 

[l  -  exp(-ikl)] 

+  (WQ  Gq  +  W  G/2)( B  9/4 )[ 1  -  exp(ikl) ] 

-  (Wq  G/2  +  W  Gq/2)  exp(ikl)  (1  -  BA  9Q/2) 

+  (W  G/4)  exp(2ikl)  (8  9/4)[l  -  exp(-ikl)]  }  (105) 

After  expanding  the  exponentials  to  their  second  order  (reasonable  given 
assumption  5): 

exp(ikl)  =  1  +  ikl  -  (kl)2/2  (106) 

and  making  the  substitutions: 

w  =  Normalized  charge  wave  amplitude 
-W/U0 
m  =  G/Gq 

Feinberg  gives  the  final  form  as  (5:1299): 

dw/dt  =  -DG0  [(w  +  m)(k2  l2  -  ikl  BA  9Q)  +  9B  k2  l2] 

=  -r[(w  +  m)(a2  +  iaf )  +  w]  (107) 

where 
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*0 


Characteristic  hopping  rate 


Characteristic  wave  vector 
[(p  W Q  q2/(e  e0  kg  T)1/2] 
Normalized  grating  wave  vector 


=  k/kQ 

f  =  Uniform  e-field  strength  normalized  by  the  characteristic  field 

■  -‘V(1  f0> 

fg  =  Characteristic  field 
=  kQ  kB  T/q 

The  solution  to  this  differential  equation  (107)  is: 
dw 

- - - - - =  -p  dt 

v(a  +  1  +  iaf)  +  m(a^  +  iaf) 

2  2 

exp[-r  (a  +  1  +  iaf)  t]  -  m  (a  +  iaf) 

w  . . 2 . 

a  +  1  +  iaf 

(108) 


or  solving  for  the  potential  amplitude,  $,  using  equation  (102): 


W, 


$  = 


0  D 


q  [exp[-F  (a  +  1  +  iaf)  t]  -  m  (a  +  iaf)] 


e  Eq  k2  (a2  +  1  +  iaf) 


(109) 


Electric  Field  Modulating  the  Index  of  Refraction 
In  the  linear  electro-optic  effect,  the  index  of  refraction  of  a 
crystal  may  be  modulated  by  an  internal  dc  electric  field.  In 
photorefractive  crystals,  this  (space  charge)  E-field  is  the  (negative 
of  the)  gradient  of  the  quasistatic  potential  due  to  internal  charge 
migration. 

From  equation  (109),  the  steady  state  solution  for  the  potential 
amplitude  is: 
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(110) 


2 

-Wq  p  q  m  (a  +  iaf) 

- *---* - 

e  Eq  k“  (a  +  1  +  iaf) 

-fn  m  a  +  if 

~g~  (-J . ) 

k  a  +  1  +  iaf 


The  space  charge  field,  is  (5:1304): 

Ec;c(x)  =  -V  <t>(x) 

=  -V  [<3>q  +  Re  <t>  exp(ik*x)] 
a  +  if 

=  i  It  L  m  (-* - )  exp(ik’x)] 

U  a  +  1  +  iaf 


(111) 


Note  that  the  space  charge  field  should  be  in  the  direction  of  k  given 
the  form  of  the  potential  (equation  (100)). 

Since  this  is  a  complex  expression,  it  may  be  expressed  in 
magnitude-phase  form.  The  significance  of  this  form  is  evident  in  later 
chapters  where  the  diffraction  efficiency  and  phase  conjugate 
reflectivity  are  dependent  on  the  magnitude  of  and  the  energy 
transfer  in  steady  state  two  beam  coupling  is  dependent  on  the  phase 
shift  between  the  grating  and  the  intensity  fringes. 


Eq  =  Limiting  value  of  the  space  charge  field 
=  f0/a 

=  p  Wq  q/( e  Eq  k)  (Proportional  to  the  trapped  charge  density) 
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Ep  =  Diffusion  field  (assumed  parallel  (or  antiparallel)  to  k) 

=  fQ  a 
=  kB  T  k/q 

Eq  =  Applied  field  (also  assumed  parallel  (or  antiparallel  to  k) 

=  f  f0 

To  find  the  phase  shift,  >1,  between  the  space  charge  field  and  the 
intensity  interference  pattern,  the  space  charge  field  is  separated  into 
its  real  and  imaginary  parts: 

-f  +  ia 

Ecr(x)  a  ———————————— 

ct  +  1  +  iaf 

f  +  ia  (1  +  a2  +  f2) 
a  —  * - * - *— 

(a2  +  1 )  +  (af)2  (113) 

i  =  tan  ^  a  (1  +  a2  +  f2)/f 

.  -1  ^  Eq  +  ^  +  E02 

E0  Eq  (114) 

With  no  applied  field  (Eq  =  0),  tan  4'  =  °°,  and  the  phase  shift  is  V 2 . 

As  the  applied  field  grows  large  (Eq  =  °°) ,  tan  4<  =  00  ,  and  the  phase 
shift  is  again,  tt/2.  Again,  these  forms  are  used  in  the  papers 
mentioned  above. 

The  space  charge  field,  then,  can  be  written  in  the  form: 

2  2 

EL  +  £ 

F,r(x)  =  £  m  E  [ - - L^2  ke  exp(i4)  exp(ik*x) 

q  (ED  +  V  +  E0  (H5) 

This  is  the  form  for  the  space  charge  field  found  in  other  papers 

such  as  those  by  Fischer  et  al.  (7:520),  Valley  and  Klein  (8:115),  and 
Gunter  (9:252). 

Modulated  Optical  Susceptibility 
The  space  charge  field  causes  a  modulation  of  the  optical 


susceptibility,  X(x),  of  the  form: 


X(x)  =  Re  X  exp(ik’x) 


(lib) 


The  derivation  for  the  form  of  the  amplitude  of  the  modulated 
susceptibility  follows: 

p  =  p  +  p 

*L  *NL 

=  e0  x(1^  ^  +  2  d*Esc-E('lJ') 


t  (1)  +  2  d-Esc/e0]-E( 


The  electric  displacement  vector  is: 
D  =  £q  E  +  P 


■c0E(“l)[0[/1b2d-5.c/i(1]-EM 
■■ol"  XU)  -  2  d-^c/£0]  E(“> 

-  c0[c*2d-ESc/£0]-E(“) 


(117) 


(118) 


=  I  +  x(1) 


(119) 


where 


D  =  e0  (b  +  Ab)"1  *E^ 


(120) 


b  =  Impermeability  tensor 

Ab=  R-E^ 

R  =  Linear  electro-optic  tensor 


Therefore: 


e  +  2  d  'EgC/ £g  (b  +  R'E^,) 


The  first  order  series  expansion  of  the  right  side  is: 

(b  +  R*ESC)_1  -  b-1  -  b-]  •R-Esc;b'1 


(121) 


(122) 


e  =  b 


(123) 


Therefore: 

2  d,Esc  •  -"o  e*R,Escc  (m) 

As  shown  earlier,  the  space  charge  field  is  related  to  the  potential  by: 


1 


Eg^(x)  =  -i  k  i  exp(ik*x)  (125) 

so  the  amplitude  of  the  modulated  optical  susceptibility,  X,  takes  the 
form  (5:1299): 

X  =  2  d>ESC 

=  i  Eq  e*R‘k*  e  0  (126) 

Summary 

Feinberg  et  al.  developed  a  model  for  the  photorefractive  effect 
based  on  the  probability  that  a  certain  site  in  the  crystal  lattice  was 
filled  by  a  charge.  The  fruition  of  this  model  in  the  case  of  spatially 
periodic  illumination  (interference  fringes)  shows  the  space  charge 
field  set  up  in  the  medium  for  small  modulation  fringes  is  spatially 
sinusoidal  with  the  same  frequency  as  the  illumination,  amplitude  given 
by  equation  (113)  and  shifted  in  phase  from  the  illumination  with  shift 
given  by  equation  (114).  This  field  modulates  the  optical  susceptibility 
of  the  medium  as  given  by  equations  (116)  and  (126). 


Kukhtarev  Solid  State  Model 


Kukhtarev  et  al.  (10:949-960)  developed  a  model  for  the 
photorefractive  process  based  more  in  conventional  solid  state  physics 
than  Feinberg 's  hopping  model.  Other  papers  such  as  Kogelnik 
(11:2909-2945),  Vahey  (12:3510-3515),  Fisher  (7:519-521)  and  even 
Feinberg  (5:1297-1305)  develop  theories  in  which  an  initial  assumption 
is  that  the  spatial  distribution  of  electrons  (and  thus  the  index 
grating)  is  proportional  to  the  incident  intensity  pattern.  Kukhtarev 
clearly  points  out  he  makes  no  such  assumption  to  begin  with,  yet,  in 
the  end  this  is  shown  to  be  the  case.  (10:950). 


ground  to  the  Model 


The  same  interference  pattern  as  Feinberg  considered  is  set  up  in 


the  crystal  as  in  Figure  8.  The  incident  waves,  2±1q>  are  assumed  to  ^ 
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polarized  in  the  plane  of  incidence  (10:951). 


Figure  8.  Geometry  for  formation  of  an  index  grating  in 
photorefractive  crystals  (10:951). 


Probabilities  considered  when  developing  the  model  are: 

Electron  transitions  from  donors 

to  conductivity  band  =  (si  +  B)(N^  -  Nrj  )  (127) 

where 

s  =  Cross  section  of  photoionization 

I  =  Light  intensity 

=  c  |ET|2  e1/2/ (4tt) 

c  =  Velocity  of  light 

e  =  High  frequency  dielectric  constant 

£.(1)  =  i*"^1  component  of  e  which  is  a  function  of  intensity  through 
the  linear  electro-optic  effect 

=  e.(l  +  Ac.) 

Ac  =  rik£iLk^1^  (i>k  =  x>  y  or  z) 

2  Ak 

_ l 

K 

1 

r^  =  Component  of  electro-optic  tensor,  R 
E^(I)  =  component  of  space  charge  field 
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=  i  component  of  refractive  index  (written  as  n  in  other 
places  and  elsewhere  in  this  paper) 


B  =  Rate  of  thermal  generation 
Np  =  concentration  of  donors 
Np+  =  Concentration  of  ionized  donors 


and 

Ionized  donors  capturing  electrons  =  'YpnNp+  (128) 

where 


yR  =  Recombination  constant 
n  =  Concentration  of  carriers 
(10:951) 


The  change  in  the  concentration  of  carriers  with  time  is  given  by: 

dn/dt  =  dND+/dt  -  Vj/e  (129) 

where 

dh'D+/dt  =  (si  +  B)(ND  -  Nd+)  -  YRnND+ 
j  =  epn  [E  -  (T/e)Vln  n]  +  ple^ 
e  =  Electronic  charge 
p  =  Electron  mobility 

E  =  Electric  field  modulating  the  index  of  refraction 
T  =  Temperature 
pi  =  Photovoltaic  current 
p  =  Photovoltaic  constant 
e^  =  Unit  vector  along  c-axis 
The  wave  equation  is: 

V2E  +  (1/c2)  d2(eE)/dt2  =  0  (130) 


(10:951) 


The  electric  field  of  the  two  interfering  beams  may  be  written  as: 
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where 


E  =  (Wc e1//2)  {e+1  I+11,/2  exp[i(4>+1  +  kxx  -  kzz)] 

+  e_1  I^1^2  exp[i(<J>_1  -kxx  -  kzz ) ] }  (131) 


e±^  =  Polarization  vectors 

k  =  Wave  vector  components 
x,z  r 

<J>  +  j  =  Phases  of  interfering  beams 

(10:951) 

The  electron  generation  rate,  including  both  thermal  and 
photo-generation,  is  given  by: 

f  =  B  +  si 

=  (B  +  sIq)[1  +  MQ  cos  |2kxx  +  <Kz)|] 

where 


(132) 


X0  X+1  +  L-1 

M0  -  XhiKO1'2  <e+1-e.,)/[I0(l  *  VV] 

=  m/(l  +  °d/°p) 
m  =  Modulation  index 

=  2(I+11_1)1/2  cos(0+1  +  O.jJ/Iq 
o^/Cp  =  Ratio  of  dark  and  photo-conductivities 
=  B/ ( s  1 0 ) 

<t(z)  =  Phase  difference  of  interfering  beams 

=  <6  -  <6 

>1  -1 

(10:951-952) 

Steady  State  Solution 

The  space-charge  field,  AE(x,I)  is  now  solved  for  and  shown  to  be 
sinusoidal  and  proportional  to  the  incident  interference  pattern. 

In  the  case  of  linear  recombination,  there  are  more  donors  than  free 
carriers  and  more  free  carriers  than  the  generation  rate-recombination 
ratio: 
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where 


Nd,  nd  “  Na  »  n  »  f Ar 


(133) 


*  Trap  concentration 
(10:952) 

The  spatial  change  in  the  E-field  is  given  by: 

dE/dx  =  [ p/(eM ) ]  d(l/n)/dx  +  (J/eu)  d(l/n)/dx  +  (T2/e)  d2ln  n/dx2 
=  (4tre/e0)(NA  -  fND/yRn) 

E  =  [J/(epn)j  +  (T/e)dln  n/dx  +  pl/(eun)  (134) 

where 

J  =  Constant  current  associated  with  the  averaged  e-field,  <E>, 
and  applied  voltage, V. 

<E>  -  1/d/  0d  E  dx 

-  JR  +  V/d 

d  =  Length  of  crystal  in  x-direction 
R  =  Resistivity  of  external  circuit 
(10:952) 

Case  1 

If  the  space  charge  screening  length,  Lp  and  the  length  of  electron 
tightening  by  the  external  field,  Eg,  1^,  are  much  greater  than  the 
fringe  spacing,  L: 

Lp  1E  »  L  (135) 

where 

Lf  -  [e0  tt  T  /(e2  Na)]1/2 

h  =  C0  E0  ^2eNA> 

Eg  =  V/d 

then  the  carrier  concentration  is: 


n  =  Nd(S  +  sI0)/(Nayr) 

=  nd  +  np 


(136) 
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and  the  space  charge  field  is  given  by: 

E  =  <E>  +  (4ire/e0)[NA  -  fND/(-yRn) ]dx 

=  <E>  -  [2tieNA/(£Q  |  kx | ) ]  Mq  sin|2kxx  +  <J>(z)| 

=  <E>  -  E^  Mq  sin  |  2kxx  +  4>(z)|  (137) 

where 

E  =  Limiting  value  of  the  space  charge  field  (comparable  to 
^  Feinberg's  limiting  value). 

=  2^eNA/( c0 | kx|) 

(10:952) 

Thus,  if  the  charge  transport  lengths  are  much  longer  than  the 
fringe  spacing,  the  space  charge  field  is  sinusoidal  and  of  the  same 
frequency  as  the  interference  fringes  but  shifted  tt/2  from  them.  This 
condition  is  analogous  to  Feinberg's  n/2  shift  for  a  very  large  applied 
external  field. 

Case  2 

If  the  space  charge  screening  length,  Lp,  and  the  length  of  electron 
tightening  by  external  field,  1^.,  are  much  less  than  the  fringe  spacing, 
L: 

Lp  1E  «  L  (138) 

then  the  carrier  concentration  is: 

n  =  <n>  (l  +  M  cos  [2kxx  +  <Kz)]}  (139) 

where 

<n>  =  Average  carrier  concentration  (from  crystal  neutrality) 

=  Nd  (B+  sI0)/(YrNa) 

M  *  Mq  (1  -  &J.  "  %) 

^  -  (Lj/L)2 

^E  *  ^E^ 

and  the  space  charge  field  is  given  by: 
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E  =  (J/e^n)  +  (T/e)  (d  In  n/dx)  +  [pl/(epn)] 

[j/(e|j<n>)]  -  Ej.  M  sin  [2k  x  +  $(z)]  +  E  {1  +  m  cos  [2k  x  +  $(z)]} 

_  U  X  ^  X 

1  +  M  cos  [2k  x  +  <J>(z)] 

X  (140) 

where 


Ep  =  Diffusion  field  (comparable  to  Feinberg's  diffusion  field 

=  2k  T/e 
x 

E  =  Photovoltaic  field 
P 

=  pl0/(e<n>) 

(10:953) 

The  average  field,  <E>,  is  found  by  averaging  over  one  period  after 
making  the  variable  substitution,  £  =  2kxx  +<J>(z): 

<E>  =  1/A  /0A  E  dx 

=  l/(2rr  )  r  [j/(6  p;n>)]  Jq2tT(1  +  M  cos£  )‘!d  ^ 

-  EpM  [(sin  £ )/(l  +  M  cos  £ )]d  £ 

+  Ep  /q2"  [(1  +  m  cos  £;  )/(l  +  M  cos  ^  )]d  £  } 

=  J/[eu<n>(l  +  M2)1/2]  +  Ep  [m/M  +  (M  -  m)/[M(l  -  m2)1/2]  (141) 

The  current,  J,  is  then: 

J  =  eu<n>(l  -  M2)1/2  (<E>  -  Ep  {.m/M  +  (M  -  m)/[M(l  -  m2)1^2]}} 

=  ey<n>(l  -  M2)1/2(<E>  -  Ep0) 

=  o±  Eq  (142) 

where 

Ep°  =  Ep  {  m/M  +  (M  -  m)/[M(l  -  m2)1/2]  } 

=  E 

P 

o  =  Transverse  conductivity 
=  ey<n>  (1  -  M2)1^2 
-  <o>  (1  -  M2)1^2 


<o>  =  Average  conductivity 


Eq  =  Effective  field  causing  macroscopic  current 
=  <E>  -  E 

P 

(10:953) 

Equation  (140)  takes  the  form: 

Eq  (1  -  M2)172  -  Ep  M  sin  (2kxx  +  4>)  +  E  [1  -  m  cos(2kxx  +  $)] 

E  - - - - ^ - - - 

1  +  M  cos(2k  x  +  $) 

x  (143) 

The  Fourier  series  representation  of  this  equation  is  in  the  form: 

E  =  <E>  +  £  E^  exp(i  n  2kxx)  (144) 

where 

n  =  ±  1 ,  2,  ... 

For  which  Kukhtarev  gives  the  Fourier  coefficients  as: 

En  =  [Eq  +  E  (1  -  m/M)  -  iEp]  {[(1  +  M2)172  -  1  ]/M }  n  exp(i  n4>  ) 

=  AE  ([(1  +  M2)172  -  l]At)n  exp[i(n'l)  +  i»]  (145) 

when  written  in  magnitude-phase  form  where 

AE  =  Amplitude  of  electric  field  grating  (comparable  to  Feinberg's 
space  charge  field  (equation  (112))) 

-  {[Eq  +  Ep  (1  -  m/M)]2  +  Eq2}172 

ip  =  Phase  shift  between  the  space  charge  field  and  interference 

pattern  (comparable  to  Feinberg's  phase  shift  (equation  (114))) 

=  tan'1  {Eq/CEq  +  Ep  (1  -  m/M)]] 

(10:953) 

If  the  Fourier  series  can  be  limited  to  first  order,  only,  the 

E-field  will  be  spatially  sinusoidal.  This  can  be  achieved  if  the  term 
2 

[(1  +  M  )  -  l]/M  is  small  so  raising  it  to  any  power  other  than  unity 
adds  an  insignificant  contribution  to  the  series. 

For  small  modulation  conditions,  m,  Mq  and  M  are  all  sma  ’  and: 

limm_>0  [(1  +  M2)  -  1]/M  =  0  (146) 

so  limiting  the  Fourier  series  to  the  first  order  is  a  good 
approximation.  As  the  modulation  index,  m,  increases  toward  unity, 
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however,  M  could  also  approach  unity  under  the  right  conditions.  In  this 
case: 

limm_>l  [(1  +  M2)  -  l]/M  =  .41421  (147) 

which  raised  to  the  second,  third  and  fourth  powers  are  .17157,  .07163 
and  .029437,  respectively.  These  are  by  no  means  relatively 
insignificant  contributions. 

It  is  possible  to  have  small  M  even  if  the  modulation  index,  m,  is 
unity,  however.  If  the  ratio  of  dark  and  photo-conductivities, 

0^/0^  =  B/sIq,  is  large,  or  if  the  space  charge  screening  length,  1^,,  or 
the  length  of  electron  tightening  by  external  field,  1£,  (which  are 
still  assumed  here  to  be  much  less  than  the  fringe  spacing,  L)  begin 
approaching  the  fringe  spacing,  M  can  be  reduced. 

Electro-Static  Field  in  Terms  of  Incident  Intensity 

The  first  order  approximation  of  the  Fourier  series  (equation  (144)) 
will  be  considered  appropriate  by  assuming  small  modulation,  M  <<  1, 
exists.  Equation  (145),  in  terms  of  the  incident  intensities,  becomes: 

Ej  =  [EQ  +  Ep  (1  -  m/M)  -  iEp]  {(1  +  M2)1/2  -  l]/M2}  exp(i4>)  2  cos  20 

(I+1I_1)1/2  (1  -  %  -  5e)/I0 

=  [(I+1I_1)1/2/I0](A  +  iB)  exp(io)  (148) 

where 

limm->0  [(1  +  m2)V2  ~  "  V2 

A  =  [Eq  +  Ep  (1  -  m/M)]  (1  -  -  C£)(cos  20)/(l  +  od/op) 

B  =  -Ejj  (cos  20)  (1  -  Sj.  -  CjO/d  +  od/op) 

(10:954) 

Index  Grating 

For  certain  crystal  orientations,  the  index  of  refraction  modulated 
by  the  linear  electro-optic  (Pockel's)  effect  is  given  by  (this  is  by  no 
means  a  general  form,  see  Appendix  A  for  specific  examples): 
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n  =  n  .  .  .  ±  (1/2)  n  .  3  r  .  .  E  (149) 

original  ong  appropriate 

If  the  electric  field,  E,  is  in  the  form  of  equations  (144)  and  (148), 

the  change  in  the  refractive  index  takes  the  form: 

An  =  ±(l/2)n  ■  3  r  |  <E>  +  [ I . , I  1)1/2/Ir.](A  +  iB) 

ong  approp  1  L  +1  -1'  0J  ' 

exp[i(0  +  2kxx)]| 

=  ±(l/2)n  •  3  r  { <E>  +  [(I. .  )1/2/In](A2  +  B2) 

ong  approp  1  LV  +1  -1'  0JV 

cos(2kxx  +  0  +  il) 
(150) 

(10:954) 

Summary 

Without  assuming  it  to  be  the  case,  Kukhtarev  has  shown  the 
space-charge  field  to  be  sinusoidally  varying  and  proportional  to  the 
incident  interference  pattern  for  small  modulation  index  or  charge 
transport  lengths  greater  than  the  fringe  spacing.  Assuming  this  field 
to  be  sinusoidal  may  not  be  correct  for  conditions  other  than  these  and 
must  be  looked  at  on  an  individual  basis. 

For  the  form  of  the  modulated  index  of  refraction  given  in  equation 
(149),  the  index  grating  is  shown  to  be  spatially  sinusoidal  and 
proportional  to  the  interference  fringes,  as  well. 

Summary 

Both  Feinberg  and  Kukhtarev  developed  models  for  the  phot ore tractive 
effect  which,  for  spatially  sinusoidally  varying  incident  intensity  and 
under  correct  conditions  such  as  small  modulation,  concluded  the 
electro-static  field  set  up  in  the  medium  is  also  sinusoidal  with  the 
same  frequency  as  the  interference  fringes  but  shifted  in  phase  from 
them.  This  field  then  modulates  the  index  of  refraction  of  the  medium. 

It  will  be  shown  in  a  later  chapter  that  diffraction  off  this  index 
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grating  produces  the  phase  conjugate  of  one  of  the  input  beams.  The 
usefulness  of  phase  conjugation  by  the  photorefractive  effect  is  that 
conjugation  is  possible  with  lower  intensity  incident  light  than  by 
degenerate  four  wave  mixing. 

In  the  photorefractive  effect,  the  crystal  is  driven  into 
nonlinearity  by  the  space  charge  field  set  up  in  the  medium  rather  than 
by  the  electric  fields  of  the  incident  light  as  in  four  wave  mixing. 
Both  Feinberg  (equation  (126))  and  Kukhtarev  (equation  (151))  show  that 
the  magnitude  of  the  index  grating  is  dependent  on  this  field  and  the 
elements  of  the  linear  electro-optic  tensor,  R.  Table  2  lists  the 
elements  of  the  electro-optic  tensors  of  several  media. 
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Table  II 

Linear  electro-optic  coefficients,  r 


Wavelength 

EJcciro-opnc  Coefficients 

\ 

'a 

Substance 

Svmmcirv 

Um> 

(10  i;  m/V> 

CdTe 

43m 

1  0 

<n  -  4  5 

3  39 

( T)  r4l  -  6  8 

106 

(  f)  rt,  -  6  8 

33  35 

I  T)  r„  -  5  47 

27  95 

IT)  -  5  04 

GaAs 

43m 

09 

~  I  ' 

l  15 

Jl  -  1  43 

3  39 

lT>  r„  -  (-34 

106 

1  7 )  -  1  51 

G»P 

43m 

0  55-1  3 

<n  o,  -  -io 

0633 

(51  (.,  -  -097 

1  15 

(5)  r„  -  -  1  10 

3  39 

(5)  r„  -  -097 

0-ZnS 

43m 

04 

in  r„  -  1  1 

< sphalerite) 

0  5 

»  r  1  („  -  1  HI 

06 

(7*)  r4l  -  2  I 

0  633 

(St  '.,--16 

3  39 

(Si  '41  -  -14 

ZnSe 

43m 

0  548 

in  -  :o 

0633 

ist  -  :o 

10  6 

in  r„  - 

ZnTe 

43m 

0  589 

in  -  4  5i 

0  61ft 

( 7 1  -  4  3’ 

0  633 

(  r )  r41  -  4  04 
(Si  (.,  -43 

0  690 

( n 

3  41 

in  -  43 

106 

<7-10  -  3  9 

(2:230-233) 


Index  of 
Refraction 

it, 

71  V 

( 10'  m/V> 

Dielectric 

Constant" 

r  unl 

n  -  3.84 

103 

(Si  (  -  94 

it  -  3.60 

120 

n  -  2  58 

94 

*  -  3 .53 

82  , 

n  -  360 

51 

(S)  (-133 

n  -  3  43 

58 

<n  /  -  13.3 

n  —  3  3 

45 

n  —  3  3 

54 

r.  -  3  66-3.08 

IS)  (  -  10 

n  -  3  33 

35 

*  -  3  10 

33 

n  -  3.02 

*  -  3  53 

18 

1  7  1  (  -  16 

n  -  3  43 
*  -  3.36 
n  -  3  35 

(Si  i  -  13  5 

it  —  3  66 

If  -  3  60 
it  -  3  39 

35 

in 

iSi 

r  -  9  1 
7-91 

n  -  306 
it  -  3  01 
it  -  399 

108 

in 

(S) 

i  i 

o  © 

it  -  3  93 
it  -  3.70 
it  -  3  70 

83 

Bi  .CieO.„ 

33 

0  666 

in  -  3  33 

71  -  2.54 

<  t 

23 

0  633 

r„  -  5  0 

tt  - :  54 

<2 

CdS 

6mm 

0  589 

1  7 1  -  3  7 

(7,1.4 

7i ,  -  :  5oi 

7ir  -  :  519 

t  T  )  r ,  -  9*5 

0  633 

1  7 )  (. ,  -  1  ft 

i  7)  r  -48 

7i  -  2  460 

71,  -  2.477 

(ft  r,  -  !0  3? 
iS)  t,  -  9  02 

1  15 

<n  -  3 1 
<n  -  3  3 

(  7 1  -  3  0 

(  7)  (  -  6  2 

71  -  2  320 

71,  -  2.336 

<5l  f,  -  9  5? 

3  39 

(71  („  -  3  5 
(7)  r,,  -  3  9 
|7)  -  30 

(  7|  r  -  6  4 

n  .  -  2  276 

7i,  -  2.292 

106 

1 7i  -  2  45 

IT)  -  2  ’5 
(7)  r,,  -  1  7 

(  7  1  r  -  5  2 

n  -  2  226 

7i,  -  2.239 

CdSt 

6mm 

3  39 

(Si  (  ,  -  1  8 

(  7 1  (, ,  -  4  3 

7i  .  -  2  452 

71,-  2  471 

(7*)  r,  -  9?0 
( D  1 1  -  1065 
lS)  r,  -  9  ?? 
(5)  ty  -  1020 

ZnO 

6mm 

063? 

( Si  -  1  4 

ISI  (.,-2  6 

n  -  1  990 

7i  *  2  006 

r,  -  r:  -  8 

*  1  T 

3  39 

1 5 1  (, ,  -  0  % 

ISI  („  -  1  9 

n  -  1  902 

7i,  -  I  916 

«  ZnS 

IwurtzK?) 

6mm 

0  633 

(SI  (,,  -  0  9 
(Si  („  -  1  8 

71  -  2.347 

71,  -  2  360 

(  T  *  r,  ■  f.  •  8 

<S>  e,  “  8  7 

Pb,,.,.Li„ 

(Ti„,Zr„,  iO, 

x  m 

0  546 

-  «  V,,  -  2320 

n..  -  2  55 

(PLZT) 
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Table  II  (continued) 


Substance 

Ssmmetrs 

^  avclengih 

X 

<  um) 

EJccircvopuc  Coefficients" 
rIM 

( 10  :  m  'Vj 

Index  of 

Refraction 

n 

Dielectnt 
#7  #  Constant* 

( 10  m  \  i  f  1 1,, ) 

LilO. 

6 

0  63? 

(Si  r,-.  -  4  l 
(  S  )  V  ■  14 

(Si  rxx  -  6  4 
(Si  rM  -  3  3 

n,  -  1  8830 
n.-l  7?6'T 

Ag. ,  AsS 

3m 

063? 

(Si  n'r  —  7(i 
<  S  »  n  '  f::  -  29 

n,  -  3  019 
#1,  -  2  7?9 

Li\bO. 

(7  -  1230* 

3m 

Ci 

0  633 

<7l  r ,  •.  -  9  6 
t  T  1  r::  -  PS 
(7|  r\\  -  309 
(  7* )  rK  -  32  6 
\T  i  T  -  21  1 

( .S  )  r , :  -  8  6 

(5)  r..  -  ?4 

(Si  r, ,  -  30  8 
(Si  /v  -  28 

#i.  -  2  286 
#i,  -  2  200 

(  7  i  f  -  » .  -  78 
in#,  —  ?': 

(Sir  -i-4? 
f  S  )  i  ,  -  28 

1  15 

(  7")  r;;  -  5  4 
(7")  r  -  |9 

#?.  -  2  229 
#»,  -  2  150 

3  39 

<  T )  r;:  -  3  1 
(  T  I  r  -  18 

(Si  #■„  -  28 
(S)  r::  -  3  1 

1 S)  r,,  -  6  5 
(5)  #■„  -  23 

n.  -  2  136 
n.-  2  07? 

LiTaO, 

3m 

063? 

(7)  -  8  4 

(7)  r, .  -  30  5 

c  r  i  -  -0.2 

(  7  i  '  -  22 

(S>  rlt  -  75 
(Sir..  -  ?? 

(Si  r,  -  20 
(  S  1  r:;  -  1 

#7.-2  176 

#7,-2  180 

<  7"  i  e  -  #.  -  51 

<  r  i  r ,  -  45 

(Si#  -  f.  -  41 

( 5 )  # .  -  43 

3  ?9 

(S»  -  r 

1  S  1  r  -  4  5 
<S)  V  -  15 
(S»  r:;  -  o  3 

#i  -  2  060 

#i,  -  2  065 

AgGaS 

42m 

0  63? 

(  7  i  rt  -40 
<  7  I  V  -  3  0 

#7  -  2  55? 
n,-  2  50* 

CvK-Av04  42m 

iCDA 

0.55 

iT  t  r4  -  14V 
i  P  i  ...»  1 V  2 

#7,  -  1  5": 

#7.  -  !  550 

KHPO,  42m 

(KDPi 

0  546 

(  P  1  r,  -  (. 

1  Pi  “  10  1 

#i  -15115 
#7,-1  4698 

(  T  1  #  -  #  •  -  42 
(  7 1  r,  -  21 

0  633 

(T  )  r,  •  8 

1  T  1  r, .  -  1  1 

#7,-1  50^4 
#7,-1  4669 

(  S  1  f  ,  —  #  •  —  44 
(  S )  # .  -  21 

3  39 

iT)  r,  -  9  ■ 

(  T  1  n  V„,  “  11 

KD  P04  42m 

iKD'Pi 

('  546 

(Pi  v  *  1b  V 
(  P  t  r4  »  V  V 

#7  -1  50"’9 
#7.-1  468? 

(  T )  1 ,  -  50 
(  S  )  #  —  # -  —  58 

0  633 

ill',  •  24  ! 

#7—1  502 
#7,  -  l  462 

(Si  e •  —  48 

(NH.iH  P04  42m 

(ADP 

0  546 

1  r  1  r„  -  2’  If 
(Pi  o,,  -  8  S6 

#7  -  1  5266 
n,  -  14808 

(T )  t  -  f ;  -  56 
<n  #  -,  -  15 

0  633 

(Pi  V  ”  1141 

1  P  1  ",’Pr.  “17  6 

#i.-J  5220 
#7,-1  477? 

( S )  i ,  -  r  -  -  58 
(Si  #x  -  U 

(NH4>D.P04  42m 

(AD-P, 

0633 

iPl  r(  “40 

(Pi  r,  ~  |(i 

#i  -  1  516 
#7,-1  4^5 

BaTiO  4mm 

IT  -  19?  Ki 

0546 

(Pi  r,  -  164(1 

(Pi  -  10V 

(Sir,  -  820 
(Si  r  -  2? 

#7,.  -  :  43n 
#7f  -  :  365 

( 7" )  t,  -  e .  -  3600 
(7)  «,  -  135 

KTa.Nh,  O, 

(KTNi.  a  -  0.1? 

<P  -  40-60®C l 

0  633 

(Pi  r.  -  8000  (P  -  28l 
(Pi  r  -  500(P  -  28i 
(Pi  r„  -  1000  (P  -  16l 
(Pi  r  -  700 (P  -  I6i 

#7,  -  2  318 
#?,  -  2.27*7 
n[  -  2  318 
#7,  -  2.281 

Ba,  ..Sr  ..Nb-O,  4mm 
if  -  19?  Ki 

0  633 

(P  1  r,.  “  6-r 
(Pi  r,  -  1140 

(7  i  r,  -  42 
(Si  r  -  1090 

#7.  -  2  31P 
#7,  -  2  298- 

1  x  -  3400 
(15  MH21 
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V.  Steady  State  Two-Beam  Couplin; 


The  previous  chapter  described  two  theoretical  models  for  the 
photorefractive  effect.  In  this  chapter  and  the  one  following,  various 
figures  of  merit  for  the  photorefractive  effect  are  presented.  These 
figures  of  merit  are  important  because  they  are  a  quantification  of 
those  aspects  of  the  photorefractive  effect  which  play  a  significant 
role  in  optical  phase  conjugation. 

The  first  of  these  figures  of  merit  is  "two  beam  coupling".  It  has 
been  observed  experimentally  that  when  two  beams  write  an  index  grating 
in  a  photorefractive  crystal  and  propagate  through  this  grating,  one  of 
the  writing  beams  may  experience  a  gain  in  intensity  while  the  other 
experiences  a  corresponding  loss. 

Initially,  this  coupling  effect  appears  not  to  have  a  place  in 
phase  conjugation.  However,  when  self-pumped  phase  conjugate  mirrors 
are  considered,  two-beam  coupling  becomes  an  important  issue. 

Self-pumped  phase  conjugation  occurs  when  a  single  beam  is  incident  on  a 
photorefractive  crystal.  This  beam  is  defocused  (13:46-51)  or  "fanned" 
in  the  crystal  and  reflected  back  to  interfere  with  the  original  beam  by 
either  external  mirrors  (14:689-691)  or  internal  reflections  off  the 
crystal  faces  (15:486-488).  These  reflected  beams  b.  :ome  the  pumps  for 
four  wave  mixing  and  get  their  energy  from  the  original  beam  through 
two-beam  coupling  (14:688-691)(15:486-488)(16:548-553). 

Two  theories  of  two  beam  coupling  are  presented.  The  first  is 
developed  from  Feinberg's  Hopping  Model.  The  second  assumes  a  sinusoidal 
index  grating  exists  in  the  medium  due  to  modulation  by  the 
electro-static  field  in  an  "index  ellipsoid"-type  argument.  (See 
Appendix  A  for  details  on  the  index  ellipsoid.)  Both  theories  show  that 


f 
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as  long  as  there  is  a  phase  shift  between  the  interference  fringes  and 
index  grating,  4»,  increased  intensity  can  occur  for  one  of  the  writing 
beams,  with  decreased  intensity  of  the  other,  provided  the  crystal 
orientation  is  correct. 


Feinberg's  Hopping  Model 

This  derivation  of  an  expression  for  two  beam  coupling  is  a 
continuation  of  the  development  of  the  hopping  model  in  the  previous 
chapter.  Assumption  6  made  in  developing  the  model  was  that  small 
modulation  conditions  would  only  be  considered  (m  <<  1).  The 
implication  of  this  is  that  one  writing  beam  is  much  weaker  than  the 
other.  For  this  development,  then: 

Beam  2  >>  Beam  1 

The  nonlinear  polarization  of  the  medium  due  to  the  modulated 
susceptibility  (equation  (153))  and  the  superposition  of  the  two 
electric  fields  is: 


pnl  =  x*e 


(151) 


=  X  -[(E1/2)e1  +  (E2/2)e2] 

=  X-e2  E2/2 

where,  from  equations  (110)  and  (126),  the  susceptibility,  X,  is: 

W0  P  q  2  Ej  E2"  a,-e2*  a2  *  iaf 

0w  w  e  eQ  k2  lEji2  +  |E2|^  a  +  1  +  iaf  (152) 

(5:1300). 

From  Maxwell's  equations,  with  o  =  0  and  the  SVEA  similar  to 
previous  development,  the  wave  equation  is  in  the  form: 
k^x  3E^/3x  +  kjy  3Ej/3y  +  k^  3EL^/3z 

+  i|kj|2  E^/2  -  i  Mq  e  w2  Ej/2  =  (i  o2  E^A)  X’e2 


exp[i(k1-k2+k2-k1)*x]  (153) 
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(154) 


If  there  is  no  driving  term  (the  nonlinear  polarization)  then: 


klx  aE^/ax  +  kly  aE^/ay  +  klz  aEj/az  =  o 
But  with  the  driving  term,  the  modulated  susceptibility,  X,  alters  the 
dispersion  relation  for  k^  as  (5:1300): 

=  kj(u)  +  6kj 


kj2  =  k^ta)  +  2  ^(^•fikj  + 


=  (oj )  +  2  kj(u))#6kj  (1 

if  6k^  <<  ^(w). 

Two  assumptions  are  made  in  order  to  achieve  a  solution: 

1.  Assume  that  the  changes  in  k^,  6k^ ,  are  much  smaller  than  the 
unperturbed  portion,  k^(u)  (already  used  in  equation  (155))  (5:1300): 


(155) 


6^  «  k^u)  (156) 

2.  Assume  k^x  aEL^/ax  +  k^y  aE^/ay  +  k^z  aE^/az  =  0  as  in 
non-driven  case  (equation  (154)). 

Equation  (153)  then  becomes: 

il^2  (u)  +  2  k^(u>)  *6k^j  E^/2  -  i  cm2  E^/2 

=  (i  P()  w2  E^/4)  (X  .%2)  (157) 

2 

Using  pq  cq  =  1/c  and  expanding  X  where  |  E?|  >>|  EJ  so  that 
|E2j2/(|E1|2  +  |E2|2)  =  1,  then  (5:1300): 

kj-fikj  =  [m2/(4  c2Cq)]  (E2/Ex)  ^*.X-fc2 
i  u)2  fn  a  +  if 

2  c2  a2  +  1  +  iaf  1  w  u  2  1  2  (158) 


If  k^  is  approximately  parallel  to6k^  then: 


i  u  2  f  n  (l+a2)a+af2  +  if  ...  ... 

— - -  - - - - —  (e1  •  e  *R>fce  •/bo)C^i,/fe?  ) 

„2  i  .  i  /i  .  Z\4  .  /_  f\Z  1  to  to  z  1  Z/. 


A  2  c  |  kjJ  (1  +a  )  +  (af)  u  w  (159) 

and  the  propagation  vector  alteration,  6k^,  causes  either  growth  or 
decay  in  E^  described  as: 
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E^Hx)  =  E1  cxp[i(k1  +  dk^x] 

=  Ej  exp[  i(ki  +  Re  6kj)*x]  exp(-Im  <$kj*x)  (160) 

The  intensity  becomes  (5:1300,1301): 


I1'(x)  =  |E^p  exp(-2  Im  <$k^*x) 

I^(x)  +  AI^(x)  =  |EjJ^  (1  -  2  Im  fik^’x)  (161) 

using  the  first  order  series  expansion  for  the  exponent. 

The  fractional  intensity  gain  for  E^  over  some  interaction 
distance,  L,  can  be  written  as  (5:1302): 

Al^(x)/I^(x)  =  -2  Im  6k^*r  L 

-  fp  L  a(l  +  +  f^) 

"  —  r” . ;:“"rr:'r;:r  <V,vE‘k‘Vs2)(v52") 

c  k1  (1  +  a  )  +  (q-f) 

1  (162) 

Note  that  from  equation  (114),  if  there  is  no  phase  shift  between 
the  interference  pattern  and  the  index  grating,  ip ,  the  numerator  of 
equation  (162)  is  also  zero.  Thus,  there  is  no  energy  transfer  between 
the  beams  under  this  condition. 

A  comparison  of  the  functional  dependence  of  equation  (162)  to 
experimental  data  is  made  in  the  applications  chapter  (VII). 

Examples 

The  following  are  examples  of  the  form  the  relative  intensity 
gain,  equation  (162),  would  take  for  two  specific  crystals,  BaTiO^  and 


BaTiO^ 

(®l*’VR,k'V82KV82*)  *  t(n04  r13  kz  e2x  v  n02  ne2  r51  kx  e2z)2 

.  /  4  ,  .22  .  x2 

+  (n  r,-,ke  +n  n  r,-,  k  e9  ) 

o  13  z  zy  o  e  51  y  2z 

2  2  2  2 
+  (n  n  rc,  k  e-,  +  n  n  r,-,  k  e0 

o  e  51  x  2x  o  e  51  y  2y 

4  ,2i 

+  ne  r33  kz  e2z>  1 

<elx  e2x  +  ely  e2y  +  elz  e2z)  (163) 
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where 


e^  is  a  mathematical  expression  used  here  to  clear  dot  product 
of  directionality  and  leave  only  its  magnitude.  This  is  done  because  it 
was  assumed  the  stimulus  (X*^)  and  the  response  (e^ ' )  were  of  the  same 
polarization. 

e^  =  Unit  polarization  vector  of  input  wave, 

=  e^x  +  e^  y  +  eiz z  (see  Figure  9  for  orientation  of  directional 
cosines) 


=  Unit  polarization  vector  of  input  wave,  E2 

=  e0  x  +  e0  v  +  e0  z 
2x  2y-  2z 

e  =  Dielectric  tensor  for  BaTiO-, 

(J  J 


R  =  Contracted  electro-optic  coefficient  tensor  for  BaTiO^ 


k  =  Direction  vector  of  grating 

=  \  ~  *2 

=  kx  +  ky  +  kz 
x  y7  z 

k^  =  Propagation  vector  of  input  E^  wave 

*  k.  x  +  k.  y  +  k.  z 
lx  ly7  lz 

k.  =  k,  cos  a.  cos  0, 
lx  1  1  1 

kly  =  cos  sin 
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k,  =  -  k.  sin  a, 
lz  1  1 

k2  =  Propagation  vector  of  input  E,  wave 


=  k2xx  +  k2yy  +  k2zz 
k2x  =  k2  cos  a  2  cos  02 


k2y  =  k2  cos  a 2  sin  02 
k2z  =  k2  sin  a  2 


a)  Vertical  Polarization  b)  Horizontal  Polarization 
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BS0(Bi12Si02Q) 


c 

OJ 


'R’k* 


e 

u> 


n  c41^kz  e2y  +  ky  fc2z' 

+  <kz  e2x  +  kx  e2z)2 

+  (k  e~  +  k  e~  )2}1^2 
y  2x  x  2y  J 

(elxe2x  +  elye2y  +  Hz^  (164) 


where 

c  =  Dielectric  tensor  for  BSO 

u) 

n2  0  0  ^ 

=  0  n2  0 

\0  0  u2f 

R  =  Contracted  electro-optic  coefficient  tensor  for  BSO 

'o  0  0  \ 

0  0  0 

0  0  0 

r41  0  0 

0  r4l  0 
\°  0  r4l/ 


Two  Beam  Coupling,  Assuming,  a  Sinusoidal  Index  Grating  Exists 

Many  authors  such  as  Kukhtarev  (10:954-956),  Vahey  (12:3510-3515) 
and  White  (17:10)  assume  an  index  grating  has  been  set  up  in  a 
photorefractive  medium  and  derive  expressions  for  the  growth  or  decay  of 
the  writing  beams  propagating  through  this  grating.  Basically,  each 
development  is  similar  with  each  author  including  his  own  details  in  the 
picture.  The  derivation  done  here  closely  follows  Vahey's  derivation 
but  is  aimed  at  attaining  White's  solution  which  is  more  readily 
applicable. 


Derivation 


Given  a  sinusoidal  index  grating  exists  in  the  medium  as  derived  in 
the  previous  chapter,  in  the  correct  crystal  configuration,  that  grating 
may  be  represented  as  in  equation  (149)  where  the  E-field  is  given  by 
Feinberg's  space-charge  field  (equation  (115))  (7:519, 520)(17:10): 

Appropriate  ”  ^^2^  napprop  rapprop  m  ^ 

{(Ed2  +  Eq2)/[(Ed  +  Eq)2  +  Eq2]}  1/2  cos(k* r  +  40 


=  n^  +  m  cos  (k*r  +  40 


where 


nT  =  ±(l/2)n  r  E  ( (En  +  Enz)/[(En  +  E  )Z  +  E^]} 

I  approp  approp  q  u  0  T)  q'  0  J 

(see  Appendix  A  for  definition  of  ± ) 
m  =  2  (I1l2)1/2/(I1  +  I2) 

4  =  tan‘1[(ED  Eq  +  ED2  +  E02)/(E0  E  )]  (equation  (114)) 

(The  geometry  follows  Figure  7.) 

Wave  Equation 

The  scalar  wave  equation  for  this  geometry  is: 

V2E  +  k2E  =  0 

where 
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2 i i 1/2 


(166) 


E  =  A^ej  exp(ik^’r)  +  A^'e^"  exp(-ik^'r)  +  A2e2  exP 

+  A2"e2“  exp(-ik2*r) 

=  A^(x)  exp[i(k  cos  0x  -  k  sin  0z)]  +  c.c. 

+  A2(x)  exp  [i(k  cos  0x  +  k  sin  ©z)]  +  c.c. 

c.c.  =  Complex  conjugate 

,2  2  2„,w  2 
k  =  w  n  (E)/c 

2  2  2  2  2 
n  (E)  =  nQ  +  2nQ  n^  m  cos  (k*r  +  4*)  +  n^  m  cos  (k*r  +  40 

^  2 

■  n^  +  nQ  rij  m  {exp[i(2k  sin  0z  +  4  )] 

+  exp[-i(2k  sin  0z  +  4)]} 

After  applying  the  Slowly  Varying  Envelope  Approximation  (SVEA), 
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equation  (166)  becomes: 

2(3A^/3x)(ik  cos  0)  exp[i(k  cos  0x  -  k  sin  0z)]  +  c.c. 

+  A^(-kz  cosz  0)  exp  [i(k  cos  0  x  -  k  sin  0  z)]  +  c.c. 

+  2(3A2/9x)(i  k  cos  0)  exp[i(k  cos  0  x  +  k  sin  0  z)]  +  c.c. 

+  A2(-k  cosz  0)  exp[i(k  cos  0  x  +  k  sin  0  z)J  +  c.c. 

-k2  sin20  E  +  (u2  nQ2/c2)  E 

+  (w2  n^  nj.  m/c2){A^  exp[i(k  cos  0  x  +  k  sin  0  z  +  40]  +  c.c. 

+  Aj  exp[ i(k  cos  0  x  -  3  k  sin  0  z  -  4')]  +  c.c. 

+  A 2  exp[i(k  cos  0  x  +  3  k  sin  0  z  +  40]  +  c.c. 

+  A?  exp[i(k  cos  0  x  -  k  sin  0  z  -  <!']  +  c.c.}  =  0 

(167) 


Grouping  terms  in  equation  (167)  synchronous  with 
exp[i(k  cos  G  x  -  k  sin  0  z)]  yields: 

2(3A^/3x)(i  k  cos  0)  +  A^-k2  cos2  0)  -  k2  sin2  QA^  +  to2  n^2  A^/c2 

2  2 
+  u)  nA  nT  m  A0  exp(-i40/c  =  0 

0  1  2  (168) 

Since: 

u2  n^2  A^/c2  -  k2(sin20  +  cos2  Q)Aj  =  0  (169) 

after  rewriting  ,ie  exponential,  equation  (169)  becomes: 

3Aj/3x  -  i  k  nj  m  A2  exp(-i40/(2nQ  cos  0)  =  0 
3A^/3x  -  k  n^  Hi  A^  sin  <\i/(2uq  cos  0) 

-  i  k  nj  m  A2  cos  4>/(2nQ  cos  0)  =  0  (170) 

The  real  part  of  equation  (170)  yields: 

SA^Sx  =  knx  A2  (sin  (I’i  +  i2)  cos 

=  (l/A1)(dI1/dx)  (171) 

by  definition.  Then: 

dlj/dx  =  k  nj  (sin  40/[nQ  (1^  +  ^)  cos  0] 

dlj/dr  =  2  ri1I2/(I1  +  I2)  (172) 

where 
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r  =  x/cos  0 
T  =  it  rij  sin  1|,/X 

A  similar  development  for  exp[i(k  cos  0  x  +  k  sin  0  z)]  yields: 

dl2/dr  =  -  2  ri^/Clj  +  I2)  (173) 

White  (17:10)  points  out  that  linear  absorption,  a,  is  crucial  to 
the  photorefractive  effect  so  his  coupled  equations  are  in  the  form: 
dlj/dr  =  -alj  +  2  ri  I2/(l  +  l2) 

dl2/dr  =  -al2  -  2  rI1I2/(I1  +  I2)  (174) 

Solution  to  the  Coupled  Wave  Equations 

White  gives  the  solutions  to  equations  (174)  as  (17:10)  (10:954): 

Il(r)  =  IjCO)  exp(-ar)[l1(0)  +  l2(0)]/[l1(0)  +  I2(0)  exp(-2rr)] 

I?(r)  =  I?(0)  exp(-ar)[l1 (0)  +  l2(0)]/[l1 (0)  exp(2rr)  +  1,(0)] 

12  1  2  (175) 

Note  that  the  term  responsible  for  the  gain  or  loss  of  intensity,  r, 
contains  the  change  in  the  index  of  refraction,  n^,  which  is 
proportional  to  the  magnitude  of  the  space-charge  field  and  the  phase 
shift  between  the  incident  interference  fringes  and  the  space-charge 
field,  4>.  If  there  is  no  phase  shift,  4*  =  0,  only  the  linear  absorption 
loss  remains  in  equation  (175)  so  no  energy  transfer  between  the  two 
beams  occurs. 

Sumnary 

'*Two  beam  coupling"  is  one  of  the  figures  of  merit  for  the 
photorefractive  effect  as  applied  to  optical  phase  conjugation.  It 
describes  the  phenomenon  of  intensity  growth  (or  corresponding  loss)  in 
one  of  the  incident  beams  writing  an  index  grating  in  a  photorefractive 
crystal. 

Two  theoretical  expressions  were  developed  to  describe  two  beam 
coupling.  The  first  (equation  (162))  is  an  extension  of  Feinberg's 
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Hopping  Model  for  the  photorefractive  effect  developed  in  the  previous 
chapter.  The  second  (equation  (175))  assumed  an  index  grating  existed 
and  could  be  described  by  an  index  ellipsoid- type  expression.  It  solved 
the  wave  equation  for  light  propagating  through  the  grating. 

Although  expressions  (162)  and  (175)  are  unique,  they  both  indicate 
that  if  no  spatial  phase  shift  exists  between  the  interference  fringes 
and  the  index  grating,  no  beam  coupling  will  occur.  A  comparison  of 
these  two  expressions  occurs  in  the  applications  chapter  (VII)  for 
specific  crystal  orientations  of  Bismuth  Silicate  (BSO). 


VI.  Diffraction  Efficiency  and  Phase  Conjugate  Reflectivit 


The  second  and  third  figures  of  merit  quantifying  the  aspects  of  the 
photorefractive  effect  important  to  optical  phase  conjugation  are 
diffraction  efficiency  and  phase  conjugate  reflectivity.  Diffraction 
efficiency  is  the  amount  of  power  in  a  wave  diffracted  off  a  grating  at 
the  Bragg  angle  compared  to  the  amount  of  power  in  the  wave  incident  on 
the  grating.  Phase  conjugate  reflectivity,  in  the  photorefractive  phase 
conjugation  geometry  shown  in  Figure  10,  is  the  power  in  the  conjugate 
wave  compared  to  the  power  in  the  probe  wave.  The  probe  wave  and  one  of 
tne  pump  waves,  E^,  write  an  index  grating  in  the  crystal.  The 


counterpropagating  pump  wave,  E^,  diffracts  off  the  grating  as  E^,  the 
phase  conjugate  of  the  probe. 


Figure  10.  Geometry  of  optical  phase  conjugation  in  photorefractive 
crystals. 


As  for  two  beam  coupling,  two  theories  for  both  diffraction 
efficiency  and  phase  conjugate  reflectivity  are  presented.  The  first  for 
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each  is  again  an  extension  of  Feinberg's  Hopping  Model.  The  second  for 
each  again  assumes  an  index  grating  exists  in  the  medium  due  to 
modulation  by  the  electro-static  field. 


Feinberg's  Hopping  Model 

Development 

The  interference  of  two  optical  waves  in  a  photorefractive  crystal 
(E^  and  E-,  in  Figure  10)  causes  charge  migration  and  the  steady  state 
potential  distribution  as  described  in  Chapter  IV.  When  a  third  optical 
wave  (5:1300) 

E^  =  e^  Re  exp[i(kyx  -  wt)]  (176) 

is  scattered  off  the  modulated  optical  susceptibility  set  up  by  the 
potential  distribution  and  E^  is  counterpropagating  to  E^,  the  scattered 
wave,  E^,  will  be  the  phase  conjugate  of  E2. 

The  Wave  Equation 

The  solution  for  the  four  wave  mixing  will  be  derived  through 
Maxwell's  wave  equation  with  the  expression  for  the  nonlinear 
polarization  (nonlinear  susceptibility)  given  by  (equation(151)).  The 
wave  equation  is  (equation  (15)): 

v2  e  =  u0  a  a  E/at  +  m0  e  a2E/at2  +  pQ 

As  in  the  second  and  third  order  mixing  derivations,  a  number  of 
simplifying  assumptions  will  be  made  to  make  the  equation  soluble: 

1.  The  media  is  non-conducting  (c  =0).  Then: 

"2  exp[i(kz>*x  -  ut)]  -  p^  e(32/3t2)E4  exp[i(k^*x  -  ut)] 

=  p0(  32/3t2)P3(x,t)  (177) 

where: 

(32/3x2)E^(x,t)  =  (3 2E/f/  x2)exp[i(k4' x  -  ut)] 

+  i  2k^x  (  3E4/3t)exp[i(k4*  x  -  ut)] 

-  k/+x2  E4  exp[i(k4*x  -  ut)]  (178) 
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Similarly  for  y  and  z: 


(32/3t2)E4(x,t)  =  -CO2  E4(x,t) 

(32/3t2)P3(xft)  =  -to2  P3(x,t)  (179) 

and  the  nonlinear  polarization  is  in  the  form  (5:1300): 

P3(x,t)  =  X(x)*E3(x,t) 

=  E^/2  Re  X*e3  exp[ i[ (k^-lc^+k^) ’X  -  cot]]  (180) 

2.  Assume  the  Slowly  Varying  Envelope  Approximation  (SVEA): 

7  7  7  7  7  7 

3  /3x  =  3  /3y  =  3  /3z  =  0 


The  wave  equation  becomes: 


2i  k4x  (  3Ea/ 3x)exp[i(kA  -x  -  tot)]  -  kAvZ  EA  exp[i(kA*x  -  cot)] 


4x  4 
2 


+  21  k4y  (3E4/ay)exp[i(k4*x  -  cot)]  -  k4y  E4  exp[i(k4*x  -  cot)] 

2 

+  2i  k4z  ( 3EA/3z)exp[i(kA *x  -  tot)]  -  kA?  EA  exp[i(kA*x  -  cot)] 


4z 


+  Uq  e  co  E4  exp[i(k4*x  -  cot)] 


=  (- 


^0 


CO 


'/2)E7  X(x)  *e,  exp[i[  ( k-  -k~+kO  *x  -  cot]] 
3  3  1  2  3  (lgl) 


By  multiplying  both  sides  by  (-i/2)  exp[-i(k4*x  -  cot)],  using  the 
2  2 

identity,  k4  =  Uq  e  co  and  if  =  u,  then: 

k  3E./3x  +  k,  3E,/3y  +  k,  3E,/3z 
4X  4  4y  4  7  4z  4 

=  (i  Uq  co2/4)E3  X(x)*e3  exp[iAk*x)  (182) 

where: 


Ak=k3-k2  +  k3~k4 


h  Ihe  scattered  beam,  E-,,  is  not  depleted  by  the  mixing,  then: 

E3(x)  =  E^ 

Solution 

Ihe  form  of  the  solution  is  assumed  to  be: 

E4  =  K  exp(iAk’x)  +  C 
■  :  i s  the  amplitude  of  E4(x),  only.) 


(183) 


(184) 


(185) 
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(186) 


3 E^/9 x  =  i  Akx  K  exp(iAk*x) 

with  similar  forms  for  y  and  z.  Equation  (182)  becomes: 

i  K  Ak»k^  exp(iAk-x)  =  (i  p^  <j /4)E3  X*t:3  exp(iAk*x)  (187) 

Therefore: 

K  =  uq  u  E3  X*e3/(4A  k-k4)  (188) 

To  solve  for  the  constant,  C,  if  there  is  no  interaction  length  in 
the  medium,  no  nonlinear  polarization  occurs  and  no  new  wave,  E^,  is 
produced : 

E4(0,0,0)  =  0 

=  K  +  C 

C  =  -K  (189) 


The  conjugate  wave  which  is  the  fourth  wave  of  the  four  wave  mixing 


process  can  then  be  expressed  (after  minor  manipulation)  as: 

Pq  u>^  E^  X*e3  exp(iAk*x/2)  -  exp(-iAk*x/2) 


E4  = 


4  Ak*k, 


exp(-iAk-x/2) 


(190) 


or,  over  interaction  length,  L,  as: 


e4(l)  = 


iL  pQ  (*)  E3  X 


sin  (Ak  L/2) 

4  k4  exp(-i  Ak  L/2)  Ak  L/2 
The  power  per  unit  area  radiated  by  the  nonlinear  polarization 


(191) 


takes  the  form: 
P^/Area  =  1/2  (c/pq) 


1/2 


^4' 


1 


c  a/2 


L^  pq^  J*  |E3  |e4" *X*e3 sin  ( Ak  L/2)  3 

- - . .  [ - ] 

16  k, 


Ak  L/2 


(192) 


=  -  (?-)J 

2  Pq  *v'  ^4 

Under  conservation  of  momentum,  Ak  =  0  so  that  (sin  x)/x  =  1  in  the 
preceding  expression  for  any  significant  power  transfer. 

Diffraction  Efficiency 

The  power  per  unit  area  of  the  pump  wave,  E^,  is: 

►"O'-'  "  it,-, i-  (193) 


P^/Area  =  (l/2)( e/pn)l/2  |E^|2 
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The  ratio  of  the  power  radiated  by  the  nonlinear  polarization  to  the 


input  power  of  is  the  diffraction  efficiency,  q,  and  from  equations 
(192)  and  (193)  is  given  by: 

n  =  |L  m0  w2  e^"-X'e^/(4  k4)|2 


L  CO  f  r 


2  ,  f2 
a  +  i 


=  ( _ ^o_)2 _ . 

4  n4  c  k  (Ij  +  I  2  +  I3)2  (a2  +  l)2  +  (af) 


x  I  e,  •  g  •  R*k*e 

Z  1  4  gj  a)  J1 


Ixj 


2I1X2 


2  2 

_  ,  N-C  ±  Z  r.2  ED+E0  ,  -s  *  -  I  2 

_  ( - )  - -  £q  - - - -  |e  .E  .R.jj.e  .e, 

4  n4  c  k  (It  +  I2  +  I 3)  (Ep  +  Eq)2  +  ZQZ  4  w 


cj  '3  * 
(194) 


using  Uq  cj/k4  =  l/(n4  c  Cq)  and  where  the  modulation  index  is: 

rn  =  2(I1I2)1/2/(l1  +  I2  +  l3)  (195) 

Note  that  Feinberg  includes  I3,  the  third  incident  intensity,  in  the 
denominator  of  this  form  of  the  modulation  index  expression  which  does 
not  appear  in  the  usual  forms  (5:1299,1300). 

Feinberg's  form  of  diffraction  efficiency,  then,  is  proportional  to 
the  square  of  the  magnitude  of  the  created  space-charge  field. 

Phase  Conjugate  Reflectivity 

The  power  per  unit  area  of  the  probe  wave,  is: 

P(2)/Area  =  (l/2)(e /uQ)1/2  |  E2|  2  (196) 

The  ratio  of  the  power  radiated  by  the  nonlinear  polarization  to  the 
input  power  of  E2  is  the  phase  conjugate  reflectivity,  R,  and  from 
equations  (192)  and  (196)  is  given  by: 


R  = 


EJ  2/|  E2|  2 


=  |  L  oj  e4"*X-e3/(4  n4  c  eQ> ]  (Iy^) 


( 


Lu  f 


0)2 


2I1I3 


2  +  f2 
a  +  1 


*  — * - * - *  I  “fez  »e  •R»k»e  ’ed 

4  n4  c  k  (Ij  +  I2  +  I3)2  (a  +  l)1  +  (afr  4  u  o  J 


L  u 


2i1T3 


2  2 
V  +  Eo 


=  ( . )  - - ^  Eq 

4  n4  c  k  (Ij  ♦  I2  +  I3)Z  (Ep  +  E  y  +  Eq 


- - -  e  .c  ♦  R.  fee  •  eJ 

Z  ,  ~  Z  1  4  U)  Ll)  31 


(197) 
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Feinberg's  form  of  the  phase  conjugate  reflectivity,  like  the 
diffraction  efficiency,  is  proportional  to  the  square  of  the  magnitude 
of  the  space-charge  field. 

Examples 

Ihe  following  are  examples  of  the  form  the  diffraction  efficiency 
and  phase  conjugate  reflectivity  expressions  (equations  (194)  and  (197)) 
would  take  for  two  specific  crystals,  BaTiO^  and  BSO. 

BaTiO-, 


a  ju  2  A  2  2  2 

|e/'#e  *R*k*e  *e,|  =  (n  r., ,  k  e,  +  n  n  r,-.  k  e,  ) 

'4  to  to  3  ol3z3x  o  e  51  x  3z 


A  2  2  2 

+  (n  r,  k  e,  +  n  n  r..  k  e0  ) 

o  13  z  3y  o  e  51  y  3  z' 

+  (n  n  r,-,  k  e0 

o  e  51  x  3x 

2  2 

+  n  n  r,-.,  k  e0 
o  e  51  y  3y 


+  n  ^  r00  k  e0  ) 
e  33  z  3z/ 


(198) 


where 

e^',  as  used  here,  is  a  mathematical  expression,  as  in  the  two 
beam  coupling  examples,  to  clear  the  dot  product  of  directionality  and 
leave  only  its  magnitude.  It  was  assumed,  here,  that  the  stimulus 
(X*e^)  and  the  response  (e^)  were  of  the  same  polarization. 

e^  =  Unit  polarization  vector  of  pump  wave, 

=  e,  x  +  e^  y  +  e,  z  (see  Figure  11  for  orientation  of 
3X  ^  z  directional  cosines) 

k  =  Direction  vector  of  grating 

e  =  Dielectric  tensor  for  BaTi0o 
w  3 

R  =  Contracted  electro-optic  coefficient  tensor  for  BaTiO^ 

See  two  beam  coupling  examples  for  details  on  directional  cosines, 
the  dielectric  tensor  and  the  electro-optic  coefficient  tensor  for 
BaTiOy 
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a)  Vertical  polarization 


b)  Horizontal  polarization 


e3x  Sln  “l 


e3x  =  s:''n  ® 


e.,  =0 

3y 

e0  =  cos  a. 
3z  1 


e^  =  cos  0 
e3z  ‘  0 


Figure  11.  Orientation  of  angles  used  in  directional  cosines  of 
polarization  vector,  e^. 


BSO  (Bi^ 2Si02p) 


’e  'R*k*e  • 

CJ  OJ 


®3l2  *  "8  r412((kz  e3y  +  ky  e3z)2  +  (kz  e3x  +  kx  e3-/ 

+  <ky  e3x  +  kx  e3y)2] 
(i99: 


where 

e  =  Dielectric  tensor  for  BSO 
w 

R  =  Contracted  electro-optic  coefficient  tensor  for  BSO 
See  two  beam  coupling  examples  for  details  on  the  dielectric  tensor  and 
the  electro-optic  coefficient  tensor  for  BSO. 


Diffraction  Efficiency  Assigning  a  Sinusoidal  Index  Grating  Exists 
This  derivation  is  aimed  at  attaining  the  form  of  diffraction 
efficiency  used  by  Kukhtarev  (10:955)  and  Gunter  (9:250).  It  is  assumed 
that  a  sinusoidal  index  grating  exists  in  the  medium  as  derived  in 


0 

i 


Chapter  IV,  the  geometry  follows  Figure  7  and  the  crystal  is  oriented  so 
equation  (165)  is  applicable  for  the  modulated  index.  Then  using  the 
wave  equation  (166),  the  same  derivation  as  in  the  previous  chapter  can 
be  made  but  with  electric  fields  replaced  by  primed  electric  fields 
representing  beams  of  the  same  frequency  as  the  writing  beams  incident 
at  the  Bragg  angle  (10:955)  but  with  small  enough  intensity  that  they  do 
not  effect  the  index  grating  (9:250).  Then  equations  (172)  and  (173)  can 
be  written  as: 

dl^'/dx  =  2  rij'Ij'/dj*  +  I2') 

dl2'/dx  =  -  2  rij'Ij'/dj'  +  I2’)  (200) 

where 


r  =  tt  n.  cos  e/x 

II  ^  =  Intensity  of  beams  of  the  same  frequency  as  the  writing 
beams,  1^  2,  incident  at  the  Bragg  angle  but  with  small 
enough  intensity  that  they  do  not  effect  the  index  grating 
But  n,  is  as  defined  in  equation  (186)  as  in  Kukhtarev's  development 
of  Chapter  IV  (10:954). 

The  solution  to  the  coupled  equations  (200)  for  diffraction 
efficiency,  n >  with  boundary  condition  I2'(0)  =  0  since  the  diffraction 
of  Ij'  is  I2',  is  given  by  Kukhtarev  as: 
n  =  I^W/I^CO) 

2  8q  exp(r'  x/2)[cosh  (r'x/2)  -  cos(yx)] 

(1  +  Bq)  [1  +  Bq  exp(r'x)]  (201) 

where 


B0  =  I2(0)/I1(0) 

l'  =  (— 2rr/X )  n  .  ^  r  cos  0  B 

x  7  ong  approp 

3 

y  =  (~n/x)  n  .  r  cos  0  A 

'  x  '  ong  approp 

with  A  and  B  as  defined  after  equation  (148)  of  Chapter  IV  (10:955). 
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(Note  the  intensities  in  8q  are  of  the  writing  beams.) 

In  this  form,  the  diffraction  efficiency  is  proportional  to  the 
exponential  of  the  diffusion  field,  only. 

Phase  Conjugate  Reflectivity  Assuming  a  Sinusoidal  Index  Grating  Exists 
The  optical  phase  conjugation  geometry  of  Figure  10  is  assumed.  At 
this  point,  no  index  grating  exists  in  the  medium  but  the  probe  and  pump 
(Ej)  will  write  one.  The  scalar  wave  equation  (166)  is  : 

V2E  +  k2E  =  0 

where  E  is  the  superposition  of  the  four  electric  fields: 

E  =  Aj(r)  exp[i(k^  cos  0^x  -  wt)j  +  c.c 
+  A2(r)  exp[i(k2  cos  Q 2^  -  ut)]  +  c.c. 

+  A^(r)  exp[-i(k^  cos  O^x  +ut)]  +  c.c 

+  A^(r)  exp[-i(k2  cos  +  ut)]  +  c.c.  (202) 

and  where 

3  E  /3  y  =  3  E/9  z  =  0 

.2  2  2^rw  2 

k  =  u  n  (E)/c 

Here,  n(E)  takes  the  same  form  as  equation  (lbr>),  however,  initially, 
each  pair  of  beams  in  Figure  1<>  will  !*■  allowed  to  interfere  with  one 
another  and  create  a  grating.  Equation  (lb5)  then  takes  the  form 
(7:519): 


no 

+  [n 

I  exp( i  d  j )/2  ]  (Aj ' 

A,  +  A  >  exp(  i  k 

[T)/I0  + 

+ 

[nII 

exp( i  v  ^ )/2  ]  (AjA^ 

♦  A  /'A,  >  expUkjj* 

r)/l„  +  c.e. 

+ 

[nm 

J 

exp(  i  j)/2  ]  Aj A-j 

exp(  i kj  |  j*  r)/I0  + 

C.C. 

+ 

[niv 

exp(i  * ,)/2  J  A^"A2 

exp(  ikjV*  r  )/I()  +  c 

.c . 

no 

+  A  n 

(203) 

where 


nI, II, III, IV  napprop  rapprop  +  E0  +  V  +  E0  ^ 

(using  Feinberg's  space-charge  field.) 

4*1  ii  jy  =  Appropriate  phase  shift  between  that  set  of 

interference  fringes  and  index  grating 

kj  =  k2  -  k^ 

kII  =  *1  +  ^ 
kIII  =  2kl 
kIV  =  2k2 

X0  =  X1  +  l2  +  13  +  rA 

If  it  is  assumed  that  the  change  in  the  index  is  small  compared  to 
the  background  index,  An  <<  nQ,  then: 

n2(E)  =  n^2  +  2  Oq  An  +  An2 

=  nQ2  +  2  nQ  An  (204) 

By  using  equations  (202),  (203)  and  (204)  in  the  wave  equation  and 
applying  the  Slowly  Varying  Envelope  Approximation  (SVEA),  the  terms 
synchronous  with  exp[i(k^  cos  0j  x  -  ot)]  in  the  wave  equation  are: 

2kj  cos  Gj  dAj/dx  =  [i  u  npn^  exp(-i  4^)  (A-p15^"  +  A^  A^)  A2 

+  i  u  nQ  nn  exp(i  4n)  +  A3  A2'  \ 

+  i  w  nQ  nni  exp(i  4m)  (AjA^  )  A^]/(c  Iq) 

(205) 

Now  it  is  assumed  that  although  all  these  gratings  are  created,  only 
one  will  give  rise  to  strong  beam  coupling.  This  grating  is  created  by 
the  interference  of  beams  1  and  2,  and,  3  and  4  and  characterized  by 
n^  exp(i4  ^)  (7:520),  then  equation  (205)  becomes: 

dAj/dx  »  i  unj“  exp(-i  lijXA^^I2  +  A^" A^)/^  c  cos  0^  IQ)  (206) 
Similarly,  for  terms  synchronous  with  exp[i(k^  cos  ©^  x  +  ut)]: 
dA^’Vdx  =  i  w  n^"  exp(-i  4j)  (Ap^  +  A^" A^)A^" /(2  c  Iq  cos  ©j)  (207) 
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with  exp[-i(k2  cos  ©2  x  +  ut)]: 

dA^/dx  =  -i  ui  rij"  exp(-i  4^)  (A.^"  +  A^A^A^/tt  c  Iq  cos  ©2)  (208) 
and  with  exp[-i(k2  cos  x  -  tot)]: 

dA2*/dx  =  -i  wn^."  exp(-i  iJijXA.^  +  A2%A^)A^"/(2  c  Iq  cos  ©2)  (209) 
If  in  Figure  10,  ©^  =  ©2  =  ©  so  the  grating  vector  is  in  the  z 
direction  only,  then  in  equations  (206), (207), (208)  and  (209): 

y  =  -i  unj"  exp(-i  Uj)/(2  c  cos  0)  (210) 

To  further  simplify  the  expressions  toward  achieving  a  solution,  the 
following  constants  are  defined  (18:313): 

c  =  A-jA^  + 

di  «  | Aj| 2  +  | A2I  2 

d2  =  |A3|2  +  |A4|2  (211) 

The  coupled  equations  (206),  (207),  (208)  and  (209)  become  (18:314): 
dA^/dx  =  (-7/IQ)[A1d1  -  A1(I1  +  I3)  +  A2  c] 
dA^'/dx  =  (-y/I0)[A1c"  -  A^'Xlj  +  I3)  +  A-^X^] 
dA^/dx  =  (t /Iq) [ A4d2  ~  A4(I4  +  I2)  +  A2,>c] 
dA2*/dx  =  (y/I0)[A4c*  -  A2*(I4  +  I2)  +  A2"d1]  (212) 

Equations  (212)  may  be  rearranged  as  (18:314): 

(d/dxKAj/A^")  =  ( -Aj/Aq"2 )  dA^’/dx  +  (1/A^")  dA^/dx 

=  (-y/Iq)[c  +  (dx  -  d2)  (Aj/Aj*)  -  c*  (Aj/Aj*)2] 
(d/dxXA^/V^’)  =  (-A^/A^'")  c^’/dx  +  (I/A2’ )  dA4/dx 

=  (y/Iq)[c  +  (d2  -  dj)  (A4/A2*)  -  c*  (A4/A2*)2] 

(213) 

By  partial  fraction  expansion,  the  solutions  to  equations  (213)  are 
(18:314): 
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Al/A3*  =  +  +  AK|2)172]  J  exP(~Mx) 

+  [A  +  (A2  +  4|c|2)1/2]  (1/J)  exp(yx) 

2c"  [J  exp(nix)  -  (1/J)  exp(px)]  * 

A4^A2~  =  +  ^2  +  4jo|2)1/2]  K  exp(-px) 

+  [-A  +  (A 2  +  41c!2)1/2]  (1/K)  exp(ux) 

2c"  [K  exp(-yx)  -  (1/K)  exp(yx)]  ^  (214) 

where 

A  d2  _ 

J,K  =  Constants  of  integration 
M  =  y(A2  +  4]  cj  2)1/2/(2  Iq  c*) 

The  phase  conjugate  reflectivity,  R,  can  now  be  solved  using 
equations  (211)  and  (214)  with  boundary  citations:  1^(0) ;  12(0);  13(1); 
and  1^(1)  =  0  where  1  is  the  crystal  length  in  x.  From  these  equations: 

A  =  I3(l)  -  1^0)  -  I2(0)  (215) 

and  Cronin-Golomb  gives  the  solution  for  phase  conjugate  reflectivity  as 
(18:314): 

R  =  I4(0)/I2(0) 

I3(l)  -  I3(l)  A1(l)/A3*(l)  +  I2(0)  a4(0)/a2"(0) 

i2(o)  a1(o)/a3*(o) 

=  4|  c|  2  I  T|  2/|  A  T  +  (A2  +  4|  c|  2  )1/2|2  (216) 

where 

T  =  tarih  U1 
2 

and  j  c |  is  *.ound  by  the  equation: 

[jc|  2  -  It(0)  I3(l)]  |  AT  +  (A2  +  4|  c|  2)1/2j  2  +  4|  c|  2 1  T|  2  I2(0)  I3(l) 

+  2 j  c | 2  I2(0)  (A2  +  4|c|2)1/2  (T  +  T")  =  0 

(217) 

The  functional  relationships  of  this  form  of  the  phase  conjugate 
reflectivity  have  been  buried  in  the  derivation  and  are  not  nearly  as 
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obvious  as  in  Feinberg's  form. 


Summary 

Diffraction  efficiency  and  phase  conjugate  reflectivity  are  the  last 
two  figures  of  merit  for  the  photorefractive  effect  as  applied  to 
optical  phase  conjugation.  Two  theoretical  expressions  were  derived  for 
each  of  these  figures  of  merit. 

The  first  expressions  for  each  (equations  (194)  for  diffraction 
efficiency  and  (197)  for  phase  conjugate  reflectivity)  were  an  extension 
of  Feinberg's  Hopping  Model.  They  showed  that  both  these  figures  of 
merit  were  proportional  to  the  square  of  the  magnitude  of  the 
electro-static  field  set  up  in  the  crystal.  (Thus,  proportional  to  the 
square  of  the  magnitude  of  the  index  grating,  as  well.) 

The  second  expressions  for  each  (equations  (201)  for  diffraction 
efficiency  and  (216)  for  phase  conjugate  reflectivity)  assumed  an  index 
grating  existed  in  the  crystal  and  solved  for  the  light  diffracted  off 
it.  The  diffraction  efficiency  was  shown  in  equation  (201)  to  be 
proportional  to  the  exponential  of  the  magnitude  of  the  diffusion  field, 
only.  The  phase  conjugate  reflectivity  expression  (216)  is  very 
complicated  and  functional  dependencies  are  not  easily  obtained. 

Although  all  three  figures  of  merit,  two  beam  coupling,  diffraction 
efficiency  and  phase  conjugate  reflectivity,  are  considered  important  to 
optical  phase  conjugation,  the  optimum  crystallographic  orientation  of 
for  each  separate  figure  may  be  different  in  a  given  crystal.  A 
specific  example  of  this  is  addressed  in  the  applications  chapter. 
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VII.  Applications  of  the  Theory 


The  previous  two  chapters  developed  theoretical  expressions  for 
steady  state  two  beam  coupling,  diffraction  efficiency  and  phase 
conjugate  reflectivity  since  these  are  the  figures  of  merit  which  will 
quantify  the  photorefractive  effect  as  it  applies  to  optical  phase 
conjugation.  In  this  chapter,  this  theory  is  applied  to  experimentally 
observed  phenoma  with  the  theoretical  predictions  compared  to 
experimental  results  to  determine  the  accuracy  of  the  models. 

Three  main  types  of  applications  will  be  considered.  The  first 
compares  the  functional  relationships  of  the  figures  of  merit,  to 
experimental  data.  The  second  investigates  the  directionality  of  energy 
transfer  in  steady  state  two  beam  coupling.  It  has  been  determined 
experimentally  which  writing  beam  will  gain  intensity  for  a  specific 
crystal  orientation.  The  two  beam  coupling  expressions  developed  here 
will  be  applied  to  predict  this  direction.  Finally,  the  theoretical 
expressions  will  be  applied  to  specific  crystallographic  orientations 
to  determine  which  orientations  are  optimum  for  two  beam  coupling  and 
diffraction  efficiency. 

Relationships  Between  Feinberg's  A I ^/I ^ ,  n  and  R,  and  Eq,  A,  k  and  PWq 

The  expressions  Feinberg  developed  for  intensity  gain  in  steady 
state  two  beam  coupling,  AI^/I^  (equation  (162)),  diffraction 
efficiency,  n  (equation  (194)),  and  phase  conjugate  reflectivity,  R 
(equation  (197)),  are  all  functions  of  several  va/iables.  Ihe 
functional  behavior  of  these  expressions  with  respect  to  some  of  these 
variables  is  investigated  in  this  section. 
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Al^/Ij  as  a  Function  of  the  External  E-field,  Eq 

From  equation  (162),  the  intensity  gain  in  two  beam  coupling  can  be 


written  in  the  form: 

AIj/Ij  -  (1  +  a2  +  f2)/[(l  +  a2)2  +  (af)2]  (217) 

where  f  =  Ep/fQ-  Setting  the  first  derivative  of  equation  (217)  equal 

to  zero,  d(Al^/I^)/df  =  0,  yields  an  inflection  point  at  f  =  0. 

Evaluating  the  second  derivative  of  equation  (217)  at  f  =  0, 

2  2 

d  (Al^/I^)/df  if_Q>  yields  a  positive  result.  Therefore,  the  point 
f  =  0  (Eq  =  0)  is  a  minimum  of  the  intensity  gain  equation. 

In  his  paper,  Gunter,  uses  an  expression  for  two  beam  coupling 
intensity  gain,  Yq  (9:252),  which  can  be  defined  as: 

Y0  =  1  +  A^/Ij  (218) 

Figure  12  plots  experimental  points  of  Yq  versus  Eq  for  BSO  crystals  and 
shows  minima  at  Eq  =  0. 


Figure  12.  Electric  field  dependence  of  two  beam  coupling  intensity  gain 
for  BSO  (9:263) 

n  and  R  as  a  Function  of  the  External  E-field,  Eq 
From  equations  (194)  and  (197),  both  the  diffraction  efficiency  and 
the  phase  conjugate  reflectivity  can  be  written  in  the  form: 

n,  R  *  (a2  +  f2)/[(a2  +  l)2  +  (af)2]  (219) 

The  minimum  of  equation  (219)  with  respect  to  f  is  at  f  =  0.  Therefore, 


DIF  F  RACT  ION  EFTICIENCY  (ARBITRARY) 


Ep  =  0  is  a  minimum  of  both  the  diffraction  efficiency  and  phase 
conjugate  reflectivity.  Figure  13-a  is  an  experimental  plot  by  Feinberg 
of  n  versus  Eq  for  BaTiO^  and  shows  a  minimum  at  Eq  =  0.  Gunter's 
symbol  for  phase  conjugate  reflectivity  is  p  (9:257).  Figure  13-b  has 
experimental  plots  by  Gunter  for  p  versus  Eq  in  BSO  crystals.  They  have 
minima  at  Eq  =  0. 


a.  b. 

Figure  13.  a.)  Diffraction  efficiency  as  a  function  of  applied  voltage 
in  BaTiOj  (5:1304) 

b. )  Phase  conjugate  reflectivity  versus  applied  field  for  BSO 
(9:281) 


i /"> 

as  a  Function  of  External  E-field,  Eq 
Because  Gunter's  plot  for  diffraction  efficiency  versus  external 
field  in  BSO  appears  as  n^2  versus  Eq,  equation  (1941  is  now  writ  Ur  r 
the  form,: 


r1^2  (a2  +  f2)1/2/[  (a2  +  ll2  +  ( .if  )2  1 1 ' 
The  minimum  of  equation  (220)  with  respect  to  f  is  at  f  = 
Eq  =  0  is  a  minimum  of  the  square-root  of  diffract  ion  »  r : ; 

Figure  14  is  an  experimental  plot  b\  -untt  r 
BSO  and  has  minima  at  F()  =  0. 
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Figure  14.  Electric  field  dependence  of  diffraction  efficiency,  n 
for  BSO  (9:263) 


n  and  R  as  a  Function  of  Fringe  Spacing, A 


From  equations  (194)  and  (197),  both  the  diffraction  efficiency  and 
phase  conjugate  reflectivity  can  be  written  as: 

n,  R  a  (a2  +  f2)/[(a2  +  l)2  +  (af)2] 

«  (f2  A4  +  A  A2)/[A4  +  a(2  +  f2)  A2  +  A2]  (221 

where 

«  =  2n/kQ  A 
A  =  (2Vk0)2 

A  minimum  of  equation  (221)  with  respect  to  A  occurs  at  A  =  o.  If 
there  is.no  applied  field  (Eq  =  0,  f  =  0)  in  equation  (221),  an 

inflection  point  occurs  at  A  =  7VkQ  which  is  a  minimum  if 

2  2  1/2 
A  =  4 n  c  kgT/pWQ  q  <  1/5  '  and  a  maximum  if 

4tt2c  c0  kgT/pW0q2  >  1/51/2. 

Figure  15-a  plots  experimental  points  of  n  versus  A  and  Figure  15-b 
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Figure  15.  a.)  Fringe  spacing  dependence  of  diffraction  efficiency  for 
BSO  (9:262) 

b.)  Phase  conjugate  reflectivity  versus  fringe  spacing  for 
BSO.  Solid  lines:  theoretical.  Points:  experimental 
(9:280) 


A Ij/I^  as  a  Function  of  Grating  Vector,  k 

From  equation  (162),  the  intensity  gain  for  two  beam  coupling  may  be 


written  in  the  form: 


A I  f/I  i  “  q2(1  +  a2  +  f 2 )/[  ( 1  +  a2)2  +  (af)2] 


(222) 


where  a  =  k/kQ.  Equation  (222)  has  a  minimum  with  respect  to  a  at  a  =  0 
(k  -  0). 

Figure  16  is  Feinberg's  theoretical  plot  of  equation  (162)  and 
experiment  plot  of  AI^/I^  versus  k  for  BaTiO^.  It  should  have  minima  at 
k  =»  0. 
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Figure  16.  Intensity  gain  for  two  beam  coupling  in  BeTiOo  as  a  function 
of  grating  wave  vector,  k.  Triangles  and  circles  are 
experimental  data,  solid  lines  are  theoretical  fit  to 
equation  (162)  (5:1302) 

n  as  a  Function  of  Grating  Vector,  k 

From  equation  (194),  the  diffraction  efficiency  can  be  written  in  the 
form: 

n  «  (  a  +  f^)/[(a  +  l)^  +  (af)^]  (223) 

where  a  =  k/k.Q.  Equation  (223)  has  an  inflection  point  with  respect  to 
a  at  a  =  0  (k  -  0).  This  point  is  a  maximum  if 
|f  |  =  |Eq  q/(kQkgT)|  >  .64  and  a  minimum  if  |Eq  qAkpkgT)!  <  .64. 

If  there  is  no  applied  field  in  equation  (223)  (Eq  =0,  f  =  0),  it 
has  a  maximum  at  a  =  1  (k  =  [  p  Wq  q^/(  e  Eq  kg 

Figure  17  is  a  theoretical  plot  of  equation  (194)  for  BaTiO^  and 
contains  the  inflection  points  mentioned  here. 

If  the  experimental  geometries  in  Figures  (7)  and  (8)  are 
considered,  the  grating  vector,  k,  may  be  written  as: 

k  =  (4  m  nv/c)sin0  (224) 

if  the  angles  of  incidence  are  equal.  (0^  =  0£  -  0  in  Figure  (7), 

0  ~®-i  =  ©  *n  Figure  (8))-  The  maximum  of  the  diffraction 

efficiency  curve  in  Figure  17  at  k  =  k^  (a  =1),  then  would  imply 
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Figure  17.  Theoretical  plot  of  diffraction  efficiency,  n,  as  a  function 
of  grating  wave  vector,  k  (5:1304) 


maximum  diffraction  efficiency  at: 

sin  ©  =  [pWQ/(e  eQ  kgT)]1,/2  [q  c/(4tr  nv  )]  (225) 

Feinberg  states  in  his  article: 

In  general,  the  charge  density  available  for  migration  will  depend 
on  the  optical  frequency  of  the  writing  beams.  (4:427) 

From  equation  (225),  if  the  angle  of  maximum  efficiency  is  constant  for 

various  optical  frequencies  then  the  density  of  filled  trap  sites,  PWq, 

is  proportional  to  the  square  of  the  optical  frequency,  v.  Feinberg 

continues: 

However,  the  following  rather  strange  fact  is  noted:  for  many 
materials,  studied  in  different  laboratories,  the  diffraction 
efficiency  always  seems  to  peak  for  the  same  full  crossing  angle  of 
about  40  in  air.  According  to  the  model  used  here,  this  implies  a 

i 

16  -3 

similar  nunber  density  of  charges  (10  cm'  )  for  all  these  quite 
different  materials,  a  rather  suspect  coincidence. (4:427) 

The  following  subsection  investigates  diffraction  efficiency  as  c 

function  of  charge  density. 
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n  as  a  Function  of  Charge  Density,  ptJq 

From  equation  (194)  for  diffraction  efficiency,  n  as  a  function  of 
charge  density  can  be  written  as: 

x2 

n  a  - 5 - 5 

(A  +  xr  +  d  (226) 

where 

x  =  Charge  density 
=  pW0 

A  =  k2e  c  0  kBT/q2 
B  =  k  A<f q  ee0/(lq) 

2  2 

Equation  (226)  has  a  minimum  at  x  =  0  and  a  maximum  at  x  =  -(A  +  B  )/A. 
It  is  reasonable  to  have  minimum  diffraction  efficiency  when  there  is  no 
charge  density.  The  maximum  occurs  at: 

pWQ  =  [~ce0/(kB  T)][ (kl  kB  T  +  i  q  A  4>0)/(ql)]2  (227) 

1/2 

R  as  a  Function  of  Spatial  Frequency,  1/A 

Because  Gunter's  plot  for  phase  conjugate  reflectivity  versus 

1/2 

spatial  frequency  in  BSO  is  in  the  form,  p  '  versus  1/A,  equation 
(197)  will  be  written  in  the  form: 

R1^2  *  (a2  +  f2)1//2/[(a2  +  l)2  +  (af)2]1//2  (228) 

where  a  =  2m/(kQA).  Equation  (228)  has  a  minimum  with  respect  to  a  at 
a  =  0  (1/A  =0).  If  there  is  no  applied  field  in  equation  (228) 

(Eq  =0,  f  =  0),  it  has  a  maximum  at  a  =  1 
(1/A  =  [l/(2Tt)][pW0q2/(ee0  A)]1^2). 

Figure  18  has  both  minima  and  maxima  at  1/A  =  0  depending  on  the 

applied  external  field,  Eq. 
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Figure  18.  (Square  root  of)  phase  conjugate  reflectivity  as  a  function 
of  spatial  frequency,  1/A  (9:280) 

Can  the  Direction  of  Energy  Transfer  in  Steady  State  Two-Beam  Coupling 
be  Predicted?  ~  ~  ‘  '  . 

The  direction  of  the  energy  transfer  in  the  steady  state  two-beam 
coupling  described  in  Chapter  V  has  been  observed  to  be  dependent  on  the 
orientation  of  the  crystal's  c-axis  relative  to  the  incident  beams 
(4:431).  Figure  19  shows  the  direction  of  energy  transfer  in  BaTiO^.  The 
theory  developed  in  Chapter  V  will  be  applied  to  two  crystals,  BaTiO-^ 
and  BSO,  to  determine  if  coupling  can  be  predicted. 

Feinberg's  Model 

The  fractional  intensity  gain  expression  developed  by  Feinberg  in 
Chapter  V  (equation  (162))  took  the  form: 

&I1(x)/I1(x)  a  (e1'',ew*R*k’e^.*e2)  (229) 

The  contracted  electro-optic  coefficient  tensor,  R,  is  referenced  to  the 
primary  coordinate  system  (PCS)  of  the  crystal  so  k  must  be  referenced 
to  the  same  coordinate  system  in  order  to  perform  the  dot  product. 

Initially,  3ince  k  is  defined  from  the  interference  expression 
(equation  (98)),  it  appears  k  can  be  defined  as  either  +k  without  any 
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Figure  19.  Directional  two-beam  coupling  as  observed  in  BaTiO^  (a 
photoref ractive  crystal  with  positive  charge  carriers.) 
(3:432) 


consequences  (k  =  -  k-2  or  k  =  -  k^ ).  However,  from  the  defined 

space-charge  field  expression  (equation  (115)),  a  unique  direction  is 
given  to  k  and  that  is  in  the  same  direction  as  a  positive  phase  shift, 
il>.  The  space-charge  field  nay  be  written  in  the  form: 

Esc(x)  =  fc  [m  pWQ  q/(ce0  k)]  |l/[l  +  PWQ  q2/(ee0  kgT  k2)]| 

cos  (k'x  +  ti/2)  (230) 

if  there  is  no  applied  field  (Eq  =  0)  so  the  phase  shift  is  ti/2 
(y  =  tt/2).  In  this  form,  it  can  be  seen  that  the  direction  of  the 
grating  shift  is  determined  by  the  sign  of  the  carrier,  q. 
Mathematically,  the  phase  shift  is  tt/2  for  both  positive  and  negative 
carriers.  However,  a  negative  carrier  will  change  the  sign  of  the 
magnitude  of  Egc  (or  give  it  a  tt  phase  charge)  and  give  it  the 
appearance  of  a  -m/2  phase  shift.  Figure  20  exemplifies  this. 

From  Figure  20,  it  is  seen  that  k  is  always  in  the  direction 
opposite  that  of  the  positive  E-field.  Since  R  and  k  are  the  only  terms 
in  the  gain  expression  (equations  (163)  and  (164))  which  contain 
directionality,  when  the  space-charge  field  is  reference  to  the  same 
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coordinate  system  as  R,  R*k  is  always  a  negative  expression  and  AI^/I^ 
(equation  (162))  positive. 


a)  Intensity  b)  Charge  density 


-negative  carrier 
(BSO) 


c)  E-field 
(positive 
referenced 
up) 


d)  E-field 
with  ip  =  0 


Figure  20.  Illustration  of  the  photorefractive  effect  for  positive  and 
negative  carriers  under  diffusion  only  (Eq  =  0)  showing 
direction  of  phase  shift. 


Feinberg's  model  for  steady  state  two-beam  coupling  may  define  the 
magnitude  of  the  coupling,  but  it  does  not  describe  the  directionality. 

Application  of  Two  Beam  Coupling  Assuming  a  Sinusoidal  Index  Grating 

Exists 

The  expressions  for  the  intensities,  1^  and  I2  developed  by  White  in 
Chapter  V,  equations  (175),  can  be  written  in  the  forms: 

Ij(r)  a  1/[I1(0)  +  I2(0)  exp(-2rr)l 
I2(r)  a  1/[I1(0)  exp(2 Ir)  +  I2(0)]  (231) 

From  these  forms  it  can  be  seen  that  Ij  is  directly  related  to  T  and  I2 
is  inversely  related  to  r. 

r  is  proportional  to  n^  (after  equation  (172)),  the  modulated  index 
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(after  equation  (165)).  For  a  space-charge  field  in  z  only,  is 

inversely  proportional  to  that  field,  E^,  for  all  polarizations. 

Therefore,  1^  is  inversely  related  to  and  ^  is  directly  related  to 

E  . 
z 

In  Figure  19,  positive  E is  taken  to  be  in  the  direction  of  the 
c-axis  for  both  configurations  (a)  and  (b).  However,  the  orientation  of 
Ez  to  the  propagation  vectors  and  is  the  opposite  in  (a)  and  (b). 
Since  configuration  (a)  is  used  in  the  development  of  Chapter  V,  a 
direct  substitution  for  BaTiO^  from  Appendix  A  into  equations  (231) 
shows  that  Ij  should  be  decreasing  as  it  propagates  and  ^  increasing. 

In  Figure  19-b,  the  effects  of  flipping  the  crystal  is  that  Ez  now 
looks  like  -E  in  the  Appendix  A  expressions.  Substitution  into 
equations  (231)  shows  that  1^  should  increase  as  it  propagates  while  ^ 
decreases. 

The  same  argument  can  be  applied  to  the  BSO  expressions  in  Appendix 
A  though  not  as  simply  since  the  PCS  is  rotated  and  the  index  ellipsoid 
is  expanded  along  one  new  axis  but  contracted  along  the  other. 
Polarizations  must  be  tracked  to  determine  the  index  a  beam  encounters 
but  the  direction  of  coupling  is  still  attainable. 

Summary 

The  expression  for  steady  state  two  beam  coupling  developed  from 
Feinberg's  hopping  model  gives  the  magnitude  of  the  gain  or  loss  of 
intensity  by  the  beams  but  it  does  not  indicate  which  beam  is  gaining 
or  losing.  The  direction  of  the  coupling  can,  however,  be  predicted  from 
a  development  in  which  an  index  grating  is  assumed  to  exist  and  the 
solution  for  propagation  through  this  grating  is  found  through  the  wave 
equation.  By  using  this  solution  and  the  proper  index  ellipsoid 
relationship,  the  direction  of  coupling  as  well  as  its  magnitude  can  be 


found. 


Why  is  There  Better  Two  Beam  Coupling,  and  Diffraction  Efficiency 
Crystallographic  Orientations  in  BSO? 

It  has  been  observed  that  there  are  different  crystallographic 
orientations  of  Bismuth  Silicate  (BSO)  which  are  optimum  for  energy 
transfer  in  two  beam  coupling  and  diffraction  efficiency  for  vertically 
polarized  incident  light.  Figure  21  presents  these  optimum 
orientations.  For  maximum  beam  coupling,  the  space  charge  field  is 
formed  perpendicular  to  the  001  face,  while  for  maximum  diffraction 
efficiency,  it  is  perpendicular  to  the  110  face  (21:131-133). 


Figure  21.  Optimum  crystallographic  orientations  in  BSO  (19:132). 


The  theory  developed  in  Chapters  V  and  VI  will  be  applied  to  BSO  to 
determine  if  these  optimum  orientations  can  be  predicted. 

Feinberg’s  Model 

From  crystallography,  the  primary  coordinate  system  (PCS)  of  BSO 
with  no  applied  electric  field  is  shown  in  Figure  22  (19:30-40). 

From. equations  (162)  and  (164),  the  intensity  gain  for  one  of  the 
beams  in  steady  state  two  beam  coupling  in  BSO  can  be  written  in  the 
form: 


V1! 


[(k  e~  +  k  e0  )^  +  (k  e0  +  k  e„  )^ 

'  z  2y  y  2z  z  2x  x  2z 

+  (ky  e2x  +  kx  e2y)2]  (elx  e2x  +  ely  e2y  +  e 


lz  +  e2z} 


(232) 
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Figure  22.  PCS  in  BSO  with  no  applied  field. 


From  equations  (194)  and  (199),  the  diffraction  efficiency  in  BSO  can  be 
written  in  the  form: 


2  2  2 
0  «  (k  e,  +  k  eQ  )  +  (k  e-.  +  k  e,  )  +  (k  e-j  +  k  e,  ) 

z  3y  y  3z  '  z  3x  x  3z'  y  3x  x  3y' 


(233) 


In  the  optimum  two  beam  coupling  orientation  (Figures  21-a  and  22-a), 

the  grating  direction  vector,  k,  is  in  the  z  direction  only.  Vertically 

polarized  light  with  unit  polarization  vectors,  e ^ ,  e?  and  e^  has  only  x 

1/2 

and  y  components  which  are  1/2  each.  Equations  (232)  and  (233),  for 
this  orientation  become: 


AI1/I1  *  k 

n  «  k2  (234) 

In  the  optimum  diffraction  efficiency  orientation  (Figures  21-b  and 

22-b),  the  grating  vector,  k,  has  k^>  and  k  y  components  where  kx  and  k 
1/2 

are  k/2  ’  Vertically  polarized  light  has  ezz  components  only,  so 

equations  (232)  and  (233)  become  the  same  as  equations  (234). 

The  relationships  for  two  beam  coupling  and  diffraction  efficiency 
from  Feinberg's  hopping  model,  then,  do  not  adequately  predict  optimum 
crystal  orientations  for  vertical  input  polarizations. 


An  Index  Ellipsoid  Approach 


The  PCS  of  BSO  without  an  applied  field  was  shown  in  Figure  22.  The 
direction  of  the  k  vectors  and  thus  the  space  charge  fields  are  shown  in 
Figure  21.  For  Figures  21-a  and  22-a,  the  E-field  is  applied  in  the  z 
direction  only.  From  Appendix  A,  the  new  indices  of  refraction  are: 

n'  =  n±  (1/2)  n3  Esc  (235) 

with  new  PCS  resulting  from  a  rotation  about  the  z-axis  of  -135°  (Figure 
23-a)  (20:132). 


In  Figures  21-b  and  22-b,  the  space-charge  field  has  x  and  y  components 
so  the  index  ellipsoid  is: 

(x2  +  y2  +  z2)/n2  +  2r41  E^yz  +  2r41  Eyxz  =1  (236) 
in  the  PCS  of  the  crystal  without  an  applied  field.  The  eigenvalues  can 
be  found  by  solving  the  determinant: 
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(1/nj  -  X 
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(237) 


where 


X  =  Eigenvalue 
=  1/n’2 

n'  =  New  indices  of  refraction 

With  this  applied  field,  the  crystal  becomes  biaxial  with  indices  of 

3 

refraction  on  the  axes  of  the  new  PCS  being  n  and  n±  (1/2)  n  r4^  E^. 
The  axes  rotation  between  the  old  and  new  PCS  is  found  by  solving: 
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where 
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0. . 

ij 

The  axes 


=  Eigenvalues  solved  for  in  equation  (237) 

=  Elements  of  rotation  matrix 

=  cos  0. . 
ij 

=  Angle  of  rotation  between  old  i^‘  axis  and  new 
rotation  is  shown  in  Figure  23-b  (19:132) 


axis 


Figure  23. 


New  primary  coordinate  system  and  indices  of  reiracton  in 
BSO  with  applied  space-charge  fields  (a)  Eg^  =|  Eg J  "2 

(b)  EgC  =  (|  Esc|  /21/2)x  +  (|Esc|/21/2)y  (19:132) 


For  diffraction  efficiency  which  is  dependent  on  the  magnitude  of 
the  index  grating,  in  the  Figure  23-a  orientation, .vertically  polarized 
light  Is.  polarized  along  x'  only  and  sees  index,  n  -  (1/2)  n  r^j  Eg^. 

In  the  Figure  23-b  orientation,  vertically  polarized  light  has  both  §x' 

A  3 

and  e  '  components  which  see  n  +  (1/2)  n  r^  Eg^  and 

n  -  (1/2)  n  r^  Eg^,  respectively.  At  points  where  the  maxima  or 

minima  of  each  component  coincide,  the  summation  of  the  components  gives 
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a  modulated  index  of  ±  n2  r^  (19:132). 

Therefore,  the  Figure  23-b  orientation  can  be  predicted  to  be  the 
better  of  the  two  for  diffraction  efficiency. 

For  two  beam  coupling,  the  Figure  23-a  orientation  is  precisely  the 
configuration  assumed  for  the  Chapter  V  derivation  so  significant  energy 
transfer  can  be  assumed.  In  the  Figure  23-b  configuration,  however, 
since  the  polarization  components  are  split  between  two  indices,  one  of 
which  has  a  +n^  term  (as  described  in  equation  (165))  and  the  other  -n^, 
one  component  would  tend  to  gain  energy  (as  in  equations  (175))  while 
the  other  would  lose  energy.  Each  component,  then,  tends  to  compensate 
for  the  other  and  energy  transfer  is  not  efficient  (19:132). 

Therefore,  the  Figure  23-b  orientation  can  be  predicted  to  be  the 
better  for  energy  transfer  in  two  beam  coupling  of  the  two. 

AFIT  Crystal 

Another  pair  of  crystal  orientations  which  has  one  orientation 
better  for  two  beam  coupling  and  another  for  diffraction  efficiency  can 
be  analyzed  by  considering  beams  incident  on  the  110  face.  The  incident 
light  is  vertically  polarized  and  forms  a  k  vector  oriented  as  shown  in 
Figure  24. 

In  the  Figure  24-a  orientation,  the  space  charge  field  is  in  the  z 
direction  so  the  analysis  is  identical  to  the  Figures  21-a,  22-a,  and 
23-a  developments  in  the  preceding  subsections.  The  new  PCS  and 
appropriate  indices  of  refraction  are  shown  in  Figure  25-a. 

For  the  Figure  25-b  orientation,  the  space  charge  field  has  both  x 
and  y  components,  however,  these  components  are  of  opposite  sign.  The 
index  ellipsoid  in  the  PCS  of  the  crystal  without  an  applied  field  is  in 
the  form: 

(x2  +  y2  +  z2)/n2  +  2r^  Eg^/2^2  xz  -  2r^  Eg^/2^2  yz  =  1  (239) 
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Figure  24.  Potential  incident  light  and  grating  vector  orientations  when 
only  the  110  face  of  BSO  is  available  for  illumination. 


(a) 


=  n  +  (1/2)  n^ 
r41  ESC 


Figure  25.  New  PCS  and  appropriate  indices  of  refraction  for  a  BSO 

crystal  illuminated  on  the  110  face  such  that  the  applied 

2  and 


space-charge  field  is  (a)  E^  =  E, 

(b)  E^c  =  (-|Esc|/21/2)x  +  (|Esc|/2i/^)y 


'SC 

1/2  NO 


The  new  indices  of  refraction  and  angles  of  rotation  can  be  found  by  a 
derivation  similar  to  that  in  the  previous  subsection.  Tbe  new  indices 
and  PCS  are  shown  in  Figure  25-b. 

By  similar  analysis  as  in  the  previous  subsection,  optically 
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polarized  light  incident  on  Figure  24-a  orientation  sees  index 
2 

n  +  (1/2)  n  r^  and  should  provide  significant  energy  transfer. 

The  same  light  incident  on  Figure  24-b  orientation  has  components  which 
3 

see  n  ±  (l/2)n  r^  Eg^.  At  points  where  the  maxima  and  minima  of 
these  components  coincide,  therefore,  the  total  modulated  index  is 

3 

±  n  r^  Eg£  and  this  configuration  should  produce  better  diffraction 
efficiency.  However,  because  of  the  change  in  sign  of  the  modulated 
indices,  the  energy  transfer  in  two  beam  coupling  should  be  minimal. 
Summary 

Just  as  in  the  preceding  section,  Feinberg's  hopping  model  appears 
to  be  fine  for  predicting  the  space-charge  field  and  phase  shift  between 
the  index  grating  and  interference  fringes  but  the  expressions  he 
develops  from  the  model  for  diffraction  efficiency  and  steady  state  two 
beam  coupling  do  not  adequately  predict  these  effects  for  various 
crystal  orientations.  The  directionality  of  these  effects,  however,  can 
be  accurately  predicted  by  using  an  index  ellipsoid  approach. 

Summary 

Three  applications  of  the  three  figures  of  merit,  two  beam  coupling, 
diffraction  efficiency  and  phase  conjugate  reflectivity,  were  performed 
in  this  chapter.  The  first  application  compared  the  functional 
relationships  of  the  figures  of  merit  developed  by  Feinberg  to 
experimental  data.  These  comparisons  showed  Feinberg's  expressions  to 

adequately  predict  the  functional  relationships  of  empirical  data. 

/ 

The  second  application  was  to  predict  the  direction  of  energy 
transfer  in  two  beam  coupling.  Feinberg' s  expressions  proved  inadequate 
in  this  respect.  However,  the  expressions  obtained  in  Chapter  V  by 
assuming  an  index  grating  existed  in  the  medium  did  accurately  predict 
the  experimentally  observed  directionality. 
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The  final  application  was  to  determine  why  one  crystallographic 
orientation  of  BSO  was  better  for  two  beam  coupling  while  another  was 
better  for  diffraction  efficiency.  A^ain,  Feinberg's  expressions  could 
not  predict  this  difference.  This  time,  however,  an  index  ellipsoid 
approach  was  used  and  the  experimental  observations  were  accurately 
predicted. 


VIII.  Conclusions  and  Recommendations 


Conclusions 

Optical  phase  conjugation  is  a  nonlinear  optical  process  which,  from 
a  set  of  incident  light  beams,  produces  a  beam  which  acts  as  the  time 
reversal  of  one  of  the  incident  beams.  This  process  has  application 
potential  because  phase  aberrations  added  to  the  incident  wavefront  as 
it  propagates  through  some  medium  would  then  be  removed  from  the  time 
reversed  wavefront  propagating  through  that  same  medium  (provided  the 
medium  remained  unchanged). 

The  first  type  of  optical  phase  conjugation  was  by  degenerate  four 
wave  mixing.  Three  electro-magnetic  waves  all  of  the  same  frequency, 
are  incident  on  a  nonlinear  material  oriented  properly.  They  cause  the 
polarization  of  the  medium  in  both  the  linear  and  the  nonlinear  orders. 
The  third  order  nonlinear  polarization  generates  a  new  E-M  wave  whose 
amplitude  is  given  through  Maxwell’s  wave  equation.  By  the  conservation 
of  energy  and  momentum,  this  wave  must  be  of  the  same  frequency  as  the 
incident  light  and  counterpropagating  to  the  probe  if  the  two  pump  beams 
are  counterpropagating  in  order  for  it  to  attain  significant  energy. 

Since  it  is  the  electric  fields  of  the  incident  light  which  drive 
the  medium  into  nonlinearity  and  large  E-fields  are  required  to  do  this, 
high  intensity  incident  light  is  required  to  produce  optical  phase 
conjugation  by  degenerate  four  wave  mixing.  A  technique  for  producing 
phase  conjugation  with  lower  intensity  incident  light  would  then  be 
desirable.  This  technique  is  through  the  use  of  photorefractive 
crystals. 

Photorefractive  crystals  have  various  sites  in  their  lattice.  When 
under  illumination,  charges  trapped  at  these  sites  are  free  to  move  from 
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site  to  site  and  finally  be  rebound  at  a  location  which  is  not 
illuminated.  When  two  incident  coherent  light  beams  are  interfered  in  a 
pho tore fr active  crystal,  the  fringe  pattern  causes  a  depletion  of 
charges  at  the  interference  maxima  and  a  high  concentration  of  charges 
at  the  minima.  This  separation  of  charge  produces  an  electro-static 
field  in  the  crystal  which  modulates  the  index  of  refracion  of  the 
crystal  by  the  linear  electro-optic  (Pockel's)  effect. 

When  a  photorefractive  crystal  is  oriented  properly,  a  third  light 
beam  which  is  counter  propagating  to  one  of  the  incident  (interfering) 
beams  diffracts  off  the  index  grating.  The  diffracted  beam 
counterpropagates  to  the  other  incident  (interfering)  beam  and  is  the 
phase  conjugate  of  that  beam.  Even  if  low  intensity  incident  light  is 
considered,  their  interference  still  produces  a  fringe  pattern  which 
induces  an  electro-static  field.  Since  it  is  this  induced  field  (rather 
than  the  E-fields  of  the  incident  light)  which  drives  the  crystal 
nonlinearity,  optical  phase  conjugation  may  occur  in  photorefractive 
crystals  with  very  low  intensity  incident  light. 

Two  models  of  the  photorefractive  effect  are  Feinberg's  Hopping 
Model  and  Kukhtarev's  Solid  State  Model.  Feinberg's  Model  is  based  on 
the  idea  that  a  charge  "hops"  from  site  to  site  in  the  crystal  and 
develops  the  probability  of  this  charge  movement.  Kukhtarev's  model  is 
a  more  standard  development  of  charge  transport.  The  significance  of 
Kukhtarev's  model  is  that  he  derives  the  conditions  for  which  spatially 
sinusoidal  incident  light  (interference  fringes)  produces  an 
electro-static  field  which  is  also  sinusoidal  with  the  same  spatial 
frequency  as  the  fringe  pattern. 

Both  Feinberg  and  Kukhtarev  concluded  that,  for  incident  interfering 
waves  with  a  small  modulation  index,  the  electro-static  field  would  be 


sinusoidal  of  the  same  spatial  frequency  as  the  interference  pattern  but 
spatially  shifted  in  phase  from  it. 

The  spatial  phase  shift  between  the  interference  fringes  and  the 
electro-static  field  is  important  for  two  beam  coupling,  one  of  the 
figures  of  merit  for  the  photorefractive  effect  as  applied  to  optical 
phase  conjugation.  Two  beam  coupling  is  the  coupling  of  energy  from  one 
of  the  incident  (interfering)  beams  which  writes  the  index  grating  to 
the  other  as  they  propagate  through  this  grating  they  have  written. 

Two  theoretical  expressions  for  two  beam  coupling  were  developied. 

The  first  is  an  extension  of  Feinberg's  Hopping  Model.  The  second 
assumes  an  index  grating  exists  and  solves  the  wave  equation  for  light 
propagating  through  its  grating.  Although  each  expression  is  uniquely 
developed,  they  both  indicate  that  if  no  spatial  phase  shift  exists 
between  the  interference  fringes  and  the  electro-static  field,  no  beam 
coupling  will  occur. 

Two  other  figures  of  merit  for  the  photorefractive  effect  as  applied 
to  optical  phase  conjugation  are  diffraction  efficiency  and  phase 
conjugate  reflectivity.  Two  theoretical  expressions  for  each  of  these 
figures  of  merit  were  also  developed.  Diffraction  efficiency  is  the 
amount  of  power  diffracted  off  a  grating  at  the  Bragg  angle  compared  to 
the  amount  of  power  in  the  wave  incident  on  the  grating.  Phase 
conjugate  reflectivity  is  the  amount  of  power  in  the  conjugate  wave 
compared  to  the  amount  of  power  in  the  wave  which  has  been  conjugated. 

The  first  expressions  for  each  of  these  figures  of  merit  were 
developed  from  Feinberg's  Hopping  Model.  They  showed  that  these  figures 
were  proportional  to  the  square  of  the  magnitude  of  the  electro-static 
field  set  up.  The  second  expressions  for  each  assumed  an  index  grating 
existed  and  asked  for  light  diffracted  off  this  grating.  No  clear  cut 
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functional  dependencies  on  the  electro-static  field  could  be  made  on  the 
expressions  derived  in  this  way. 

When  the  figures  of  merit  developed  by  Feinberg  were  compared  to 
empirical  data,  they  were  found  to  accurately  predict  the  functional 
dependencies  of  this  data  with  respect  to  external  electric  field 
applied  to  the  crystal,  interference  fringe  spacing  and  charge  density. 
However,  when  Feinberg 's  expression  for  two  beam  coupling  was  used  to 
predict  the  direction  of  the  coupling  for  a  particular  crystallographic 
configuration,  it  proved  to  be  inadequate.  Also,  when  Feinberg' s 
expressions  for  two  beam  coupling  and  diffraction  efficiency  were  used 
to  determine  the  magnitude  of  each  with  respect  to  various 
crystallographic  orientations  of  Bismuth  Silicate  (BSO),  they  proved  to 
be  inaccurate  when  compared  to  empirical  data. 

The  conclusion,  then,  is  that  the  functional  dependencies  of 
Feinberg' s  figures  of  merit  appear  to  be  correct  when  applied  to 
empirical  data.  However,  if  the  absolute  magnitude  of  these  figures  or 
even  relative  magnitudes  between  the  different  orientations  are  desired, 
there  are  crystallographic  orientations  for  which  Feinberg's  figures  of 
merit  are  inaccurate. 

By  using  the  expression  for  two  beam  coupling  developed  assuming  a 
sinusoidal  index  grating  existed  in  the  medium,  correct  predictions  for 
the  directionality  of  the  coupling  can  be  made.  Predictions  for  the 
relative  magnitudes  of  two  beam  coupling  between  different 
crystallographic  orientations  can  also  be  accurately  made  using  this 
expression.  The  intensity  gain  (or  loss)  of  one  beam  is  a  function  of 
the  scaled  exponential  of  the  magnitude  of  the  index  grating  while  the 
loss  (or  gain)  of  the  other  is  a  function  of  the  negative  of  this 
exponential.  Solutions  are  made  by  determining  the  modulated  index  with 


an  index-ellipsoid  approach. 

The  index-ellipsoid  approach  is  also  useful  in  determining  the 
relative  diffraction  efficiencies  between  different  crystallographic 
orientations.  Feinberg's  expression  showed  the  diffraction  efficiency 
to  be  proportional  to  the  square  of  the  magnitude  of  the  electro-static 
field  (and  thus  to  the  square  of  the  magnitude  of  the  modulated 
refractive  index). 

Recommendations  for  Further  Research 

This  thesis  is  not  a  complete  compilation  of  all  aspects  of  the 
photorefractive  effect  applied  to  optical  phase  conjugation.  Some 
questionable  areas  are  listed  here. 

Experimental  Data  Versus  Theory 

Most  published  papers  in  optical  phase  conjugation  and  other 
photorefractive  phenomena  include  experimental  data  points  to  which  the 
theory  is  curve  fit.  Evaluations  of  the  theory  derived  in  this  thesis 
were  made  compared  to  this  data.  Since  AFIT  can  now  take  photorefractive 
measurements  in  house,  an  attempt  should  be  made  at  applying  the  theory 
to  in-house  data.  Such  research  should  be  very  insightful  in  gaining  an 
understanding  of  the  theory. 

Self-Pumped  Phase  Conjugate  Mirrors 

Since  the  fundamentals  of  photorefraction  and  phase  conjugation  have 
been  laid  out  in  this  thesis,  the  next  logical  step  seems  to  be  to 
extend  the  tutorial  to  self-pumped  phase  conjugate  mirrors  (SPPCMs). 
SPPCMs  greatly  broaden  the  spectrum  of  PCM  applications  because  they  do 
away  with  the  need  for  external  coherent  pump  sources.  Several  papers 
have  been  written  in  the  area  which  are  open  to  interpretation.  For  a 
listing  of  these  papers  and  an  overview  of  self-pumped  phase 
conjugation,  refer  to  Appendix  B. 
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Transient  Solutions  to  the  Photorefractive  Effect 

All  the  analyses  reported  in  this  thesis  are  steady  state  solutions. 
Further  research  is  needed  looking  into  the  transient  solutions.  Because 
of  the  nature  of  the  photorefractive  effect,  it  is  reasonable  to  suspect 
that  incident  intensity  and  time  might  be  inversely  related  in  both 
creating  and  erasing  the  diffraction  grating.  Not  only  does  the  theory 
behind  this  suspicion  require  further  analysis,  but  it  would  be  an 
excellent  area  in  which  to  make  empirical  measurements  to  compare  to  the 
theory . 

Appendix  B  presents  Feinberg's  analysis  of  photorefractive  time 
decay  and  some  empirical  data  from  other  sources. 

How  Much  Does  Optical  Activity  Effect  BSO  Analyses? 

Optical  activity  is  the  rotation  of  the  plane  of  polarization  of 
linearly  polarized  light  as  it  passes  through  a  medium.  The  polarization 
state  of  light  is  considered  in  nearly  all  the  analyses  reported  in  this 
thesis.  Further  research  is  needed  to  include  optical  activity  in  the 
tracking  of  polarization  as  light  passes  through  the  crystal  and  to 
determine  the  effect  it  has  on  the  analyses  and  derivations,  themselves. 
Appendix  B  contains  an  overview  of  optical  activity. 

Can  Boundary  Conditions  be  Ignored  in  Photorefractive  Analyses? 

In  the  theory  presented  in  this  thesis  and  in  most  published  papers 
it  appears  boundary  conditions  (i.e.  Snell's  Law-type  diffraction  at 
surfaces)  are  omitted  and  assumed  to  be  inconsequential.  A  superficial 
analysis  of  this  reveals  that  in  a  Figure  9- type  geometry,  the  grating 
wave  vector  is  given  by: 
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k  kj  -  k^ 

=  (2-rrn^A)  cos  a^'  cos  0^'  -  2nn2A 

=  [(2TTn^/X)cos  a^'  cos  0^'  -  ^t^A^os  a ^  cos  ©2']  * 

+  [(2ttii^A)  cos  a^'  sin  0^'  -  (21^  A)  cos  02'  sin  ©2 '  ]  y 

-  [(2-rrn^A)  sin  a^'  +  {l-nn^/v)  sin  ©2 * J  z 

=  [ ( 2-rrn ^/X )  cos  a^'  cos  0^'  -  (2nn2A)  cos  02'  cos  02*]  x 

+  [(2ttA)  cos  a^'  sin  0^  -  (2nA)  cos  02'  sin  y 

-  (2tt/A )  (sin  +  sin  02)  z  (240) 

where 

Op  a2,  Op  ©2  =  Angles  of  incidence  in  air 
a^',  02',  ©i',  ©2'  =  Angles  of  interference  in  crystal 
Qj'  =  sin_1[(sin  a1)/n1] 

02'  =  sin  *[(sin  02)702] 

©j'  =  sin  1  [(sin  ©1)/n1) 

©2'  =  sin  [(sin  ©2 ) /n2 1 

From  this  rough  look  at  the  problem,  there  is  a  difference  between  using 
the  interface  angles  in  the  crystal  and  the  angles  of  incidence  outside 
the  crystal.  Further  research  is  required  to  determine  the  magnitude  of 
this  difference  or  any  other  problems  this  difference  might  lead  to. 
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Appendix  A:  The  Index  Ellipsoid 

In  any  crystal,  the  index  of  refraction  along  an  axis  of  the 
principal  coordinate  system  (PCS)  may  be  used  to  create  an  "index 
elliposoid"  as  follows: 

x2/n  2  +  y2/n  2  +  z2/n  2  =  1  (241) 

x  ’  y  z 

Figure  26  represents  equation  (241)  pictorially.  In  an  isotropic  crystal, 
the  three  indices  are  the  same  so  the  ellipsoid  is  a  sphere.  In  a 
uniaxial  crystal  two  indices  are  the  same  and  one  is  differents  so  the 
axis  of  the  different  index  is  an  optic  axis.  In  a  biaxial  crystal,  all 
three  indices  are  different  and  the  optic  axes  are  off  the  principal 
axes. 


x 


Figure  26.  Illustration  of  the  index  ellipsoid. 

The  polarization  of  light  passing  through  the  crystal  is  broken  up 
into  components  along  the  PCS  to  determine  which  index  each  component 
sees.  For  a  propagation  direction,  s,  not  along  a  principal  axis,  the 
index  ellipsoid  must  be  sliced  perpendicular  to  s  to  determine  the 
indices  each  polarization  component  sees  (Figure  27). 
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X 


Figure  27.  Appropriate  indices  on  index  ellipsoid  for  a  propagation 
direction,  s. 


When  an  external  dc  electric  field,  E,  is  applied  to  the  crystal, 
the  index  ellipsoid  is  deformed  so  its  representation  in  the  original 

PC'  i s : 

x"(,l/n  L  +  Ab, )  +  y^(l/n  ^  +  Ak  )  +  z^(l/n  ^  +  Ab,) 
x  1  y  2  z  3 

+  2  Ab^  yz  +  2  Ab-  xz  +  2  Ab^  xy  =  1  (242) 

where 

th 

Ab^  =  icn  component  of  Ab  tensor  (in  contracted  form) 

Ab  =  R-E 

R  =  Contracted  electro-optic  coefficient  tensor 
The  new  ellipsoid,  of  course,  also  can  be  written  in  a  form  similar  to 
equation  (241)  in  a  new  coordinate  system  which  is  rotated  from  the 
original  PC.5. 

To  find  the  indices  on  axes  of  the  new  coordinate  system  the 
following  determinant  is  solved  for  the  eigenvalues,  : 

1/n  1  +  Ab,  -  a  Ab.  Abr 

x  1  A  s 

Ab6 
Ab^ 

where 

H  *  VV2 


1/n  +  Ab.,  -  X 

>'  2 

Ab. 


Ab, 


1/n  *  +  Ab3  -  X 


(243) 
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n .  '  =  Index  in  new  PCS 
1 


The  rotation  of  the  coordinate  system  is  found  by  solving: 


l/nx  +  Abx 
Ab, 


Abr 


X  . 
1 


4b6 

l/ny2  *  Sb2  -  X, 
Ab, 


Abc 


Ab, 


1/n  +  Ab,  -  X . 

z  3  1 


“1 

“  — 

ail 

ai2 

_ai3 

«  0 


(244) 


for  a. .  where 
ij 

a^j  =  Element  of  rotation  matrix 


=  Cos  0. 


ij 


.  th 


.  th 


0^  =  Angle  between  new  i  axis  and  old  j  axis. 

The  following  are  specific  examples  of  the  new  indices  of  refraction 
given  a  space-charge  field  is  set  up  along  one  of  the  original  principal 
axes: 


BaTiOy 


■51 

0 


0 

0 

0 

r51 

0 

0 


13 

r13 

r33 

0 

0 

0  / 


If  ' 


4bi\ 

ib2 

Ab, 

AbA 

Ab5 

Ab, 


(245) 


If  the  grating  is  set  up  along  x,  the  index  ellipsoid  in  the  original 
PCS  is: 

n  o  o  o  o 

(246) 


(x2  +  y2)/nQ2  +  z2/ng2  +  2r^Evxz  =  1 


-  ca  r 

51  x 

The  rotation  is  in  the  x-z  plane  so  n  '  =  n  and  the  eigenvalues  are: 
X  =  {l/n„2  +  l/no2  ±  [(1/n^2  +  l/nD2)2  +  4r512Ex2/(n^2)]1/2}  /2 


o  e 


If  the  grating  is  along  y,  the  index  ellipsoid  is: 

+  y2)/no2  +  z2/ne2  +  2r51  Ey  yz  =  1 


(247) 

(248) 
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Rotation  is  in  the  y-z  plane  so  n^'  =  nQ  and  the  eigenvalues  are  the 
same  as  equation  (247)  replacing  with  E  .  If  the  grating  is  along 

z,  the  index  ellipsoid  is  : 

x2(l/no2  +  r13Ez)  +  y2(l/nQ2  +  r13Ez)  +  z2(l/ne2  +  r33Ez)  =  1  (249) 
The  principal  axes  remain  the  same  with  the  new  indices  of  refraction 


being: 


n  '  =  n  '  =  n  /[l  +  n  2  r.,  ,E  j^2 
x  y  o  L  o  13  zJ 

s  n  -  (1/2)  n  2  r1  E 
o  o  13  z 

n  '  =  n  /[  1  +  n  2  r.,.,  E  j^"22 
z  e  L  e  33  zJ 

s  n  -  (1/2)  n  2  r^  E 
e  e  33  z 


(250) 


(Note  that  this  is  the  form  of  the  index  assumed  in  Chapters  IV,  V  and 
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If  the  grating  is  along  x,  the  index  ellipsoid  is: 

(x2  +  y2  +  z2)/n2  +  2r41  Ex  yz  =  1 
Rotation  is  in  the  y-z  plane  so  n  1  =  n,  the  eigenvalues  are: 


X  =  lAi  t  r, ,  E 
41  x 


(252) 


(253) 


and  the  new  indices  are: 


n  ’  =  n/[l±  n2  r,  1  E  J^2 
y,z  L  41  xJ 

=  n  *  (1/2)  n3  r51  Ex 

(Again,  this  is  the  form  assumed  in  Chapters  IV,  V  and  VI.) 


(254) 


If  the  grating  is  along  y  or  z,  the  index  ellipsoids  are: 


(255) 


(x2  +  y2  +  z2)/n 2  +  2  r^  E^  xz  =  1 

(x  +  y  +  z  )/n  +2  E£  xy  =  1 

respectively.  Rotation  is  in  the  x-z  and  x-y  planes,  respectively. 

n  '  =  n  and  n  1  =  n,  respectively.  The  eigenvalues  are  the  same  as 
y  z 

equation  (253)  replacing  Ex  with  E^  or  Ez,  respectively.  New  indices 
are  the  same  as  equation  (254)  replacing  Ex  with  E^  or  Ez  and  n  '  with 
nx*  or  n  '  with  n  ' ,  respectively. 
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Appendix  B:  Recommendations  for  Further  Research 


Self-Pumped  Phase  Conjugate  Mirrors  (SPPCMs) 

Papers  on  the  theory  and  experimentally  observed  phenomena  of  SPPCMs 
include  Cronin-Golomb  et  al.  (14:689-691),  Feinberg  (15:486-488), 
MacDonald  and  Feinberg  (16:548-553),  Gower  (21:196,197),  Rodriguez  et 
al.  (22:1732-11736),  Smout  et  al^  (23:77-82),  Gunter  (24:210-214),  and 
Eason  and  Smout  (25:51-53). 

Among  the  mechanisms  of  interest  in  this  effect  is  the 
self-defocusing  of  the  beam  or  '*beam  fanning"  (13:46-51)  by  which  part 
of  the  beam  is  redirected  as  it  passes  through  the  crystal  (Figure  28-a). 
The  reflection  back  to  the  original  beam  of  this  redirected  part  is  what 
sets  up  the  geometry  for  four  wave  mixing  (Figure  28-b). 


Figure  28.  Geometry  for  self-pumped  optical  phase  conjugation  showing 
(a)  beam  fanning  and  (b)  four-wave  mixing. 


Sinde  a  major  drawback  to  making  an  ordinary  phase  conjugate  mirror 
practical  is  the  required  coherence  of  the  probe  and  pumps,  the 
self-contained  geometry  of  a  self-pumped  phase  conjugate  mirror  appears 
to  be  the  solution. 
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Transient  Solutions  to  the  Photorefractive  Effect 


Feinberg  defines  the  grating  erasure  rate  as  follows:  If  an  index 
grating  is  written  in  the  crystal  by  the  interference  of  two  beams,  then 
these  two  writing  beams  are  turned  off  and  a  third  very  weak  beam  is 
scattered  off  the  grating  while  the  crystal  is  illuminated  with  strong 
uniform  light,  a  measure  of  the  time  decay  of  the  scattered  beam  is  a 
direct  measure  of  the  time  decay  of  the  amplitude  of  the  index  grating. 
The  modulation  index  under  uniform  illumination  is: 

m  =  0 

so  the  steady  state  portion  of  the  potential  is  zero  (5:1300). 

i  1 2 

Since  the  diffraction  efficiency,  n,  is  proportional  to  |w|  ,  the 

2 

decay  is  proportional  to  j w j  ,  as  well.  From  equation  (109): 

| w | 2  “  exp[-2F  (a2  +  1)  t]  (256) 

so  the  decay  rate,  A  ,  is  given  by: 

Ae  =  2  T  (a2  +  1)  (257) 

(F  as  defined  in  Chapter  IV.)  (5:1300) 

Gunter  presents  a  set  of  graphs  for  diffraction  efficiency  and 
intensity  gain  for  two  beam  coupling  versus  time  for  fixed  parameters, 
incident  beam  intensity  ratio  (his  Pq),  applied  field  and  fringe  spacing 
(9:262). 
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Figure  29.  Time  dependence  of  a.)  diffraction  efficiency,  n,  and 
b. )  two  beam  coupling  gain,  for  BSO  (9:262) 


How  Much  Does  Optical  Activity  Effect  BSO  Analyses? 


Optical  activity  is  the  rotation  of  the  plane  of  polarization  of 
linearly  polarized  light  as  it  passes  through  a  medium.  Figure  30 
represents  this  phenomenon.  As  light  propagating  in  z,  linearly 
polarized  in  x,  passes  through  the  medium,  its  polarization  vector  is 
rotated  until  it  is  linearly  polarized  in  y  at  some  length,  L  (2:94). 

Gunter  lists  the  optical  activity  of  BSO  to  be  45° /mm  for  514. 5nm 
wavelength  light  (9:240).  Another  existing  paper  on  optical  activity  is 
Simon  and  Bloembergen  (26:1104-1114). 
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Figure  30.  Illustration  of  optical  activity  (2:94) 
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